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Ab initio quantum mechanical calculations at the self 
consistent field level have been performed for interacting 
water molecules to predict their infrared spectra. These 
calculations have been made using Gaussian 76 and a 4-31G 
basis set. The calculated ab initio force field has been 
used to predict the fundamental absorption frequencies and 
the corresponding normal coordinate transformation. This 
normal coordinate transformation has then been used in con- 
junction with Atomic Polar Tensors (APTs) that were also 
found in the ab initio calculation to predict the absolute 
integrated intensities for these fundamental normal 
vibrations. 

The role that the chemical environment plays in the 


predicted spectra for interacting water molecules has been 


investigated by analyzing the predicted spectra for different 
assumed orientations of the "complex" at various assumed. inter- 
molecular distances, and by varying the electron donating 
properties of the electron donor. The former effect was in- 
vestigated by predicting spectra for water pairs in the linear, 
cyclic, and bifurcated configurations, at varying intermolecular 
distances. The latter effect was studied by calculating the 
spectra for water complexed with two oxygen atoms which have 
varying charge co! . H50 . of where M = 0, -l, or -2) where 
the oxygen atoms with the larger negative charge are stronger 
electron donors. Finally, the spectrum for ionized water, 
HAO" is also predicted. 

These calculations indicate that both the frequencies 
and intensities of the vibrations for water interacting with 
its environment are predicted to be strongly dependent upon 
the nature of this environment. The predicted spectra for 
these water molecules in different chemical environments are 
then compared with the predicted spectrum for isolated mono- 
meric water, and also experimental spectra when available. 
The vibrational changes from isolated water to the interact- 
ing water molecule in these different environments are then 
analyzed and discussed to provide some insight into the effect 


of intermolecular interactions on the vibrational spectrum of 


the interacting molecule. 


CHAPTER 1 
INTRODUCTION 

The infrared spectrum of a molecule is an extremely rich 
and valuable source of information. For many years it has 
been used in both the qualitative and quantitative analysis 
of chemical compounds. The concept of group frequencies (1), 
and its use in the identification of molecular compounds is 
one of the most fundamental and fruitful applications of 
infrared spectroscopy. As a result, most experimental spec- 
troscopists, and theoreticians working with them, tradi- 
tionally placed great emphasis on the locations of the 
band centers, or absorption frequencies. In recent years, 
however, there has also been an increased interest not only 
on the fundamental infrared absorption frequencies, but also 
on the absolute integrated intensities (2,3). Specifically, 
there has been a great deal of interest in the parameteriza- 
tion and interpretation of infrared intensities. The use of 
electro-optical intensity parameters (EOPs) (4,5) has been 
pioneered by the Soviet workers (6), while the use of Atomic 
Polar Tensors (APTs) was first introduced by Biarge, Herranz, 
and Morcillo (7), and later reformulated by Person and Newton 
(8) and used in conjunction with quantum mechanical calcula- 
tions of signs of (3P/3Q)s. Both methods have had good success 


in the parameterization of infrared intensities (9). The APT 


method is particularly useful in quantum mechanical calculations 
and therefore has been adopted here. 

The success of the quantum mechanical calculations of APTs 
offers hope that in the future the prediction of infrared in- 
tensity parameters may be almost routine. Using these inten- 
sity parameters, the absolute integrated intensities can be 
predicted, provided that the normal coordinate transformation 
is well known. The normal coordinate transformation can be 
easily obtained from the force field of the molecule by per- 
forming a normal coordinate analysis (10,11). Unfortunately, 
the available experimental data are usually not sufficient to 
determine completely an accurate vibrational force field. . For 
this reason quantum mechanical calculations of vibrational 
force fields have become increasingly more popular (12-16). 
These calculations, coupled with the calculation of APTs, have 
made possible the completely ab initio prediction of infrared 
spectra (including both band frequencies and absolute 
intensities). 

These predictions are useful for several reasons. They 
can be extremely helpful in assigning complicated spectra 
(17), especially spectra for relatively large (> 10 atoms) 
molecules. A second application of these types of calcula- 
tions concern chemical systems which cannot be isolated in 
pure form and packaged in a bottle. These systems include 
radicals and ions (18) frozen in rare-gas matrices, reaction 
intermediates observed by time-resolved spectroscopy or 


observed on catalytic surfaces, molecular fragments existing 


in flames and plasma, small molecular compounds that are be- 
lieved to exist in space, and dimers or other intermolecular 
complexes that are found in rare-gas matrices or supersonic 
expansions. In recent years there has been an increased 
interest in the application of vibrational spectroscopy 
towards characterizing these unstable and hard to isolate 
molecular species. However, this task is made much 
more difficult by the absence of spectroscopic standards 
and standard spectra. Therefore, any spectral determination 
of the concentration of these transient species must rely on 
a separate estimate of its intensities. The development of 
ab initio calculations of these intensities appears to be 
well suited for making these estimates. Additionally, 
the spectral patterns predicted by these calculations (includ- 
ing both wavenumbers and intensities) can be qualitatively 
used as a direct confirmation of the spectral assignment for 
the species, especially in the absence of standard spectra. 
The calculated intensities would be especially valuable in 
the case that the observed wavenumbers of the unstable species 
could not be compared with those of a structurally similar 
stable molecule; or in the opposite situation where the ob- 
served wavenumbers matched equally well with more than one 
structure. In situations such as these the correct assign- 
ment must match the predicted wavenumbers as well as the 
predicted intensity pattern. 

There are very encouraging developments in the ab initio 


predictions of infrared spectra of unstable or hard-to-isolate 


molecular species. However, these methods must be further 
tested if one is to be able to routinely predict reliable 
spectra for these types of species. Therefore, we have used 
the ab initio method to predict and interpret spectra of 
some of these types of species. 

The water dimer has been the subject of many experimental 
(19-29), and theoretical (30-41) investigations. One reason 
for this interest is that it is thought to be a prototype 
model for liquid water, which is of paramount importance in 
the life of biological systems. A second reason is that the 
water dimer is a classic example of a donor-acceptor complex, 
in this case, specifically a hydrogen bonded complex. The 
water dimer represents a unique example of a hydrogen bonded 
system because one water molecule is the electron donor and 
another is the electron acceptor. 

The importance of vibrational spectroscopy in the study 
of H-bonded systems is well documented (42,43). The infrared 
spectra of such complexed species exhibit characteristic 
spectral changes when compared with that of the corresponding 
isolated molecules and with the spectra of other complexes. 
The comparisons among these spectra are valuable sources of 
information concerning the properties and characteristics 
of these systems. 

The hydrogen bond interaction can be described as an 
intermolecular interaction through a hydrogen atom. Typi- 
cally, the interaction is of the type RAH--B, where RAH is 
an electron acceptor, B is an electron donor molecule, and 


the hydrogen bond is the H--B interaction. When water is 
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both the electron acceptor and the donor, A-H is the O-H bond 
of one water molecule, and the electron donor B is the oxygen 
atom in the other water molecule. In the water dimer, this 
interaction is slightly more complicated by the presence of 

two O-H bonds. In the isolated molecule, these two bonds are 
symmetrically equivalent, and therefore are strongly coupled 

to form two stretching vibrations, a symmetric stretch and an 
asymmetric stretch. In the complex the symmetry is reduced in 
relation to the isolated monomer, the O-H bonds are no longer 
symmetric, and the coupling between these two bonds is reduced. 
In general, the forces governing intramolecular coupling are 
much stronger than those governing hydrogen bonding. Therefore, 
this reduction in coupling between the two 0-H bonds is important 
to consider in analyzing complexes involving water molecules. 

For H-bonded complexes, the vibration most perturbed by 
the interaction is the symmetric stretch. This band undergoes 
a shift to lower energy, and an increase in intensity. The 
Similarities and differences between these two changes in the 
spectra of H-bonded systems are vital to our understanding of 
these complexes. 

The frequency shift is usually thought to originate ina 
weakening of the A-H bond as the H-bond is formed. This in 
turn is believed to be due primarily to two effects: the 
polarization effect, where the donor molecular polarizes elec- 
tron density from the A-H bond into the non-bonding region of 
the A atom, where it no longer compensates for the A-H nuclear 
repulsion; and the electron transfer.effect, where the donor 


molecule donates electron density into the antibonding 


orbitals of the acceptor molecules, thereby weakening the bond. 
Regardless, the net result is that the A-H bond is weakened, 
and so the frequency shifts towards lower energy. This effect 
can be analyzed by calculating the internal force constants 

for the complex. 

The intensity enhancement-is well known experimentally, 
but not as well understood. The electrostatic polarization 
and charge-transfer effects just mentioned cannot individually 
account for the great enhancement in the intensity. Zilles 
and Person (39) have recently applied a quantum mechanical 
model of the APTs in order to interpret the intensity enhance- 
ment for the water dimer. They have found that the largest 
single most important contribution to the intensity enhance- 
ment is from the increase in charge flows in the complex 
compared to the monomer. This interpretation can be taken 
one step further by attempting to compare the nature of these 
charge flows in the linear dimer with those in other systems 
involving water as the electron acceptor. 

Furthermore, the reduction in coupling between the two 
O-H bonds usually creates some mixing between the symmetric 
stretching vibrations and both the bending and assymetric 
stretching vibrations. So, in effect, it is no longer strictly 
a symmetric stretching vibration, but a mixture of all three 
vibrations. Naturally, a consequence of this mixing is that the 
normal vibrations in the complex are no longer the same vibra- 
tions as those in the isolated molecule. The intensities, as 


well as the frequencies, are strongly dependent upon the form 


of the normal vibrations. The effect of these changes on both 
the frequencies and intensities are also important to 
understanding the spectral changes present in the complex, 

and these effects will be examined in detail. 

The experimental studies of the infrared spectrum of the 
water dimer are quite varied with respect to the interpretation 
of the results. Most of these results have been previously 
summarized (41 and references citedtherein). The most gener- 
ally accepted interpretation of the infrared spectrum is that 
of Tursi and Nixon (19). They have assigned the six high fre- 
quency modes of the water dimer in a nitrogen matrix. Their 
assignment is consistent with a linear structure (see Fig. 
1.1) for the dimer, and is accepted by many as the correct 
interpretation for the dimer. Other workers prefer the cyclic 
(see Fig. 1.2) structure (20,21,24), while still otherspostu- 
late the existence of the bifurcated structure (see Fig. 1.3) 
(26). 

Most theoretical calculations support the linear struc- 
ture as the most stable geometry for the water dimer (31,34). 
Zilles and Person (39) have used the experimental force field 
of Tursi and Nixon (19) for this linear structure and quantum 
mechanical APTs to predict and interpret the infrared intensi- 
ties of the linear dimer. To expand and further develop our 
understanding of the infrared spectra of interacting molecules 
we have used ab initio quantum mechanical calculations to pre- 
dict the vibrational force fields and APTs for not only the 


linear configuration, but also the cyclic and the bifurcated 
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structures. We will examine the spectra predicted for these 
three different configurations and discuss the differences 
among them, as well as the differences between our calcula- 
tions and those of Zilles and Person (39). In addition, we 
have tried to simulate the extremum, namely complete electron 
donation, by calculating the spectrum of H,0°. 
Finally, the infrared spectrum of a model for water iso- 
lated in a beta-alumina complex has been calculated. Water, 
as an impurity in beta-alumina, is important to the electrical 
conductivity and mechanical strength properties of beta- 
alumina (44). This in turn is important towards the use of 
beta-alumina as solid electrolyte in the sodium-sulfur battery. 
Until recently, the nature of the hydration reaction of beta- 
alumina was not known. Bates and coworkers (45) have deter- 
_mined that the water molecules diffuse into the conduction 
layers and form strong electrostatic bonds with the mobile 
cations, thereby affecting the properties of the beta-alumina. 
Infrared spectroscopy is used as a tool to probe the kinetics 
and thermodynamics of the reaction of water with beta-alumina. 
Therefore, before attempting to use the infrared spectrum as 
a probe for other properties, it is important to first under- 
stand the structure and vibrational properties of hydrated 
beta-alumina. Bates and coworkers have recently determined 
the crystal structure and vibrational spectrum of hydrated 
beta-alumina (44). The water molecules are located in the 
crystal so that they are strongly hydrogen-bonded to the 


oxygen atoms of the aluminate ion (see Fig. 1.4). 
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Consequently, this is a good system to use to expand our 
studies on hydrogen bonded water molecules. 

The quantum mechanical calculation of the spectrum of 
hydrated beta-alumina takes on importance for several reasons. 
First, it offers the opportunity to verify (or dispute) the 
experimental spectrum and the interpretation of this spectrum. 
Secondly, it represents a different example of an intermolecular 
interaction involving a water molecule. It offers the oppor- 
tunity to investigate the role that the electron donor plays 
in intermolecular interactions. And lastly, it is a challenge 
to our ability to simulate a complicated chemical system by 
using a rather simple model. 

A summary of the calculational procedure developed for 
this work, along with a review of both normal coordinate 
analysis and infrared intensity analysis,is given in Chapter 
2. The results of these calculations on isolated water mono- 
mers and the linear, cyclic, and bifurcated forms of a pair 
of interacting water molecules are presented in Chapters 3, 

4, 5, and 6 respectively. The importance of some of the 
effects just discussed, such as coupling between O-H bonds, 
weakening of the O-H bond, charges flows, and the mixing of 
normal coordinates, are discussed with respect to each of the 
systems calculated. The predicted spectrum for our model of 
water isolated in a beta-alumina complex is given in Chapter 
7, and the predicted spectrum for H,0" is presented in 
Chapter 8. A final summary and suggestions for future cal- 


culations are presented in Chapter 9. 


CHAPTER 2 

THEORY 
The theory of both normal coordinate analysis (10,11, 46) 
and infrared intensity analysis (2,3,41,47,48) is well known 
and described fully elsewhere. We have developed a slightly 
different approach to analyze simultaneously the normal coor- 

dinated and infrared intensities. This method will be dis- 

cussed in detail, including the similarities and differences 


with the more traditional procedures. 


2.1 Coordinate Systems 


A polyatomic molecule of N atoms has 3N degrees of 
freedom. Of these 3N degrees of freedom some are rigid body 
motions (translations and rotations) and some are elastic 
body motions (vibrations). The analysis of the molecular 
vibrations is best done in a coordinate system which allows 
the maximum separation of the overall motions into transla- 
tions, rotations, and vibrations. This problem has been dis- 
cussed in the past (11,49,50), and amounts to choosing a 
coordinate system such that the interactions between vibra- 
tions and the translations and rotations are minimized. The 
necessary coordinate system is a molecule fixed coordinate 
system which translates and rotates along with the molecule. 


These conditions require that there is no component of either 
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BA 


translation or rotation present in a vibration, are are known 
as the Eckart-Sayvetz conditions (1l). 

The Eckart-Sayvetz conditions provide explicit defini- 
tions for the translations and rotations of a molecule. They 
are given below in terms of mass-weighted principal cartesian 


coordinates (11). 


= | 2a 
are qy 
o* «a a 
L 45 
0, = mS =m q 
- Mo a r Yq 
Be 4s 
Po. = prim q 
3 M,% g a xy 
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4 a xX i pk ee 
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Here Par Por P3r Par Poy and Pe represent translation and 
rotation in the x, y, and z directions, respectively, Mo is 


the total mass of the molecule, mo is the mass of atom , 


° 12} 
x 


Ps oe and z. are the principal cartesian coordinates for 


atom a, q, , 4, , andq, are the mass weighted principal 
a a 


Yo 


cartesian coordinates for atom a, and I. Toy’ and a are 


x Zz 


the moments of inertia along the principal cartesian axes, 


or the principal moments of inertia. 
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The remaining vibrational coordinates should reflect these 
constraints, and at the same time form a complete, nonredundant, 
vector representation of the vibrational subspace. The calcu- 
lations can be further simplified by utilizing the symmetry of 
the molecule and by requiring that all the coordinates be 
orthogonal to each other. The resulting coordinate system 
is a mass-weighted, molecule-fixed, symmetry-adapted, prin- 
cipal cartesian coordinate system, hereafter denoted as the 
S$ coordinates. 

The transformation from mass-weighted principal cartesian 
coordinates q to S$ coordinates can be written in matrix 


notation as 


S=Uq (2a) 


where U is an orthogonal transformation matrix and q is the 
column vector representation of the mass-weighted principal 
cartesian coordinates. The matrix U can be constructed via 
a Gram-Schmidt orthogonalization (51) of the irreducible 
representations of the cartesian coordinates to the Eckart- 
Sayvetz conditions. The linear combinations of the principal 
cartesian coordinates which reduce as the irreducible repre- 
sentations of the point group of the molecule can be found 
by following the prescription given in Cotton (52). These 
linear combinations are then made orthogonal to the Eckart- 
Sayvetz conditions, giving the U matrix. 

Since U is constructed as an orthogonal matrix, its in- 


verse is simply u* where u* represents the transpose of ft. 
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Equation 2a can now be rewritten as 

q=u's (2b) 
where all the symbols have their previously defined meanings. 
Equation 2b can now be applied in the method of finite differ- 
ences (53) in order to calculate the force constants and dipole 
derivatives in S$ coordinates. The values of 5 are incremented 
by a specific amount, usually +0.01 u? A°, and equation 2b is 
used to calculate the values at q of these displaced values 
of 5. 

The Gaussian 76 ab initio molecular orbital program (54) 
is then used to calculate the energies and dipole moments for 
the equilibrium geometry and also at appropriate displacements 
from equilibrium. The force constants are derived by assuming 
a pure quadratic potential, and then calculating the second 
derivatives of this potential function. Mathematically, these 


constants can be calculated by using the following equations 


(from (53)) 
E(x_+Ax) + E(x_-Ax) - 2E(x_) 
F ROE a < Ma, Sa a RS Lear i * (3) 
xx 2 
Ax 
- . (E(x, +4x,y,tdy) = E(x,+Ax,y,) - E(x 5 +Y,tAy) + E(x,+Y,) 
xy AxAy 
(4) 


where —_ represents a diagonal force constant, P ey an off- 
diagonal force constant, E(x.) the energy at the equilibrium 
configuration, E(x +4x) and E (x,-Ax) the energies at the pos- 


itive and negative displacements of coordinate x from 
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equilibrium respectively, E(x, +Ax,y +Ay) the energy at the 
simultaneous positive displacement of both coordinate x and 

y, and Ax the numerical value of the displacement of coor- 
dinate x. The dipole derivatives are approximated by averag- 
ing the displacements at the positive and negative distortions, 


to account in part for the nonlinearity of the dipole moment. 


P os = a 
i 3 AP, e P(X, tAx) Ry Be) ‘ P(X Ax) P(X) 
ox Ax 2Ax 2Ax 
(5) 
= Py (XQtAx) -P, (x -Ax) 
2Ax 


Here P(X) represents the permanent equilibrium dipole moment 
in the a = x, y, or z direction, Pi (xp Ax) and P , (xy Ax) repre- 
sent these same components of the dipole moment at the positive 
and negative displacements of coordinate x respectively, and 


Ax is the displacement of coordinate x. 


2.2 Normal Coordinate Analysis 


After the force constants have been calculated, the next 
step is to derive the vibrational frequencies and the nature 
of the normal vibrations (or the normal coordinate 
transformation). For a harmonic oscillator the solution of 


Newton's equations of motion gives the secular equation (11), 
[GF - EA| = 0 (6) 


where F is the potential energy matrix, G is the inverse kin- 
etic energy matrix, E is the identity matrix and A is the 
eigenvalue matrix. The solution to this system of homogen- 


eous linear equations gives the eigenvalues, or normal 
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frequency paramaters A, and the eigenvectors, or column vector 
representations of the normal coordinates, also known as the 
normal coordinate transformation matrix, L. 

If the potential energy matrix is represented in mass- 
weighted cartesian coordinates, then the G matrix becomes 


the identity matrix (11). Equation 6 can be rewritten as 


RKimwLA (7) 


where K represents the matrix of force constants in mass- 
weighted cartesian coordinates, £ is the eigenvector matrix 
in these coordinates, and A is the diagonal eigenvalue matrix. 


The values of A are related to the vibrational frequencies by 


4408 md 1 = 4 
rl ) =; Rey * (5==) or v, = 1303.1*,7 (8) 


and the matrix 2 is the orthogonal normal coordinate trans- 


formation which transforms from normal coordinates to mass- 


weighted cartesian coordinates 

q=24Q (9) 
where Q is the vector of normal coordinates. Thus in mass- 
weighted cartesian coordinates solving the secular equation 


reduces to finding the matrix 2% which diagonalizes the 


potential energy matrix, K. 


27K h=A (10) 


Therefore the force constants calculated in S$ coordinates must 


24 


be transformed to force constants in mass-weighted cartesian 
coordinates K. The force constants in S$ coordinates, K, can 
be easily transformed to mass-weighted cartesians using the 
U matrix. Since this matrix is orthogonal its inverse is 
simply equal to us In mass-weighted cartesians the poten- 
tial and kinetic energies are equal to 


ar = qo 1g (11) 


2v=q°kq (12) 


respectively, where g represents differentiation with respect 
to time, and e* is actually the identity matrix, E, since 
the potential energy is now expressed in mass-weighted 
cartesian coordinates. 

In S$ coordinates the potential and kinetic energy can 


be written as 


2art=S G8 (13) 


2yv=e5°K«S (14) 


where G is the inverse kinetic energy matrix in $ coordin- 


ates and K is the potential energy matrix in S. Since 


S=uqandS* =q*u* (15) 
equation 13 and 14 can be rewritten as 
T+, + -1 = 
erm @ U'G" Ug (16) 
av=q*ukKug (17) 


25 


By comparing equations 12 and 17 the transformation from K 


to K can be written as 
k=u"Ku (18) 


The secular equation can now be solved in terms of mass- 
weighted cartesian coordinates using equation 10. 

Here a is a matrix of dimension 3N by 3N, and each row 
of this matrix represents one of the 3N normal modes. Six (or 
five for linear molecules) of these modes are the translations 
and rotations, and the remaining modes are the fundamental 
vibrations. Each row of a consists of N (3 by 1) vectors 
corresponding to the motion of each atom N in that particular 
normal coordinate. By analyzing the corresponding vectors for 
each atom we can determine what motions constitute that par- 
ticular normal vibration. This will be discussed in more 
detail later. The matrix A is a diagonal one containing 
the eigenvalues of the secular equation. Six of these eigen- 
values correspond to the translations and rotations of the 
molecule, and are therefore of value zero. The remaining 


values give the vibrational frequencies via equation 8. 


2.3 Infrared Intensity Analysis 


The absolute infrared intensity of the ith vibrational 
mode of a molecular is proportional to the square of the 
dipole moment derivative with respect to motion along the 


ith normal coordinate. In units of km/mole this becomes (2) 


A, (km/mole) = 974.9 4, | (aP/2Q,) | (19) 


26 


where dq, is the degeneracy of the ith mode, and (aP/aQ, ) is 
given in units of e-u? where e is the atomic unit of charge, 
equal to 1.602 x 10719, If the value of the dipole moment 
derivative with respect to the normal coordinate is known, 
then the intensity for that mode can be easily calculated 
with equation 19. The problem of infrared intensity analysis 
then reduces to evaluating and interpreting this derivative. 
As stated previously, the method of finite differences 
can be used to approximate the dipole moment derivatives. By 
using equation 5, the dipole moment derivatives in S$ coordin- 
ates can be calculated. Since the normal coordinate transfor- 
mation is given in terms of mass-weighted cartesians, then the 
dipole moment derivatives must be transformed to this same 
coordinate system in order to utilize this normal coordinate 
transformation to calculate the intensities. This transfor- 


mation is written as 


Fq ™ Pail a 


where Pa is the 3x3N matrix of the dipole derivatives in mass- 
weighted cartesians and Po is the corresponding 3x3N matrix 

of dipole derivatives in §. The first six (or five for linear 
molecules) columns represent the dipole derivatives with re- 
spect to the three translational and three rotational motions 
respectively. These values can be calculated from the charge 
of the molecule (qa)e the mass of molecule (Mo). the principal 
moments of inertia (Ty? Tog ang I. 


yy 
ents of the permanent dipole moment re, Py! and Po 


st and the vector compon- 


These relations are 
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given by King (55) and have recently been reformulated by 


Rogers and Hillman (56), and are summarized below in columns 
one through six. The remaining columns of this matrix, Por 
are calculated using equation 5, and represent the Aipele> 
derivatives with respect to the vibrational 5 coordinates. 
The three rowsof this matrix represent differentiation of 
one of each of the three vector components of the dipole 


moment; either Poe Py! or P 


qe - -§ AP, AP 
7 fe) fe) O —_ —% AS, tnt oe 
yy ZZ 
q =p? p? AP AP 
Po=| 0 = O —r ) —; ASe cee = 
a My ‘i r 7 3N 
xX VAY A 
(eo) [e) 
a3 P -PY AP AP, 
° ee a lad, aS 
M r I 7 3N 
° xX yy 
(29 


The matrix Pa consists of the dipole moment derivatives 
in mass-weighted cartesian coordinate space. By removing the 
mass dependence the derivatives in cartesian coordinates are 


obtained 


Pp =p mM? (22) 
2G. 


where am is composed of the juxtaposition of the Atomic Polar 
Tensors (APTs) (57) for the N atoms in the molecule, and mi 
is a diagonal matrix represented by the triplet sets of the 


square root of the masses of the atoms in the molecule. The 
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By matrix has the general form (3 x 3N) 


dP dP oP oP dP dP 
aa ae Hh ae x_ rz 
aX, oy dZy xy dYn au 
oP dP dP dP oP dP 
Pp = X = = Ce aa Re = 
Py 2 . 3X La ae 
oP. oP. on oP. oP. oP. 
OX, oy. 02, aX OX dZy 
(23) 


The APTs will be discussed in more detail later. 


The transformation from Pa tO P and in essence the abso- 


Q’ 
lute integrated intensities, can be accomplished with the L 
matrix 


PSP 2 (24) 
pee i 


where Py is the matrix of dipole moment derivatives in normal 


coordinate space. This matrix is also 3 by 3N, with the first 
five or six columns corresponding to the dipole derivatives 
with respect to the translations and rotations. The values of 
these derivatives are all zero since there is no component of 
vibration in a pure translation or Estetica and hence no in- 
tensity or frequency for these modes. Each of the next (3N-6) 
columns correspondsto the dipole derivatives with respect to 
one of the (3N-6) vibrational normal modes. The derivative of 
one of the three vector components of the dipole moment Pie 

Pj; and oe are again each responsible for one of the three rows 


Y 
of this matrix. This matrix has the following representation 
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aP oP 
0 0 0 0 0 G ‘eee e. 
/ 9Q5 9Q3n 
9P 9P 
Pp. =| 0 0 0 0 0 © vee capes aah eee 
Fo 30, 305, | (25) 
ae dP 
0 0 0 0 0 Onc see 


aQ- ol ees 


The absolute integrated intensity for the ith mode is 
simply proportional to the sum of the squares of the vector 
components of the dipole moment derivatives for the ith 
column. The actual value of the absolute intensity in km/mole 
is found using equation 19, where the following relationship 


is now true 


2 2 2 
2 aP aP aP 
iz | al eaten (26) 
8Q5 9Q5 9Q; 90; 


2.4 Transformation to Different 


Coordinate Systems 


The § coordinates are very convenient for the calculations 
of frequencies and intensities. However, the detailed analyses 
of the frequencies and intensities are best accomplished in two 
different coordinate systems. The frequency parameters are 
most easily analyzed in terms of force constants in internal 
coordinates and the intensity parameters will be discussed 
in terms of dipole derivatives with respect to space fixed 
cartesian coordinates, or APTs. Therefore, the frequency and 
intensity parameters in S$ coordinates must be transformed to 
the corresponding parameters in internal and cartesian coor- 


dinates respectively. 
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2.4.1 Internal Force Constants 

As chemists, we have a desire to interpret the absorption 
frequencies in terms of molecular structure. Therefore, the 
necessary frequency parameters should be related to the molecu- 
lar geometry. The most useful frequency parameters are thus 
the force constants in internal coordinates (10,11). 

The internal displacement coordinates are changes in 
bond lengths and bond angles of a molecule. For most appli- 
cations five coordinates -- bond stretching, angle bending, 
torsional motion, out-of-plane bending (wagging), and linear 
bending-—-are sufficient to describe completely all of the vi- 
brational motions (10,11). For more unusual situations other 
specific internal coordinates can be defined by following the 
prescription of Califano (11). More detailed information on 
internal coordinates can be found in references 10 and ll. 

The internal coordinates are related to the cartesian 
coordinates by 


he 


k = 1, 2 ««« (3N=-6) 
(27a) 


Ry = 


By wheres o oy1)343%.. (38) 


i 


or in matrix notation 


R=BX (27b) 


where B is a (3N-6) by 3N matrix with typical element defined 
as (OR, / 3%; ) and R is the column vector representation of 
internal coordinates. The B matrix is a function solely of 
the molecular geometry, and formulas for the calculation of 


its matrix elements for the five most common type of internal 
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coordinates have been derived (1l). The A matrix transforms 


the internal coordinates back to cartesian coordinates 
X=AR (28) 


where A is apparently the inverse of B. However, since the 
B matrix is not a square matrix, but rather a rectangular 
one, it does not have an inverse, and the A matrix is calcu- 
lated via the G matrix (58). 

In internal coordinates 


@ =p t-g* (29) 


where ut is a (3N by 3N) diagonal matrix of the reciprocal 


masses. The "inverse" of B is then obtained by (47) 


B* ¢=B pm? pt (30) 
ge get gt et (ai3 
A=M Bt gt (32) 


At this point is it appropriate to mention that A in equation 
32 can be calculated only if the G matrix is a square non- 
Singular matrix. This in turn means that there cannot be 
any redundancies defined in the B matrix elements. This is 
often a problem if one wishes to fully utilize the symmetry 
of the molecule. The methods of dealing with this problem 
will be discussed later; suffice it to say there that the 
problem exists and must be resolved. 

The relationship between the force constants in mass- 


weighted cartesian coordinates K and in internal coordinates 
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Fs is now given by 


FL = AT KA (33) 


+ , . 
where A is the transpose of the A matrix, and F; is the force 


constant matrix in internal coordinates. 


2.4.2 Atomic Polar Tensors 

The APT has been defined in the first three columns of 
equation 23 above, and equation 20 calculates the APTs in 
principal cartesian coordinates. However, usually the prin- 
cipal cartesian coordinate system is not ideally oriented 
for the interpretation of infrared intensities. This arises 
in part because the actual values of the polar tensor elements 
depend upon the orientation of the axis systemused. Therefore, 
the bond coordinate system is usually the best system for ana- 
lyzing the APTs. The transformation of the APT from one coor- 
dinate system to any other system is done by rotating the axis 
system. As the axis system is rotated, the values of the polar 
tensors change according to the rotation performed. The new 
values of these tensors can be calculated by performing the 
appropriate tensor rotations on these tensors. The matrices 
for the rotation of the axis system by the angle 6 about 


either the X, Y, and Z axes respectively are (from (59)) 


1 0 0) cos8 0 sind cos@ -sin@ 0 \ 
R. ={0 cos@ -siné R = 0 a 0 Rz =| siné cosé 0 | (34) 
i 


0 sin cosé -sind cos 0 0 i 
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where the sign of the angle 6 is determined by viewing along 

the rotation axis from the positive direction. If the old 

axis system must be rotated in a clockwise direction to 
transform to the new coordinate system, then the angle is 
positive, and if the old axis system must be rotated in a 
counter-clockwise direction to transform to the new coordin- 

ate system, then the angle is negative. The rotation matrices 
defined in equation 34 are the ones usually used for the rotation 
of a cartesian axis system (59). The transformation for tensors 
is given by 


' 
+ 
T = Ry TR, Q,=tX, Y,. OL Z (35) 


' 
where T is in thenew coordinate system, and T is the tensor in 


the old coordinate system. 


2.4.3 The S$ Coordinates 

We have just seen the advantages of interpreting the pre- 
dicted frequencies in terms of internal force constants and the 
predicted intensities in terms of APTs. Why, you may ask, 
have we chosen to perform the actual calculations in 5 
coordinates? This question is discussed below. 

The transformation from cartesian coordinates to S coordin- 
ates is performed by the orthogonal matrix U (see equation 2a). 
The orthogonality of U makes the reverse transformation matrix, 
from S$ coordinates to q coordinates, simply the transpose of 
U, u* (see equation 2b). In internal coordinates the analo- 
gous equation to equation 2a is equation 27. The reverse 


transformation is then performed by equation 28. However, 
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since the B matrix in equation 27 is not square, the A matrix 
in equation 28 cannot be obtained by simple inversion of B, 
but must be obtained using equation 32. Equation 32 requires 
that all of the redundancies which are often necessary to com- 
pletely describe the vibrational degrees of freedom using only 
the previously defined five internal coordinates be removed. 
The method of finite differences requires the transformation 
given by either equation 2b or equation 28. In internal coor- 
dinates equation 28 is often not easily obtained because the 
redundancy cannot be removed in any unique manner. Therefore, 
we avoid the problem with redundant coordinates by using 
equation 2b, or the S$ coordinates. 

The cartesian coordinates used in the definition of the 
APT also avoid the redundancy problem. However, these coor- 
dinates do not utilize the symmetry of the molecule, nor do 
they explicitly separate the vibrational degrees of freedom 
from the translations and rotations. This means that in 
cartesian coordinates there are 3N coordinates rather than 
3N-6, and the symmetry of the molecule cannot be used to 
simplify the problem. This greatly increases the computa- 
tional effort necessary to predict the infrared spectrum. 
Therefore, since the S$ coordinates explicitly separate 
the vibrations from the translations and rotations, and 
utilize the symmetry of the molecules, they are 


preferred to cartesian coordinates. 
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2.5 Experimental Infrared Intensities 


The experimental absolute infrared band intensity is 


given by (60) 


A, (km/mole) = (aS, In(I,/1) dv, (36) 


‘ 


where A, is the absolute integrated molar absorption coeffi- 
cient or intensity for the ith fundamental mode in k/mole, 
c is the sample concentration in moles/liter, 2 is the path- 
length in cm, I, is the intensity of the incident light in 
the absence of the sample, I is the intensity of the trans- 
mitted light after passing through the sample, V5 is the 
wavenumber of the band in cmt, and the integration is per- 
formed over the entire band region. The factor of 100 con- 
verts A; into the units of km/mole. In these units a very 
strong band would have an intensity of around 1000 km/mole, 
while a very weak band would have a value of 0.10 km/mole (57). 
The alternate definition of the experimental integrated 
absolute band intensity uses integration over inv; instead 
of fT, (60). This quantity is given as the symbol A and is 


1 


defined as 


2 big 
rT, (cm /mole) = aS, in(I./T) d(1nv,) (3:7) 


where all the symbols have their previously defined meanings 
and units. The definition of rT; is preferred for many purposes 
(61,62), but most spectrometers record spectra linear in wave- 


numbers, so integration over VG is much easier to measure 
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experimentally. As a result, A; is the more commonly reported 
value. The difference between A, and qr. is not too severe, 


however, the two quantities being related through 
A. = T. v; (38) 


Here Vi is the band center in cms; A; is the intensity 
in km/mole and ry Measures the intensity in em*/mole. 

The quantities A; measure the total vibrational strength 
for a transition from the ground vibrational state to an 
excited vibrational state (v5 =0 > ee 1) (60). Theo- 
retically, the most important quantity related to this 


strength is the transition dipole moment Minn §+) (from (57) ) 
* “aw 
Mion (+) =f ¥. fF vy dt (39) 


where the transition occurs from state m to state n, P is the 
dipole moment operator, o is the wavefunction for vibrationally 
excited state m, and ve is the wavefunction for the initial 


state n. 


2.6 Theoretical Infrared Intensities 
The absolute intensity of an infrared absorption band is 


related to the theoretical quantity inn (2) by (60) 


or a 
81 er 2 


age Men ST og (40) 


where Ns is Avagadro's number, V5 is the wavenumber of the 
absorption band, h is Planck's constant, c is the speed of 


light, and a; is the degeneracy of the ith vibration. If the 


SH 


vibrations are harmonic so that Vin and Vn in equation 39 are 
harmonic oscillator wavefunctions, then equation 40 becomes 
(60) 


av™4; oP 2 

A, (km/mole) = a cach ee 30, (41) 
where (3P/3Q;) is the dipole moment charge with respect to 
motion along the ith normal coordinate. If the proper con- 
stants are substituted into equation 41, then that equation 


can be simplified (57) to 


A, (km/mole) = 974.9 |(3P/2Q,)|* di (42) 
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, where e is the atomic 
a 


where (3P/3Q,) has units of e-u_ 
unit of charge (1 e = 1.602 x 10° 9G) and u is the atomic 
Mass unit. The units are Sea for (aP/aQ;) when the value 

of P is expressed in e A (where 1 e A = 4.8 Debye) and the 
normal coordinates Q are in units of A ut, 

Equation 42 indicates that the intensity for the ith 
mode depends on the form of the ith normal coordinate. Un- 
fortunately, the normal coordinates are strongly dependent 
upon the molecule in question. They usually involve motion 
of all the atoms in the molecule simultaneously, and the 
amplitudes of motion depend also on the masses of the atoms. 
The interpretation and parameterization of infrared intensi- 
ties in terms of these quantities are not easily accomplished 
because of the dependence of the form of the normal coordin- 


ates on the molecule in question. This dependence makes it 


difficult to compare these derivatives (3P/3Q) among different 
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molecules. Instead, it appears to be more useful to attempt 
to use the fact that the atoms in a molecule are connected by 
bonds, and that these atoms (in similar chemical environments) 
should exhibit similar properties within different molecules. 
The most natural coordinate system to use to try to inter- 
pret the absolute intensities are the cartesian coordinates. 
These coordinates are uniquely defined for every atom in the 
molecule, and are not dependent upon the geometry of the 
molecule. They can be oriented in the same direction in any 
molecule, thereby offering the possibility of directly compar- 
ing parameters among different molecules. This is necessary 
if we hope to be able to interpret intensities in terms of 
chemical structure. The transformation from normal coordin- 
ates to cartesian coordinates can be derived using the chain 


rule of differentiation 


a ash 
ax; ig | 2Q5 | | ax, 

3P aP + 

sor Sul [ S—] -G. . (44) 
an, i 9Q, 13 


where (3P/3Q; ) are derivatives of the dipole moment with re- 
spect to normal coordinates Qi, (aP/9Q,) are these same deriva- 
tives with respect to cartesian coordinates 5. and (80; /2X.) 
is the derivative of the normal coordinate with respect to a 
cartesian coordinate. This last derivative is one of the 
elements in the normal coordinate transformation matrix. This 
matrix is given the representation ae and represents the 


transformation from cartesian coordinate to normal coordinates 


Q=L <x (45) 


The concept of the dipole derivatives in cartesian coor- 
dinates, or the Atomic Polar Tensor (APT), was introduced in 
1961 by Biarge, Herranz, and Morcillo (7), and later reformu- 
lated by Person and Newton (8). The matrix of APTs for any 
particular molecule (termed P) is 3 by 3N and consists of 
the juxtaposition of N 3 by 3, sagend rank tensors, one for 
each atom in the molecule. Each element of the Nth APT is 
defined as (OP, /9%5), where - is the ath component (either 
x, y or z) of the total dipole moment, and 5 is the direc- 
tion of the displacement of the jth coordinate (also either 
X, Y, Or z). More detailed discussions on the formalism of 
the APT are available in the literature, (2,3,7,8,41,47,48). 
However, it should be remembered that the APT is simply a 
more easily interpretable representation of the dipole 
derivatives in normal coordinates. 

2.7 Characteristic Invariants of 

Atomic Polar Tensors 

The actual values of the elements of the APT depend, 
of course, on the orientation of the axis system used. For 
this reason it is imperative that the coordinate system being 
used always be reported along with the APT in order to properly 
interpret the APT. We begin our discussion of how to interpret 
the APT by starting with those properties which are invariant 


to the axis system utilized. 
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The mechanics of molecular vibrations impose two important 
conditions on the APTs. These conditions, and also the other 
invariant properties of the APT, will be briefly discussed here. 
A more thorough discussion is available elsewhere (55,63-65). 

The translational symmetry of the dipole moment requires 
that the APTs for a molecule sum to a diagonal tensor; each 
element in this tensor is equal to the total charge in the 


molecule. From (55) 


= qt _ 
: [¥,(P)]" = qi £ (46) 


where the summation runs over all the atoms in a molecule, 
Yn is the cartesian gradient operator (3/3x, 3/3y, 23/32), 
which when applied to a vector quantity (P, the dipole moment) 
gives a 3 by 3 tensor. E is the identity matrix whose only 
nonzero elements are ones along the diagonal. For a neutral 
molecule this reduces to the requirement that the elements 
of the APTs sum to zero, which is commonly expressed as the 
null condition 
z [v, (P)]* =0£ (47) 
gO ‘t 
A second important condition is an intensity sum rule 
first derived by Crawford (63). This sum rule is most 
easily derived following (55) through an analogy with the 


molecular kinetic energy T, 


() (48a) 
ax 
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ope etsy f2e-\ 4c) 282) «tee (48) 


where x, y and Zz represent differentiation with respect to 
time. The momenta in cartesian space (aT/aXx) can be trans- 
formed to moment in normal coordinate space using the chain 


rule of differentiation. 


2 aX, 


av |. (ar \(ad), / at |fad) ,/ az \/ Xo (49) 
ax aQ |/\ax aw ax dX aX 


where w represents rotational motion, and Xo translational 
motion. In these coordinate the kinetic energy takes on 
the form 


me \* 1 \foer \2 . 1 for J?] (50) 
: at 3 + : 
i | 2Qj ’ aa /\ ~ o \ 3X 


Oa 
The relationship between equation 49 and the corresponding 


equation for the dipole derivatives is obvious 


(33 c 2) (29} e [22 (2s ; 2 oe (51) 


aX 
Their equivalence is complete because in the harmonic oscil- 


lator approximation 30/aX = 90/a9X and 30/3X = 9w/dX. There- 
fore an expression analogous to equations 48b and 50 for the 
kinetic energy can be written for the dipole derivatives 


1// ap \? .[ ob ap \2 ap \? 4 Afse 3 \2.: 74 
2 Ml ax] *\ dy] ™iaecd 1 = 2looc) * 2 (cee) he 
n Mw N Yy N i\ "SE a 


ag a (e) 
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The above equation can be simplified by summation over the 
three components of the dipole moment Pos Pvt Po¢ giving 


the intensity sum rule 


= \2 = \2 = \2 ‘é = \2 
Ll ae ee ee a ap_ |" , / aP (53) 
9Q; 


ee + 
wn YN \\ O*y 9¥x | #43 


i, 
or, by using equation 42 and multiplying both sides of equa- 


tion 53 by the appropriate constants 


N_ot 
DA, +9 = ts p(t 2 


x (54). 
nM 3c N My N 


where Xn defines the "effective charge" according to King 


= (\2 a ND = 
aP . dP + dP 


dX, 8Yn 82 


2 N \+ N 
= =(Px || Px (55) 


The term 2 defines the translational and rotational correc- 


tions in equation 54 


2 2 2 2 2 2 2 
N__1\/P “+P P “+P P_ +P 3q 
— te x. y.. 4! ae Ree 2 ae e (56) 
ss fo ZZ XX yy fe) 


where the first three terms in the brackets of equation 56 
represent the contributions to the dipole moment derivatives 
due to rotation of the axis system, and the last term is the 
correction due to translation. This last term arises only 
when the permanent total charge q, on the molecule is not 


equal to zero. 
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For any molecule the LHS of equation 53 is invariant to 
a similarity transformation. Therefore the quantities on the 
RHS of this equation must also be invariant to similarity 
transform. This is equivalent to stating that the sum of 
the square of the elements of the APT for atom N is invar- 
lant to axis rotation. 

The importance of the effective charge as defined in 
equation 55 is twofold. First, it is invariant to a similarity 
transformation, and therefore is not dependent upon the axis 
system defined. Ideally, these quantities should then be 
transferable from an atom in one molecule to the same atom 
in another molecule (as long as this atom is in similar chem- 
ical environments in both molecules). Secondly, the effective 
charge is important because from equation 54 we can now see 
that the effective charge for an atom (normalized by the re- 
cripocal of the mass for that atom) measures the contribution 
of that atom to the total intensity sum for the molecule. 

A second quantity useful in characterizing APTs is the 
mean dipole moment Pye From (55) 
5s = 5 Trace (Py) = 5 = + ~ + Ta, (57) 
The trace of any tensor is invariant to asimilarity transform, 
so the mean dipole moment is also invariant. 

The last invariant property of the APT is related to the 


two previously defined invariant properties according to 


8B = F(Xy =p ) (58) 
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where By is defined as the anisotropy. It measures the devi- 


ation of the APT from the constant diagonal tensor, pYE (55). 


2.8 Quantum Mechanical Interpretation 
of Atomic Polar Tensors: The Charge, 
Charge Flux, Overlap Theory (CCFO) 


The charge-charge flux-overlap (CCFO) model was first 
introduced by King (55,66,67) and is related to the equilib- 
rium charge, charge flux concepts of Decius (68). The reader 
is referred to the original references, and also some more 
recent work illustrating the application of this model in 
interpreting infrared intensities (13,39,41,56,69) for a more 
detailed description of this model. Here only a brief outline 
of the more important concepts of this model will be discussed. 

The exact quantum mechanical definition of the dipole 


moment within the LCAO-MO method is as follows (41) 


i ea FLOR, ‘ : nie Davy tal tals? LS 
HaYB 
where Uy represents the orbitals on atom A, YB represents the 
orbitals on atom B, cs is the position vector of electron 1 


with respect to nucleus A, D is the density matrix element 


HAVB 


between orbital Ua on atom A and orbital v., on atom A, and 


B 


Q, is the gross atomic charge on atom A. The electronic pos- 


A 


ition vector may be defined as specifying the position of an 


electron, and may be expressed as 
Yr, =R, +2 (60) 


where ae is the position vector of ‘nucleus A and ry is the 


electron position vector. 
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Qn is the gross atomic charge on atom A and is given by 


Gr te 52 isd se (61) 
v 


2 
verry 
B Uy B AUB 7eAB 
The term D., s is the density matrix element which is defined as 
AB 
n 


D =2 EZ c¢ c (62) 
B iad) Px? Ye 


and § is the overlap integral 
“a5 


Ss =<u lv > (63) 
Uap A' B 


For simplicity, equation 58 can be rewritten as 


Q,R, + = F(A,B) (64) 


pes 
A A,B 


where F(A,B) is now given by 


F(A,B) = sie Pug Aa, GI, (65) 
The first term in equation 64 is obtained from the net charges 
located at the nuclear positions and can be assigned to a 
Single atom. The second term is essentially a hybridization 
term which measures the contribution due to the displacement 
of charge away from the center of the nuclear position... This 
effect is proportional to the off-diagonal density matrix 
element between orbitals centered on atoms A and B_ and is 


therefore intrinsically associated with a pair of atoms (41). 
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The APT for atom A is derived by applying the gradient 
operator (with respect to the cartesian coordinates of atom 


A) to the dipole moment vector of the molecule 


=, 4+ 
2 APTIC) = CV. (P) ] (66) 


Application of this equation to the quantum mechanical defini- 


tion of the dipole moment gives (Equation 64) 


Ct. —— ab = + 
li FLV, (Q,Rq) J 2 ELV Pee Bl (67a) 
Cc. _ -y qt - + 
Pp. =Q. 8+ ERALV(Q,) J 7 ZEA F (A,B) ) J (67b) 
= p°(c) + Po(cF) + P&(0) (67c) 


Thus the APT can be broken down into three components, which 
are represented by the three terms on the RHS of equation 67b. 
These three components are the charge, charge flux, and over- 
lap tensors respectively. The physical meaning for each of 
these tensors will be briefly outlined below. 

The charge tensor is a diagonal tensor, each element 
being composed of the net atomic charge on the atom in 


question 


Qn = Z, - N (68) 


where Qn is the net atomic charge, N, is the Mulliken popula- 


A 
tion (70), and Zn is the number of electrons in the neutral 
atom. The charge contribution to the APT is static, and in- 
volves only one atom A. The charge flux tensor is defined 


(41) as 
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+ 


Cc = + oe 
S-(CE) = -y RLV Ate-)) . = x 
ee x 2X0 


{ —o -— =O. .— 
N, (Xa-AX,) -N, (X,+AX,) 
Cc 
A Cc d eg Ot Cc (69) 


where N, (Xe + AX) is the Mulliken population for atom A after 


for atom C along the positive 


direction, N, (Xe - AX) is the Mulliken population of this 


displacement of coordinate x 


same atom after displacement along this same coordinate in 
the negative direction, Ra is the vector describing the equi- 


librium cartesian coordinates of atom A, Axe is the numerical 


=O 
Xar 


runs over all atoms in the molecule. The vector Ry is then 


value of the displacement of coordinate and the summation 
multiplied by transpose of the vector quantity given in paren- 
thesis in equation 69 to give the charge flux tensor for atom C. 

These quantities measure the contributions from each atom 
to the dipole moment change due to the redistribution of charge 
among all the atoms as one particular atom is displaced along 
one of its cartesian coordinates. In effect, these terms arise 
as a result of the changes in the atomic charges as one atom 
moves in either the x, y, or z direction. The contribution to 
the dipole moment change is obtained by "normalizing" the 
changes in the atomic charges by the respective atomic 
coordinates. This contribution can be considered as a "dia- 
tomic contribution" (41) because the quantities in equation 
69 involve a pair of atoms A and Ce 

The last component of the APT is the overlap tensor. 


This is a "triatomic contribution" (41) involving atoms A, 
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Band C. Physically, the overlap tensor corresponds to the 
changes in hybridization between atoms A and B as atom C is 
displaced along one of its cartesian coordinates. The effects 
on the dipole moment derivatives come about due to the changes 
in the displacement of charge from the nuclear centers. This 
last tensor is usually obtained from the difference of the 
total APT and the sum of the charge, and charge-flux 


contributions. 
C a AIC Cc 
P~(overlap) = = [Pp (charge) + P~(charge-flux)] (70a) 
P°(0) = P.© - ef(c) - P&(cr) (70b) 


It can be shown by analogy with the total APT that each 
of the three component tensors also has invariant properties. 
For the charge tensor this is already obvious. The charge ten- 
sor is a constant diagonal tensor, therefore the trace and the 


sum of the squares of the diagonal elements are also invariant. 


2 
N 7 eins N pe 
[x ic) a ate (aq) =, (71) 
prc) ai Trace (Charge Tensor) = ea zr PNc = Q 
3 363 ii N (72) 


where PNc) is a diagonal term of the charge tensor for atom 
if. 


N, given by Qn: The trace of the charge tensor is termed the 
mean atomic charge, while the sum of the squares of the ele- 
ments in the charge tensor is called the effective atomic 
charge. The anisotropy of the charge tensor is by defini- 


tion always equal to zero: 
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N 
B(C) = 0 (73) 


The charge flux tensor has been defined in equation 69. Using 
this definition and the chain rule of differentiation, the 


following equation can be written 


22a) _ [2a\ fag . (22a) faw) , (22a) [2% ee 
ae 3Q | |x Jo | (dx 3X, | 3x 


where all the symbols have their previously defined meanings 
and where the last two terms go to zero because atomic charges 
are not dependent upon either translation or rotation. By 
comparing this with equations 51 and 52, the relationship 


given below can be written 


safes) fa (2th leo 
n “y N Yn N i i 


where it is obvious that the RHS of equation 75 is invariant 
to similarity transform, so therefore the LHS of this equa- 


tion must also be invariant to similarity transform. If both 


sides of equation 75 are premultiplied by the quantity Ry = 
(x,*+y, 7427), then the following equation results 
N | a oN IX 9Yn 3Zy 
ae) / 
3Qn 
_ a 2; 3 
= ‘fs . (xX, +, +z") (76) 


Yq +247) is the square of the distance of 


atom A from the origin, and this distance is invariant to 


The quantity (x, 7+ 
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rotation of the axis system (similarity transform). Therefore, 
a product on the LHS of equation 76 is invariant to a similar- 
ity transform since both quantitites on the RHS of this equa- 
tion are invariant. If the LHS of equation 76 is expanded, 

the following result is obtained 

: 4: 2 3Q, 2 , 2 2Q, 2 ie 2 | 22, 2 io 2 | 32, \2 
N My A dX A 9Vn A dZy A dX 


2 { 20, \2 2 [ 22, \2 of \2 4} 2 eae 2/22, \2 
Ya Yn A | 32, A | ox, ay A 32, 
/ 
(77) 


If we now sum the RHS of this equation over all atoms A, and 

add the resulting terms, we obtain the squares of the nine ele- 
ments of the charge’flux tensor for the Nth atom. Since we are 
summing invariant quantities, then the total sum is also invar- 
iant to axis rotation. Therefore, we conclude that the sum of 
the squares of the elements of the charge flux tensor is invar- 
iant to axis rotation. This invariant quality is analogous to 
the effective charge for the APT, and is termed the effective 


charge of the charge flux tensor, (CF). 


Qf \2 O- \2 Q 
2 A Fe A * »¥| (78) 


a ae asain eae 


A Zz 

A *N Yy on] | 
Furthermore, as is true for any tensor, the trace of the 

charge flux tensor for the N atom is invariant to axis 


rotation. As a result, the trace of the charge flux tensor 


becomes the second invariant property of that tensor. 


ob 


_N ‘ , a 
P(CF) = 5 [Trace (Charge Flux)] = -_ P(CF 


Neg (79) 


uF 


aucel: is a diagonal element of the charge flux tensor 
for the Nth atom. The third invariant property of the charge 
flux tensor is derived from the previous two invariants and 
is termed the anisotropy etch). It is expressed as 


N 
(Mcef =3 ‘ol -[Pren (80) 


A similar argument can also be made for the overlap 
tensor. Each term of the overlap tensor is defined in equa- 
tion 67c, where F(A,B) is given in equation 65. An equation 

similar to equation 73 can be written for the overlap terms 
he ~ jem = - aX 
aF(A,B) _ {2F(A,B)|{29Q\ , (LF(A,B))| /aw\ , (2F(A,B) ° 

aX 3Q aX aw 3X dX, aX 


, 


(81) 
where as in the case of the charge flux tensor, the last two 
quantitites in equation 81 are equal to zero. This makes the 
correspondence between equation 74 and equation 81 exact, and 
the corresponding argument leads to the following result: 


2 
= 2 ee 2| - 2 - 
5 & Ee) . Rin 8) F de) ee oF (a8) (82) 
n ‘My N Yn N 


a a 
j f 


Here the LHS of this equation is equivalent to the sum of the 
squares of the overlap tensor. Since the RHS of this equation 
is dependent only on the form of the normal coordinate and is 


independent of the axis system used, then the LHS also becomes 


52 


invariant to axis rotation. The sum of the squares of the 
overlap tensor is then invariant to axis rotation and is 
termed the effective overlap, defined here in the LHS of 
equation 82. Again, as is true for any tensor, the trace 
over the overlap tensor becomes the second invariant property 
of that tensor. The trace of the overlap tensor, described 
in an analogous way to the traces of the previous tensors, 


is termed the mean overlap, and is defined as 


pY (0) = 


eS) 


N 
(Trace (overlap)] = u (0); (83) 


N 
Hhicce (oo); is a diagonal element of the overlap tensor. 


Finally, the anisotropy of the overlap tensor is also invar- 


lant to axis rotation and is defined as 


(so) : 4(kSoF - (5% 0}?) (84) 


The three component tensors of the APT, namely the charge 
charge flux,and overlap tensors, each has its own invariants 
analogous to the total APT. The sum of the squares of all of 
the elements for each individual tensor is invariant to axis 
rotation. For the total APT this quantity is the "effective 
charge," for the charge tensor the "effective atomic charge," 
the charge-flux tensor gives the "effective charge flux," and 
the overlap tensor yields the "effective overlap." The rela- 
tionship between the last three terms and the "effective 


charge" as defined by King and Mast (65) is given below 


\N = xX (c)? + (cr)? + x (0)? + Residuals (85) 
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where the residuals come about because the sum of the squares 
do not equal the square of the sum. These residuals are de- 


fined as follows: 


N N N N 
Residuals = : [2 P(C) ,,P(CF),. + 2 PAC) 557 (0) 55 


N N 
+ 2 Pie) PCO), (86) 


where the sum extends over all nine elements of the respective 
tensors. 

The trace of the total APT has been previously defined as 
the mean dipole moment in equation 57. The trace of the charge, 
the charge flux, and the overlap tensors are defined as the 
mean atomic charge, the mean charge flux, and the mean overlap 
and are defined in equations 72, 79, and 82 respectively. The 
anisotropy for each tensor derives from the two previous invar- 
iants for each tensor, _ and a, and is defined for the 
charge, charge flux, and overlap tensors in equation 73, 80, 
and 84 respectively. The significance of each of these 


terms will be discussed later. 


CHAPTER 3 
THE WATER MONOMER 

The frequencies and intensities were calculated for water 
interacting with a number of different chemical environments. 
The different environments were chosen to simulate different 
experimental conditions. Among those investigated were water 
interacting with another water molecule in several different 
geometries. These orientations include the linear, the cyclic, 
and the bifurcated pair of interacting water molecules, as 
shown in Figures 1.1, 1.2, and 1.3 respectively. In addition, 
the spectrum of water isolated in a beta-alumina complex was 
modelled by performing calculations on H50 interacting with 0 
atoms, out atoms, and 07? atoms. Finally, the spectrum of 
ionized water, BD" j was also calculated. The results calcu- 
lated for these species were then compared with those results 
calculated for isolated water. The changes in the calculated 
spectral parameters relative to the isolated molecule were 
dependent on both the geometrical orientation and the chemical 
environment of the water molecule. Therefore, calculations 
were done first on the isolated water molecule, and then these 
results were compared with the results obtained for a water 


molecule that was allowed to interact with its environment. 
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3.1 Calculation Procedure 

For water and all other molecules the force constants and 
dipole derivatives in § coordinates were calculated on the 
Amdahl 470 V/6-II computer of the Northeast Regional Data 
Center (NERDC) operating with OS MV/SP-JES2 Release 2 and 
using the Gaussian 76 system of ab initio Soiicilar orbital 
programs (54). This program was modified to allow the calcu- 
lation of the dipole moments of ions. This modification was 
first suggested by Rogers and Hillman (56), and permits the 
program to calculate the dipole moment for charged species in 
the same manner as done in the case of neutral molecules. The 
standard 4-31G basis set was used for all molecules (71). The 
integral threshold was set at 1.0 x tor: and the convergence 
limit at 5.0 x sO" % These values are standard for use with 
Gaussian 76. There has been some work on the dependence of 
the calculated results on the values used for these limits 
(72), but this problem was not considered here. The molecular 
geometry was input using principal cartesian coordinates, and 
the density matrix for the converged set of wavefunctions was 
stored on the IBM 3380 disk pack at the NERDC, to be used as the 
initial guess wavefunction for the calculations at the dis- 
placed geometries. This procedure reduces the amount of com- 
puter time needed for the ab initio calculations on the 
displaced geometries. Typically, a calculation on the 
equilibrium geometry would take 20 iterations before conver- 
gence was met, while a calculation on a displacement geometry 


would only take three or four iterations if the density matrix 
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for the equilibrium geometry was used as the initial guess. 
Two programs, ORTHO and PUNCH, were written for the Amdahl 
470 to set up the appropriate disk data sets at the NERDC 
for input to Gaussian 76. ORTHO writes the U transformation 
matrix to a 3380 disk pack. PUNCH uses this matrix as input 
and performs the matrix algebra described by equation 2b. 
This results in the calculation of the cartesian coordinates 
of the displaced geometries. These coordinates are then 
written to an IBM 3380 disk pack in the proper format for 
input to Gaussian 76. Finally, Gaussian 76 is used to cal- 
culate the energies and dipole moments of the equilibrium 
geometry and the displaced geometries. These results are 
then used in conjunction with equations 3, 4 and 5 to give 
the force constants and dipole derivatives in § coordinates. 

These two programswere designed to complement each 
other, each program creating an IBM 3380 disk data set for 
use in the next program. This eliminates some of the errors 
that can occur if one has to type all the data into the com- 
puter by hand (either using cards or a remote terminal). 

After the ab initio calculations were completed, the cal- 
culated force constants K, and the calculated dipole deriva- 
tives Por were transferred to an 8-inch floppy disk for use 
in the normal coordinate analysis. This analysis was performed 
on this laboratory's S-100 bus microcomputer. This system in- 
cludes the CompuPro/Godbout 8085/8088 CPU integrated to oper- 
ate with CP/M 2.2 (73) and CP/M 86. Two Shugart SA 801R disk 


drives for use with 8-inch floppy disks are used for data and 
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program storage. Other peripherals include a Televideo 910 
terminal, Centronics 739 dot matrix printer, Hayes Stack 
Smartmodem 1200 (for connecting to the NERDC) and a Micro- 
angelo Graphics Subsystem. The graphics subsystem includes 

a 15-inch high resolution graphics screen (480 x 512) detach- 
able keyboard and light pen. In addition, there is available 
128K of S-100 static memory, of which 64K is accessible by 

the 8 bit 8085 CPU. Software which is currently used operates 
with the 8 bit processor (8085) and includes Microsoft FORTRAN 
80, Microsoft BASIC 80, Microsoft BASIC COMPILER, and Micro- 
soft muMATH/muSIMP 80 (74). All of the normal coordinate pro- 
grams used here (see Appendix A for a brief description of 
each program) operates with the 8 bit 8085 CPU and 64K of 
S-100 static memory. A whole series of programs was written 
to perform the normal coordinate analysis. These programs 
were designed to be run in sequence; each program creates 
8-inch floppy disk data sets for use in the following programs. 
The proper sequence of programs necessary to perform a complete 
normal coordinate analysis is given in Appendix A, along with 
the necessary data input for each program (either from an 8-inch 
disk data set or by the operator via the Televideo terminal). 
All of the programs were written in either Microsoft FORTRAN 
(74) or Microsoft BASIC (74). The library subroutines neces- 


sary to run the main programs are also listed in Appendix A. 
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3.2 The Predicted Spectra for 

the Water Monomer 

A bond length of 0.951 A° and a bond angle of 111.2° was 
used for the calculation on the isolated water molecule. The 
numbering of the atoms and the orientation of the principal 
cartesian axis system used in these calculations is shown in 
Table 3.1. This geometry corresponds to the absolute minimum 
on the potential curve (theoretical geometry) for the ab initio 
SCF calculation of isolated water with the 4-31G basis set by 
Lathan, Hehre, Curtiss and Pople (75). There is some debate as 
to whether it is best to use the theoretically or experimentally 
determined geometry when calculating the spectroscopic proper- 
ties of molecules. However, for some of the systems investi- 
gated here there are little or no experimental data available. 
As a result, in order to be consistent, the theoretical geom- 
etries (calculated at the SCF level with the 4-31G basis set) 
were used whenever possible. When this geometry was not avail- 
able, other available information (either from other calcula- 
tions at different levels or from experimental work) was used 
to formulate a "best guess" at a reasonable geometry. This 
was done rather than searching the potential surface for the 
absolute minimum (i.e., optimizing the geometry) because the 
main interest is in spectral changes of the interacting species, 
and not absolute numbers. Therefore, calculations for the iso- 
lated species (in this case water) and then for the interacting 
species were performed at the same level (ab initio level with 
the 4-31G basis set), and then the predicted differences were 


examined. 
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The principal cartesian coordinates for the theoretical 
geometry of the isolated water molecule (R = 0.951 A°, 6 = 
111.2°) are given in Table 3.1. 

In order to do a full calculation of the frequencies and 
intensities, all the diagonal force constants, all the off- 
diagonal force constants, and all the dipole derivatives 
must be calculated according to equations 3, 4, and 5, 
respectively. Since the force constant matrix in $ coordinates 
is block factored by symmetry, some of the off-diagonal force 
constants are equal to zero and do not have to be calculated. 
However, all the off-diagonal force constants within the same 
symmetry block must be calculated. 


The water molecule has C symmetry. The 9 degrees of 


2v 


freedom reduce to 3A, + Az + 2B, 


tional coordinates reducing as 2A, + Bo: According to equa- 


+ 3By, with the three vibra- 


tion 3, two displacements are needed for each calculation of 
the three diagonal force constants, one in the positive dir- 
ection and one in the negative direction. This is a total of 
six displacements that are needed to calculate the three diag- 
onal force constants. The By symmetry block is 1 by l, and 
contains no off-diagonal force constants. The Ay block is a 
symmetric 2 by 2 block, and therefore there is one non-zero 
off-diagonal force constant. To calculate this off-diagonal 
force constant using equation 4, one additional displacement 
is needed, the simultaneous positive displacement of both 


coordinates. This makes a total of seven displacement geometry 


and one equilibrium geometry calculations that are necessary 
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Table 3.1 Principal cartesian coordinates for the isolated 
water monomer at the theoretical geometry (A°). Atom number-. 
ing is defined below. 


Z 


se 


eh 
Hy Hy 
xX ¥ Z 
Hy 0.0 -0.7847 -0.4772 
05 0.0 0.0 0.0611 
H 0.0 0.7847 -0.4772 
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for the full calculation of the frequencies and intensities 
for the isolated water molecule. 

The calculated frequencies and intensities for the water 
monomer at its theoretical geometry are summarized in Table 
3.2, where the experimental values (76,77) are also listed 
for comparison. As is usually the case with ab initio cal- 
culations, the predicted frequencies and intensities do not 
agree exactly with the experimental values (17,78). This is 
not too surprising considering the assumptions made and the 
errors inherent in these types of calculations. First, we 
have assumed that the molecule is a pure harmonic oscillator, 
and that the potentical curve can be well represented by a 
parabola. In reality, the potential curve is anharmonic, 
and the experimentally observed parameters reflect this 
anharmonicity. We have chosen to ignore this anharmonicity 
in our calculations. Strictly speaking this is not correct; 
therefore, the comparison between our calculated harmonic 
frequencies and the experimentally observed anharmonic fre- 
quencies cannot be expected to be in complete agreement. 

(The effect of anharmonicity can be quite significant, some- 
times up to 200 cm). ) Secondly, the 4-31G basis set employed 
in these calculations is not large enough to correctly describe 
the electronic wavefunction. The inadequacies of this basis 
set (such as the improper behavior of the gaussian function 

at small distances, the absence of polarization functions and 
the inflexibility in the orbital exponents) are significant 


enough that the electronic wavefunction does not describe 
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Table 3.2 Predicted wavenumbers and intensities for isolated 
water compared to the experimental wavenumbers and intensities. 


Wavenumbers v (cm -) Intensity A (km/mole) 


Mode Calc. Exp.* Cale. Exp.” 


Harmonic Observed 


ms 3983 3832 3657 3.90 2.24 
Vo 1676 1648 1595 125 54 
V3 4104 3942 3756 54 45 


@Reference 58 


PReference 41 
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correctly the electronic behavior of this molecule, leading 
to additional errors in these calculations. Finally, the 
effect of neglecting electron correlation in these calcula- 
tions is not entirely known. However, from some preliminary 
results by Bartlett (79), it appears that more sophisticated 
calculations may lead to much better agreement with experi- 
mental values. He has calculated the APTs for the isolated 
water molecule using "many body methods." The agreement be- 
tween his calculated results and the elements of the experi- 
mental APTs given by Zilles (41) is within 0.05 e. This re- 
sults suggests that the disagreement between the experimental 
frequencies and intensities and our calculated values isa 
result of the assumptions made and the methods used, and that 
this disagreement can be reduced by using more advanced theo- 
retical methods. en as a first approximation, the re- 
sults one can obtain using ab initio methods with the 4-31G 
basis set are generally accurate enough to be useful in in- 
terpreting infrared spectra (17). For this particular type 
of calculation (SCF with 4-31G basis set), the frequencies 
are expected to be within 10-20% of the experimental ones (71), 
and the intensities are predicted to be within a £actori*or 2 
of the experimental values (80). The wavenumbers and intensi- 
ties listed in Table 3.2 predicted for the water molecule with 
the theoretical geometrical parameters do indeed reflect these 
limits. 

Calculations of the frequencies and intensities were also 


done for the water monomer with varying bond angle and bond 
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lengths. Calculations were performed by varying the bond length 
while fixing the bond angle, by varying the bond angle while 
fixing the bond length, and by varying both the bond angle and 
bond length simultaneously. In addition, the frequencies and 
intensities were calculated for the experimentally determines 
geometry, where the bond length is 0.9572 A° and the bond angle 
is 104.5° (76,77). The calculated wavenumbers and intensities 
for these various geometries are summarized in Table 3.3. 

The data in Table 3.3 indicates that the values of the 
calculated frequencies are fairly sensitive to the assumed 
geometry. The calculation at the experimental equilibrium 
geometry predicts wavenumbers in much better agreement with 
the observed values than does the calculation at the theoreti- 
cal geometry. Changes as small as 0.006 A° in the assumed bond 
length for these calculations are responsible for changes of 
100-150 cmt in the predicted frequencies for the bond stretches. 
The bending frequency is less sensitive to changes in the 
assumed molecular parameters, but is nevertheless still 
affected. Surprisingly, a change of slightly less than 7° in 
the assumed bond angle results in drastic changes of nearly 
300 cmt for the two stretching modes. However, while the 
results on the frequencies indicate that close attention must 
be paid to the assumed geometry to predict reliable absolute 
numerical values, the predicted intensity values are much more - 
independent of the assumed molecular structure. The intensities 
are predicted to be nearly the same for each of the normal modes 


for all the different assumed structures. The largest 


65 


0S 8S 9p [x4 Sz es bs (®TOu/uy) 


680b Tete bbov S86E ZLep Lt0b bolp’ (,_™) f 
€2t 9zT 9Tt 86 86 szt Szt (®TOur/ury) % 
TST S99T b69T €6LT Z69T S69T 9L9T Go») ss 
08°” 00°F 08*z oe°e 06°b o9°€ 06°€ (9 Tow/ury) 7 
968E €66E£ 826E LO6E 8ZEb 968E €86E GC) : 
(a03;da00y 
uo1309eTg) (t0u0q 
e°Itt=0 4o2398T a) (Te3uewy 138dxq) (TeoTq020843) 
4 = 
i 856°0=-4 oc lI=0 of ‘OTT=0 ofS “bOT=0 ofS bOT=6 of TlT=0 of TTl=0 2POW 
ov os6‘o= "4 ov os6*="u="y ov 1S56°=Cu="y ov z156*="u="y ow 1S6°="u="y oW z156*="u="y ov 1s6"0=Cu= by 


*eTNosTOuU 19zeM 9Yy} FO 


"€ STqeL 
SOTIZOWOSH powmsse juUuSeTeFJIpP AOJF SOtT}TSusJUT pue sTSquNuSsAeM pdzZoOTperd EE 


66 


difference is the change by a factor of 2 in the predicted 
intensity of the asymmetric stretching mode as the asummed 
bond angle is change by 7°. However, it is obvious that in 
general the predicted frequencies are very sensitive to 
changes in the assumed geometries, while the calculated 
intensities are very much less sensitive to these same 
changes in geometry. 

At first glance, the strong dependence of the wavenumbers 
on the assumed geometry might seem very discouraging concern- 
ing the possibility of calculating absolute absorption fre- 
quencies because the final results are strongly dependent 
upon the initial assumed geometry. This problem has been 
examined in the past (12-16), but is not of interest here. 

In this thesis we shall be interested primarily in the 
changes that come about as a result of the interaction of 

the water molecule with its environment. Hence, the internal 
geometry of the water molecule can be chosen arbitrarily, and 
the intermolecular distance and orientation between the water 
molecule and the interacting species varied. In this way the 
resulting changes in the predicted spectrum can be attributed 
to the interaction and not a consequence of a geometry change. 
By fixing arbitrarily the internal geometries of the water 
molecules and varying only the intermolecular orientation or 
distances, we can then compare the predicted wavenumbers and 
intensities to observe their dependence only on the varying 
intermolecular parameters. Any dependence on the initial 


assumed geometry is eliminated because the internal geometry 
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Table 3.4 Predicted APTs and invariants for isolated water 
at the theoretical geometry, in e, where le= 1.602 x 10719 a. 
The coordinate system used is defined below. 


Z 


= 


oN, 


H 


482 0 0 y* 51360" e" 
H 

BL 1s 0 274 -.182 See 1072 ee 
0 ~.095 226 BY =. 50995-0° 
-.960 0 0 x? = .4754 a2 

0 2 
P. 2 = 0 ~.548 0 B* =) .4268 e* 
0 0 ~+452 gta 2187 at 
482 0 0 x? = .1292 e 
ee a 0 274 +.142 p* = .1071 27 
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is eliminated because the internal geometries of all the water 
molecules are assumed to be the same. 

The APTs (2) for the three atoms of the water molecule 
are presented in Table 3.4. The results presented here are 
for the principal cartesian coordinate system, but the coor- 
dinate system can be easily rotated to obtain APTs in any 
other coordinate system using equation 35. The invariants 
of the total APTs, Xn? pN and Bye are presented also. The 
breakdown of the APTs into their ‘three component CCFO tensors 
is given in Table 3.5. The overlap tensor is calculated by 
equation 70b as the difference between the total APT and the 
sum of the charge and charge flux tensors. All of the effec- 
tive invariants for the constituent CCFO tensors are given 
in Table 3.5 also. These invariant properties will be dis- 
siadad later. 

The internal force constants for the water monomer with 
the theoretical geometry were calculated from the force con- 
stants in S$ coordinates using equations 18 and 33. The inter- 
nal coordinates were defined as the two O-H bond stretching 
coordinates, and the H-0O-H angle bending coordinate. The B 
matrix elements were calculated using the little s vectors 
as described by Califano (11). The resulting B matrix was 
then used to calculate the A matrix using equation 32. The 
resultant force constants calculated with this A matrix are 
given in Table 3.6, along with experimental values for 
comparison. The B matrix elements used in these calculations 


are not weighted. This means that the units of the stretching 
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Table 3.5 Predicted CCFO tensors and their effective invariants 
for the isolated water molecule at the theoretical geometry, in 
e (le= 1.602 x 10719 c). The coordinate system is defined in 
Table 3.4. 


4017 0 0 ~.8034 0 0 
Hy O 
p(c) = 0 .4017. 0 P(é) = 0  =.8034 0 
0 0 .4017 0 0 ~.8034 
0 0 0 0 0 0 
P(CF) = | 0 -0550 .0370| P(CF) = 0 -.1090 0 
0 .1040 -.0450 0 0 0910 
0803 0 0 =, 1606 0 0 
wat Oo 
P(d) = 0 -.1827 -.1790 | P(d) = 0 3644 0 
0 -.1990 -.1307 0 0 2604 
2 
3 0 
x@l =}. 16168" XQ2 = .6455 e2 
2 2 
H 
: 2 O5 & Zz 
xt .0017 e X02 = .0067 e 
H.2 2 
Xl = .0190 e2 x02 = .0755 e@ 


The invariants for Hz are identical to those for H,. The com- 
ponent tensors for Hy are obtained simply by switching the 
signs on the off-diagonal terms for H}. 
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Table 3.6 Predicted internal force constants for the water 
monomer at the theoretical geometry compared to the experi- 
mental values. Units are md/A° where the weighting factor d, 
where d = 1 A°, is used to weight the force constants involvy- 
ing bending coordinates (see text and also reference 11). 


Internal Force Constants (md/A°) 


Calculated Experimental? 
2 R5 9. Lat 72750 
de3 0.714 0.668 
RR, -0.065 -0.068 
R) 83 0.300 0.270 


“References 58; fit to observed frequencies. 
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coordinates are in A°, and the units of the bending coordin- 
ates are rad (ll). As a result, the force constants calcu- 
lated using the A matrix derived from this B matrix are also 
not weighted. If the energy is expressed in md-A°, then the 
units of the stretching force constants are md-A°, the bending 
force constants are given in md-A°/rad’, and the stretch-bend 
interaction constants are given by md/rad. In order to compare 
these calculated force constants with other values in the lit- 
erature, care must be taken that the same weighting scheme is 
used in both instances (48). In the remainder of this thesis 
the weighting factor (11) d, where d = 1 A®° will be used to 
convert all force constants to units of md/A° without changing 
the numerical values of these constants. The normal coordinate 
transformation matrices obtained in these calculations, and all 


subsequent calculations, are summarized in Appendix B. 


CHAPTER 4 
THE LINEAR WATER DIMER 

The geometry for the linear dimer was taken from Curtiss 
and Pople (30), and is the minimum energy structure for 
(H50) 5 obtained by optimizing all geometrical parameters with 
respect to the total energy of the complex using the 4-31G 
basis set. This geometry will now be referred to as the 
theoretical equilibrium geometry for the linear dimer. The 
principal cartesian coordinates are given in Table 4.1. This 
molecule has Cy symmetry, having only a single symmetry element, 
the plane of symmetry. The calculated heavy atom distance 
(O--O) is 2.832 A°, slightly shorter than the distance deter- 
mined experimentally by Dyke, Mack and Muenter (81) of 2.98 A°. 
The discrepancies between the calculated value and the experi- 
mental value is due to the inadequacies of the basis set and 
the absence of electron correlation. Calculations performed 
with a larger basis set (including polarization functions) but 
without electron correlation predict a value for the O--0O dis- 
tance of 2.974 A°. This value is considerably closer to the 
experimental value than is the value calculated with the 4-316 
basis set. This suggests that the 4-31G basis set overestimates 
the strength of the hydrogen bond, but if the deficiencies of 
the basis set are eliminated then more accurate results are 


predicted. 
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Table 4.1 Principal cartesian coordinates for the linear water 
dimer at the theoretical equilibrium geometry (A°). Atom 
numbering as in Figure 1.l. 


Xx Y Z 
Hy =. 7876 -2207 1.8360 
Ho - 7876 -2207 1.8360 
03 0.0 -.0397 1«3729 
Hy 0.0 sO271 ~.4999 
Os 0.0 0613 -1.4573 
H 0.0 -.8109 -1...8337 
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The calculation performed by Curtiss and Pople (30) pre- 
dict the hydrogen bond to be linear within 0.05°, with an 
optimized energy of -151.83038 hartrees. In addition, the O-H 
bond length involved in the hydrogen bond is predicted to be 
0.958 A°. The remaining three O-H bonds are predicted to be 
all the same, and 0.950 A°. The H-O-H angles are calculated 
to be 111.3° in the electron acceptor, and 112.0° in the elec- 
tron donor. A calculation by Lathan, Hehre, Curtiss and Pople 
(75) with the same basis set gives a bond length of 0.951 A° 
and a bond angle of 111.2° for the isolated water molecule, 
with a total energy of -75.90864 hartrees. Therefore, the 
binding energy of the dimer (the energy of the dimer minus 
twice the energy of the monomer) is calculated to be 8.2 
kcal/mole. The internal geometries of the individual monomer 
units in the dimer are predicted to undergo slight changes 
upon complex formation. The bonded O-H bond length (H,-O;) is 
predicted to elongate by 0.007 A°, and the two internal bond 
angles increase slightly. While these changes appear to be 
very small, it is still surprising that the internal geometries 
of the individual monomer units are predicted to change. Most 
workers have assumed that the internal geometries of the two 
monomer units are unchanged from the values of the isolated 
molecule (81). Later in this chapter it will be shown that 
this geometry in the monomer units is very important towards 
understanding the spectral changes predicted for the dimer 
relative to the isolated monomer. 

The irreducible representations of all 3N degrees of free- 


dom for this molecule are 11A' + 7A". The 3N-6 or 12 vibrations 
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reduce as 8A' + 4A". For each of the 12 diagonal force con- 
stants two calculations at displaced geometries are necessary. 
This gives a total of 24 ab initio calculations that are needed 
to calculate all the diagonal force constants. For each non- 
zero off-diagonal force constant an additional ab initio cal- 
culation is necessary. In the A' symmetry block there are 
7+64+5+4+3+2+1 or 28 additional calculations necessary to de- 
rive all non-zero off-diagonal force constants. In the A" 
block there are an addition 3+2+1l or 6 off-diagonal force 
constants that must be calculated. This bring the number of 
ab initio calculations necessary to calculate the complete 
vibrational force field to 24+28+6 or 58 displaced geometries 
plus one equilibrium calculation, or a total of 59 calculations. 
4.1 The Predicted Spectrum for 

the Linear Dimer 

The predicted frequencies and intensities for the theoret- 
ical equilibrium geometry for the linear configuration of the 
water dimer are summarized in Table 4.2. The corresponding 
frequencies and intensities for the isolated water molecule 
are listed for comparison adjacent to those for the correspond- 
ing vibration of the linear dimer. In addition, the frequen- 
cies and relative intensities observed experimentally by Tursi 
and Nixon (19) are also listed here for comparison. The pattern 
predicted by our ab initio calculation is in good agreement with 
the experimental spectrum assigned to the linear structure by 
Tursi and Nixon (19). Behrens and Luck (27) argue that the in- 


terpretation of the spectrum of (H0) 5 by Tursi and Nixon (19) 
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is not conclusive evidence for the linear structure. Their 
arguments are based on the fact that all six of the bands 
assigned to the linear structure do not follow the same con- 
centration dependence, as they should if they all belonged to 
the same species. However, reliable assignments of these bands 
cannot be made solely on the basis of their concentration 
dependence because there are often overlapping and broad bands 
corresponding to higher aggregates which serve to hinder reli- 
able measurements of relative intensities. The calculation of 
Curtiss and Pople (30) along with other work (19,26,81) points 
to the linear structure as the most energetically favorable one, 
and using this structure in our calculations we have been able 
to reproduce reasonably well the experimental spectrum of Tursi 
and Nixon (19) (see Table 4.2). 

While it is probably true that these ab initio 4-316 cal- 
culations are still not accurate enough to give reliable abso- 
lute values, it is encouraging that the correct patterns are 
reproduced, including the spectral changes that occur upon 
complex formation. The nature of these spectral changes is 
of primary interest, and not the absolute values predicted 
for the frequencies and intensities. Since these ab initio 
calculations with the 4-31G basis set correctly reproduce the 
observed spectral changes, we believe that these calculations 
are valuable in understanding these changes. These changes 
in turn are important to understand because their nature might 
offer some valuable insights towards understanding intermolec- 


ular interactions. 
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The assignment of the vibrations as predominantly the 
electron donor symmetric stretch, or electron acceptor bend- 
ing, for example, were made with the aid of a program which 
was written to picture a vibrating molecule. The program, 
VIBRATE, accepts the equilibrium cartesian coordinates for 
each atom in the molecule and then draws a ball and stick 
model of that molecule. This portion of the program has been 
adapted for Our microcomputer from the PLUTO program of Mother- 
well (82). The cartesian coordinates corresponding to the 
displacements of an appropriate increment in both the positive 
and negative direction along the normal coordinates are also 
input from a disk data set into the program. Atoms (or circles) 
are then turned on and off at these displaced cartesian coor- 
dinates for any desired normal mode. These circles are turned 
on and off at both the positive and negative displacements, 
giving the illusion that the atoms in the molecule are actu- 
ally vibrating along that particular normal coordinate. By 
observing the motions of all the atoms in the molecule for any 
particular normal vibration, the characteristic motions of 
that normal mode can be accurately described. These vibrations 
can be frozen on the graphics screen, then sentto the printer 
for hard copy. This program was written in BASIC and uses the 
Microangelo high resolution graphics screen for the display. 
The low frequency modes were also assigned using this program. 
This type of assignment has certain advantages over the Poten- 
tial Energy Distribution (PED) (58) descriptions of normal 


vibrations because they allow one to actually visualize the 
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vibration. The PEDs can be misleading because they do not 
describe the actual normal vibration, but only the contribu- 
tions from a given set of internal coordinates that make up 
that particular normal vibration. 

At this point let us just briefly comment on the charac- 
teristics of the 12 vibrational modes for the theoretical 
linear configuration of two interacting water molecules. The 


a iL 


highest frequency modes predicted at 4126 cm ~, 4076 cm, 


1 and 3872 cm 


3992 cm (Table 4.2) are all stretching modes 
and are isolated almost entirely in only one molecule. These 
vibrations are the asymmetric stretch of the electron donor, 
the asymmetric stretch of the electron acceptor, the symmetric 
stretch of the electron donor, and the symmetric stretch of the 
electron acceptor, respectively. The intermolecular interaction 
does not couple or perturb the form of these high frequency 
stretching modes at all, and in fact they remain localized in 
one molecule. These normal modes behave as if two separate 
molecules existed without any interaction. The two modes 
occuring at 1800 om72 and 1714 om? correspond to bending 
vibration of the water units. However, these vibrations are 

no longer localized, but are coupled between both molecules. 
For instance, the vibration at 1800 om7+ is predominantly 

the electron acceptor bending motion, coupled in phase (i.e., 
both molecules open their angles simultaneously) with some 

of the bending mode of the electron donating molecule. The 


vibration at 1714 cm=+ is mostly electron donor bending 


coupled out-of-phase with some electron acceptor bending. 
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In this lower wavenumber region the vibrations are no longer 
localized on the individualized water molecules, but rather 
become vibrations of the entire complex. While it is possible 
to classify the vibration at 1800 cm and the one at 1714 cm! 
as electron acceptor and electron donor bending, respectively, 
strictly speaking this is not correct due to the presence of 
intermolecular coupling which perturbs and mixes these 
vibrations. These bending vibrations are not only coupled 

with each other, but they are also coupled to some of the low 
frequency intermolecular bending and twisting motions. 

The low frequency (or intermolecular) modes are generally 
characterized by twisting and out-of-plane (where reference is 
made to the plane of the molecule) bending motions. The vibra- 
tion predicted at 835 cm? is mainly the 0,-H,-0, out-of-plane 
(of the electron acceptor) bend (hydrogen bond bend) along with 
some electron donor twisting motion. The mode at 749 em? cor- 
responds to in-plane electron acceptor twisting in phase with a 
little out-of-plane (of electron donor) electron donor wagging 
(electron acceptor twists upwards at the same time the hydrogen 
atoms of the electron donor wag upwards). The vibration at 208 
om=+ is the out-of-phase combination of these same electron 
acceptor twisting and electron donor wagging motions just 
described. The mode at 313 cm + can be characterized as pre- 
dominately electron donor twisting motion mixed with electron 
acceptor hydrogen (free hydrogen) wagging motion. At 210 cm + 
the vibration corresponds to the out-of-plane (of electron 


acceptor) twisting of the electron acceptor mixed 
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with some in plane (of electron donor) electron donor bending. 
Finally, the vibration at 103 cm=+ is the hydrogen bond strech- 
ing mode. 

The low frequency modes are presented here for complete- 
ness, but will not be discussed at any length. There is very 
little experimental evidence about them available from studies 
in this region of the spectrum. The uncertainties in the ab 
initio calculations are magnified in this region because the 
force constants are extremely small. These frequencies are 
therefore expected to be much more sensitive to the assumed 
geometry, and the confidence that one can place on the pre- 
dicted values of these absolute frequencies is certainly 
limited. For some configurations some of the intermolecular 
modes were predicted to have imaginary frequencies. This re- 
sult is a consequence of performing the calculations at dis- 
placements that are not about the absolute minimum of the 
potential surface. It is not expected that these imaginary 
frequencies affect greatly the frequencies of interest, namely 
the high frequencies intramolecular modes, due to the large 
energy separation. However, it is certainly interesting to 
note here that the predicted intensities for the low fre- 
quency modes are in general equal to or greater than those 
for the intermolecular vibrations. This might offer some 
hope to experimentalists interested in investigating the 
bands in this region. However, until now, the difficulties 
(such as the presence of internal rotation, and a low signal 


to noise ratio) in this region of the spectrum have made 
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the study of these bands extremely difficult. Even so, the 
intensities of these modes merit a few comments. 

Four of the six low frequency modes are predicted to have 
a large intensity; greater in fact, than for the high frequency 
modes. As a result of comparing the motions involved in the 
high intensity low frequency modes, it was observed that these 
four vibrations all involved a considerable contribution from 
the atoms in the electron acceptor, while the remaining two 
low intensity modes were mainly motion of the electron donor. 
This suggests that perhaps most of the intensity for the com- 
plex originates in the electron acceptor. 

In general, the data in Table 4.2 indicates that the 
intramolecular vibrations of the electron donor are essen- 
tially not affected by the interaction with another water 
molecule. However, the electron acceptor experiences con- 
siderable change, especially in the symmetric stretch 
(v, (EA) ]. This band is predicted to shift 110 em> to lower 
frequency and to increase in intensity by a factor of almost 
72, from 3.90 km/mole to 279 km/mole. The remaining modes 
for this complex are not predicted to change very much from 
the corresponding vibrations in the isolated water molecule. 
Both asymmetric stretching vibrations are predicted to in- 
crease by slightly more than a factor of 2 in intensity, but 
this is certainly less than the factor of 72 predicted for 


the symmetric stretch of the electron acceptor. 
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4.2 Intensity Analysis for Linear 

Pair of Water Molecules 

The change in intensity predicted for the symmetric stretch 
of the electron acceptor is extreme, from a relatively weak 
band in the monomer to a strong band in the dimer. A detailed 
ab initio calculation and analysis of the linear dimer has 
recently been published by Zilles and Person (39). However, 
our equilibrium calculation here utilizes a slightly different 
approach so that our results are slightly different from those 
reported previously (39). For instance, we have used an un- 
scaled ab initio force field and unscaled ab initio polar 
tensors, while Zilles and Person (39) used the experimental 
force field of Tursi and Nixon (19) with scaled APTs. Sec- 
ondly, while our equilibrium calculation was performed at 
the theoretical geometry, Zilles and Person (39) chose to 
use the experimental geometry. While the use of an experi- 
mental force field, scaled APTs and the experimental Gecmmlees 
usually gives better agreement with experimental results, our 
calculation is based on purely theoretical methods and there- 
fore gives an indication of the quality of predictions that 
can be expected using these methods. Finally, we have ex- 
tended the work of Zilles and Person (39) by actually quan- 
tifying the predicted changes in the polar tensors to 
changes in intensities. In addition, we have investigated 
the dependence of the predicted spectrum on both intermol- 
ecular orientation and distance by performing similar calcu- 


‘lations at different geometries and intermolecular distances. 
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The results of our equilibrium calculation are very sim- 
ilar to the results reported by Zilles and Person (39). How- 
ever, there are some slight differences in the two approaches 
so that caution must be used in making a different numerical 
comparison between the two results. Zilles and Person (39) 
found that the only significant change in the APTs of the 
dimer relative to the isolated monomer was in the Py term 
of the bonded hydrogen. Therefore they attributed the large 
predicted intensity increase (in the dimer relative to the 
monomer) in the symmetric stretching mode of the electron 
acceptor as due to this change in this term of this polar 
tensor. When the CCFO analysis was performed on the APTs 
by Zilles and Person (39), they found most of the increase 
in the APTs for the dimer was due to an increase in charge 
flux as the dimer was formed. Our results are similar, and 
are briefly summarized below. 

In order to attempt to understand the large predicted 
intensity increase for the symmetric stretching mode of the 
electron acceptor, the intensity parameters were transformed 
to APTs. The APTs for the theoretical linear configuration 
in the coordinate system given in Figure 1.1, are given in 
Table 4.3, along with the properly rotated APTs for the theo- 
retical water monomer. First, it must be remembered that the 
intensities are related to the normal coordinates via equation 
24. Since the normal coordinates are mass weighted motions of 
all the atoms in a molecule, the heavier atoms will have smaller 


amplitudes and therefore contribute less to the normal coordinates. 
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Table 4.3 Calculated APTs for theoretical linear dimer alogg 
with corresponding APTs for monomer (e) (le = 1.602 x 107 Ch} 
Atom numbering and coordinate system as given in Figure l.l. 


Linear Dimer Monomer # 


3284 -.0777 -0596 -274 aay YS 074 

H 
Nes = -.0406 eaeok +1079 -.081 «295 pe bt Us 
0331 -0968 - 4388 - 040 ii3 412 
3284 +.0777 -.0596 -274 ~AZd 0.074 

H 
p 2 = - 0406 3231 «L079 O81 2295 me oe! 


~.0331 -0968 - 4388 


=.6055 0 0 -.548 0 0 

ss 7 
pS 0 1, ieee, <3 3156.1" 0 in, $9.1 #228 
- 9 ~yI633: 250 0 ~228 -.825 


u 


.4827 0 0 482 0 0 
Hy 
Ps 0 .5892 -.0688 0 148 -.089 
0 -.0803 .3700 0 -.042 See 
-1.0010 0 0 -.964 0 0 
0 
ag Se 0 -.8115 -.0078 0 oe 047 
0 .0433  -.5055 0 047 ~.483 
4670 0 0 482 0 0 
He 
= 0 36é5 .0763 0 369 044 
0 .0066  .1000 0 ~.003 2337 


@the APTs for the atoms in the monomer are for coordinate sys= 
tems rotated from Table 1 separately for each atom to corres- 
pond to the coordinate system of the dimer in Figure 1.1. 
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Therefore, when analyzing infrared intensities the lighter 
atoms are of greater importance than the heavier atoms. 
While bearing this fact in mind, the results in Table 4.3 
show reasonably large changes in the APTs of both the oxygen 
and bonded hydrogen atoms of the electron acceptor. The Poy 
terms for both of these atoms change by a substantial amount. 
The Poy term for the oxygen atom decreases by 0.308 e, while 
the Poy term for the bonded hydrogen increases by 0.441 e. 
Since the lighter atoms play a more important role in the cal- 
culation of infrared intensities, the following discussion will 
emphasize the bonded hydrogen atom when trying to explain the 
large predicted intensity increase in the symmetric stretching 
mode of the electron acceptor for the dimer. 

The polar tensors were further analyzed to determine 
which, if any, of the three component tensors was mainly re- 
sponsible for the change in the total APT. The charge, 
charge flux, overlap (CCFO) analysis was performed for all 
of the atoms in the linear pair of water molecules. The 
resultant charge, charge flux and overlap tensors are listed 
for the linear pair and the corresponding water monomer in 
Tables 4.4, 4.5, and 4.6, respectively. 

The three tables, Tables 4.4, 4.5, and 4.6 indicate that 
the charge and overlap tensors for all of the atoms are rela- 
tively constant in going from the monomer to the linear pair. 
There are slight changes, but nothing significant. However, 
_the charge flux tensors show most of the changes reflected in 


the total APTs. For the Poy term of the bonded hydrogen, the 
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Table 4.4 


10-19 gq), 
Figure 1.1. 


Linear Dimer 


.4237 0 0 .4017 
H 
1 
Gaye 0 4237 0 0 
0 0 -4237 0 
.4237 0 0 .4017 
H 
2 
i a 0 4237 0 0 
0 0 4237 0 
-.8184 0 0 -.8034 
03 
Pam 0 -.8184 0 0 
0 0 -.8184 0 


P(C) = 
+3812 


ue) 
=~ 
Qs 
i] 
a 
Oo oo 
w 
Le) 
v4 
a oOo 
Ww 
NO 
oO oO 
> 


Calculated charge tensors for the theoretical linear 
dimer compared to those of isolated water (e) 


(le = 1.602 x 


Monomer 


Atom numbering and coordinate system as given in 
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Table 4.5 Calculated charge flux tensors for the theoretical 
linear dimer compared to the corresponding §ensors for iso- 
lated water monomers (e) (le = 1.602 x 1071 C). Atom number- 
ing and coordinate system as given in Figure 1.1. 


Linear Dimer Monomer 
-0957 -0661 -.0370 ~053 +.034 -.012 
Hy 
P(CF) = «L293 -.0470 0152 _-089 -.033 -020 
-.0663 -0087 -.0132 -.054 -020 -.012 
-0957 -.0661 -0370 ' 053 -.034 ~O2ZF 
H 
2 
P(CF) = -.1293 -.0470 «OL5E -.089 -.033 -020 
- 0663 -0087 -.0132 -054 -020 =, 012 
-.1769 0 0 -.105 0 0 
RS 
P(CF) = 0 -.0760 -.0320 0 - 066 -.040 
0 50228 337 0 -.040 ~025 
-.0268 0 0 0 0 0 
Hy 
P(CF) = 0 +.4495 -.0276 0 089 -.051 
0 =.0352 -.0905 0 OL2 -.081 
0184 0 0 0 0 0 
Os 
P(CF) = 0 -.2354 -0436 0 -.043 O91 
0 -0919 -0499 0 091 029 
-.0061 0 0 0 0 0 
me 
P(CF) = 0 -.0442 -.,0141 0 -.046 -.040 
0 -.0960 -0334 0 -.104 054 
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Table 4.6 Calculated overlap tensors for the theoretical linear 
dimer compared to the corresponding tensors for isolated water 
monomers (e) (le = 1.602 x 10719 c). Atom numbering and coor- 
dinate system are given in Figure 1.1. 


-. 1910 -.1438 - 0966 -. 181 =,155 «095 
H 
P (5) aj =. 1699 -.0536 -0928 ~. 170 ='..074 093 
-0994 -0881 0283 +.094 093 022 
~. 1910 -1438 -.0966 =. 18. ak 5S =,.095 
H2 
P(O) = - 1699 -.0536 .0928 170 -.074 < 093 
-0994 -0881 -0283 -.094 -093 -022 
- 3898 0 0 
03 
P(O) = 0 -1060 -.1834 


- - 0 «1343 -.038 
0 -.0451 -.0031 0 -.054 .030 
-.1461 0 0 
Os 
P(O) = 0 .2972 -.0486 
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increase in the total APT is predominately due to the increase 
in the charge flux tensor, and not a result of a change in the 
charge or overlap tensor. Quantitatively, the increase in the 
charge flux tensor is responsible for 0.360 e, or 82% of the 
change in the Py term of the total APT for this atom. 

A program was written, PXPQ (see Appendix A), to determine 
exactly what percentage of the total intensity increase from 
monomer to dimer was due to the difference in the normal coor- 
dinates (or due to the mechanical effect), and what amount 
was a result of the change in charge, charge flux and overlap 
in the bonded hydrogen. This program uses a normal coordinate 
transformation matrix and a set of polar tensors as input. It 
allows the operator to adjust any of the elements in any of the 
polar tensors for any of the atoms as he wishes. The intensi- 
ties using the input normal coordinate transformation matrix 
and these adjusted polar tensors are then recalculated. If 
we use the rotated monomer polar tensors as the initial input 
APTs, and the calculated normal coordinate transformation for 
the linear dimer, the calculated intensities will show the 
effects due to the changing normal coordinates (mechanical 
effect), but not the electronic effects from the interaction. 
If we then add only the increase in charge flux for the Poy 
term of bonded hydrogen and recalculate the intensities, 
these intensities will reflect only the normal coordinate 
effect and the effect of the increase in charge flux for the 
bonded hydrogen atom. The results of both of these calcula- 


tions are given in columns four and fiveof Table 4.7. The 
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normal coordinate effect is minimal, as evidenced by only a 
slight increase in the intensities of the symmetric stretching 
modes, from 3.90 km/mole to 4.57 km/mole for the electron 
donor and to 11.6 km/mole for the electron acceptor. This is 
to be expected based on the fact that the stretching vibrations 
of the dimer were already seen not to couple or become per- 
turbed much by the interaction. When the charge flux increase 
for the Poy term of the bonded hydrogen is added to the monomer 
APTs, the predicted intensities for the electron acceptor ex- 
hibit considerable changes. The asymmetric stretch increases 
from 47 km/mole to 123 km/mole, and the symmetric stretch in- 
creases from 11.6 km/mole to 235 km/mole. The predicted in- 
tensities for this pair of vibrations in the complex are 110 
km/mole and 279 km/mole respectively. By adding only the 
increase in charge flux for the Poy term of the bonded hydro- 
gen of the electron acceptor (Hy) all of the predicted inten- 
sity increase for the asymmetric stretch of the electron 
acceptor for the linear pair can be reproduced. By comparing 
columns five and nine of Table 4.7, it is obvious that prac- 
tically all of the predicted intensity increase for the linear 
dimer is a direct consequence of the increase in charge 

flux in the Poy term of the bonded hydrogen. The effect 

of both the change in charge and overlap tensors on this 
intensity increase can be seen by successively adding the 


differences from dimer to monomer in the Poy term of the charge 
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and overlap tensors to those APTs used for column five. It is 
seen in columns six and seven that these contributions are 
indeed minimal; containing only 6% and 7% of the total pre- 
dicted intensity increase respectively. 

The only band of the linear dimer which is not adequately 
reproduced in column seven is the asymmetric stretch of the 
electron donor. It is still a factor of 2 too low compared to 
the predicted intensities in column nine. If we refer back to 
Table 4.3 and look again at the hydrogen atoms for the electron 
donor, there is an increase of 0.054 e in the Ry term of both 
hydrogen atoms in the dimer relative to the monomer. This 
change is relatively small, but is still large enough to be 
responsible for the increase in intensity for the asymmetric 
stretch of the electron donor. From Table 4.5, 0.0427 e or 
78% of this change is due to an increase in charge flux:.> The 
terms for both hydrogens were increased by this amount, and 
the resulting predicted intensities are given in column eight 
of Table 4.7. By adding the difference in charge flux from 
dimer to monomer to the eer term, the predicted intensity for 
the asymmetric stretching mode of the electron donor rises 
from 55 km/mole to 83 km/mole. The value of 83 km/mole is 
in closer agreement with the predicted value for the dimer 
of 110 km/mole. Therefore, by considering only the increase 
in charge flux in the dimer relative to the monomer, almost 
76% of the intensity of the symmetric stretch for the electron 


donor can be accounted for. Therefore, we conclude that the 
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intensities for the six intramolecular vibrations of the dimer 
are well reproduced by adjusting only the three terms in the 
polar tensors--the Bf term, the pul and the pu2 term--and 
that most of the change in these terms as the dimer is formed 
is due to a change in charge flux and not to changes in 
charge or overlap. 
4.3 Frequency Analysis for the Linear 

Pair of Water Molecules 

The predicted frequencies for the six intramolecular 
vibrations of the linear dimer (see Table 4.2) are in good 
agreement with the predicted values for the corresponding 
vibration in the isolated monomer, differing by only 20-30 
om: However, three of these vibrations are predicted to 
occur at wavenumbers considerably different from the corres- 
ponding vibration in the monomer. The symmetric stretching 
vibration of the electron acceptor is predicted to decrease 
110 cmt in the dimer relative to the monomer. The two bend- 
ing vibrations in the dimer are predicted to increase by 124 
cm-+ and 39 cm. In order to understand these changes the 
calculated frequency parameters or force constants were analyzed. 

The frequency parameters in S$ coordinates were transformed 
to internal force constants using equations 18 and 33. The 
internal coordinate definitions are given in Table 4.8, and 
the resultant B and A matrices were calculated with the pro- 
grams BMAT and AMATRIX respectively (see Appendix A). The 
weighting factor d, where d = 1 A® has been applied to all 


the calculated force constants and therefore the units of 
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Table 4.8 Internal coordinate definitions used for the linear 
dimer. Atom numbering as given in Figure 1.1. (d is a weight- 
ing factor where d = 1 A°.) 


Coordinate No. Definition 
Rj R (H,-0,) 
Ry R (H,-03) 
R3 R (H,-O.) 
Ry R (He-O;) 
Re d@ (H,-0,-H,) 
Re d@ (H,-O,-H,) 
Ro R (H,-03) 
Rg dé (H, -03-H,) 
Rg dé (H,-0,-H,) 
Rig Be KH 903 ~Og= By) Pee ee 
Ri dé (0,-H,-O.) (in plane) 
Ri2 dé (0,-H,-O,) (out of plane) 
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all these force constants are md/A° (48). The weighted in- 
ternal force constants for the theoretical linear dimer are 
given in Table 4.9. 

The calculated force constants for the intermolecular 
vibrations are all very small, generally having a magnitude 
about 0.05 md/A°. Therefore, the effect of these force con- 
stants on the high frequency vibrations is expected to be 
minimal, so only the force constants corresponding to the 
intramolecular vibrations will be discussed here. These 
force constants can be found in the upper left 6 by 6 corner 
of the matrix given in Table 4.9. The values for the first 
six diagonal force constants in the table represent either 
O-H stretching force constants or H-O-H bending force 
constants. These values can be directly compared to force 
constants calculated for isolated water because there exists 
a corresponding stretching or bending coordinate in the 
monomer. The calculated values for the theoretical geometry 
of isolated water can be found in Table 3.6. The most obvious 
and most significant difference is in the bonded O-H stretch- 
ing constant, which changes from 9.14 md/A° in the monomer to 
8.54 md/A° in the linear dimer. The exact nature of the effect 
of this change on the frequencies was investigated using the 
program PERTURB (see Appendix A). 

PERTURB accepts a set of internal force constants from 
an 8-inch floppy disk and calculates the corresponding wave- 
numbers using this set of force constants. It then gives 


the operator the choice of adjusting any of the force 
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constants, and recalculates the wavenumbers using these 
adjusted force constants. In this way we are able to in- 
vestigate how changes in individual force constants affect 
the predicted wavenumbers. 

The calculated internal force constants for the dimer 
were used as the initial input to the program PERTURB. The 
bonded O-H stretching force constant of the electron acceptor 
was then increased from 8.54 md/A° to 9.14 md/A°, or to the 
same value for the O-H stretch in the isolated water molecule. 
The frequencies were then recalculated with this force con- 
stant changes and the remaining force constants at their 
original calculated values. The results are listed in column 
three of Table 4.10. The frequency of the symmetric stretch 
for the electron acceptor increases from 3873 cm7+ to 3968 
cmt, and this can be compared to a frequency of 3983 cm + 
for the isolated water molecule. Almost 86% of the total 
predicted frequency shift in the dimer (relative to the 
monomer) of 110 cmt can be accounted for with the bonded 
O-H stretching force constant. As the electron donor perturbs 
the electron acceptor, the strength of the bond between the 
bonded hydrogen and the oxygen atom in the electron acceptor 
decreases, and therefore the force constant for this O-H 
stretch decreases. This decrease in turn is almost entirely 
responsible for the predicted decrease in the symmetric stretch- 
ing mode of the electron acceptor in the dimer relative to the 


same stretching mode in the monomer. 
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However, it must be remembered that the individual monomer 
units in the theoretical geometry for the linear dimer do not 
retain the theoretical monomer geometry. The calculation of the 
optimized geometry for the linear dimer (30) predicts that the 
geometry of both the electron donor and acceptor are perturbed 
(relative to the optimized geometry for an isolated water mole- 
cule) by the presence of the second molecule. The calculated 
frequencies and intensities for isolated water molecules at the 
geometry of the electron acceptor and electron donor are given 
in Table 3.3. For the electron acceptor the calculated wavenum- 
bers in the dimer are in good agreement with the same values 
calculated for an isolated water molecule with the same geometry 
as the electron acceptor. The largest difference is in the bend- 
ind mode, where the frequency is predicted to be 1800 cm + in 
the dimer, and is calculated to be 1751 cm™+ in the isolated 
electron acceptor. The asymmetric and symmetric stretching modes 


i and 23 cm respectively. This suggests 


agree within 13 cm 
that the predicted change in the vibrational frequency of the 
symmetric stretch of the electron acceptor of the dimer, relative 
to the isolated monomer frequency, is due to the change in geom- 
etry of the acceptor water molecule in the complex relative to 
the geometry in the isolated water molecule. Since we have al- 
ready shown that the predicted frequency shift for the symmetric 
stretching mode of the electron acceptor was due to the decrease 
in the stretching force constant of the bonded O-H, we conclude 


that this decrease in the force constant is due to 


the change in the geometry of the 
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acceptor molecule in the dimer relative to the isolated monomer. 
For the electron donor, the calculated wavenumbers in the com- 
plex are in fair agreement with the calculated values for the 
isolated water molecule at the theoretical geometry. However, 
when the frequencies for an isolated water molecule with the 
geometry of the electron donor are calculated, the frequencies 
of the stretching modes are in almost exact agreement with 

those predicted for the donor molecule in the dimer. The 
bending frequency for the electron donor is in somewhat worse 
agreement, still being 50 cmt too high in the complex. 

While we have been able to account for the predicted shifts 
in the stretching frequencies of the dimer (relative to the 
monomer) by the changes in geometry of’the individual monomer 
units, the same cannot be said for the predicted shifts in the 
bending frequencies. Even after correcting for the geometry 
changes, the bending frequencies predicted in the dimer are 
still about 50 cm + higher than expected. This is due to the 
coupling of these vibrations between the two molecules. If 
we recall the descriptions of the bending vibrations of the 
dimer given earlier, they were no longer isolated in one mole- 
cule but were instead mixed between both molecules. However, 
these bending vibrations are not only coupled to the bending 
vibration of the second molecule, but they are also coupled 
to the low frequency intermolecular vibrations. Column five 
of Table 4.10 gives the calculated frequencies using only the 
complete vibrational force field corresponding to the six 


intramolecular vibrations (i.e., the upper left 6 by 6 block 
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of Table 4.9). In this way, the interaction constants between 
the high frequency modes and the low frequency modes are 
assumed to be zero, and therefore no coupling between these 
modes can occur. The bending frequencies calculated in this 
way agree almost exactly with the values that one would expect 
based only on the geometry change. Therefore, we conclude that 
the difference in the bending frequencies in the dimer relative 
to the isolated theoretical monomer is a result of two effects. 
The first effect is the predicted geometry change of the two 
interacting molecules in the dimer relative to the monomer, 

and the second effect is the coupling between the bending 


vibrations and the low frequency intermolecular vibrations. 


4,4 Further Calculations on Weaker and 

Stronger Interacting Linear Pairs 

of Water Molecules 

There exists a great deal of evidence for the relationship 
between the energies of OH---O hydrogen bonds, as manifested by 
the bonded O-H stretching wavenumber, and the OH--O hydrogen 
bond length, or the intermolecular separation between the bonded 
species (43,83). The question of the effect on the wavenumbers 
and intensities of the vibrations of a hydrogen bonded dimer as 
the distance of separation of the monomer units is changed may 
be addressed by performing calculations at a series of O------ fe) 
interatomic distancesR, distributed about the theoretical 
equilibrium value of 2.832 A°. 

Hence, calculations were done for the linear structure as 
the interaction distance between the two molecules was varied 


over a region -0.30 A° to +0.30 A° relative to the theoretical 
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equilibrium interatomic distance. The intramolecular param- 
eters were left constant and equal to the parameters in the 
theoretical equilibrium geometry for the dimer. The only 
parameter that was allowed to vary was the distance between 
the two molecules. The entire calculation was repeated at 
each new interaction distance, encompassing changes of -0.30 
A°, -0.15 A° + 0.15 A°, and +0.30° relative to the theoret- 
ical equilibrium value of R(O---0O) = 2.832 A°. The predicted 
frequencies and intensities for the linear dimer at these 
varying intermolecular distances are reported in Table 4.11. 
The low frequency modes are reported here also, but will not 
be discussed. Just let us say that upon inspection of the 
predicted frequencies and intensities for these modes it be- 
comes apparent that these parameters are strongly dependent 
upon the assumed intermolecular distance. 

The high frequency intramolecular vibrations are easily 
assigned using the program VIBRATE (see Appendix A). The 
changes in both intensity and frequency as the interatomic 
distance is varied seem to follow logical trends. All three 
frequencies for the electron donor are relatively constant, 
being within 10 cm of each other across the entire series. 
In addition, it appears that as the intermolecular distance 
increases, these frequencies approach the values calculated 
for the isolated electron donor given in column six of Table 
4.10. The intensities vary a little bit more, with a varia- 
tion of 40 km/mole from the highest value (at the shortest 


intermolecular distance) to the lowest value (at the longest 


104 


Aouenbazrz Axeutbewy | 


99 TST vT c8t OVE €OT Ar We 9E~ == e el 
88 8C~ v6T 607 8°T 602 CVE £9 = e TT 
66T OLE OO€ cOE S62 OT? Dre 8VE LT¢ O07? OT 
0D £0 cET S8E Bk: VTE 68 9GE 6PT COE 6 
86E 89S USE 667 FEE 6L0 8TE £67 €bc £90 8 
£6T BEL 6672 OcL OTE Ses DOE SS8 SSE S68 L 
OFT CTLT S€T 6TLT Let STLT DoT LELt 8eT Sec 9 
8rT cO8T est TO8T SST OO8T cLt c6LT O8T LOLT S 
LET 688E S6T O88Ee 6L7 ELBE €1b S98E LSs9 €VBE v 
ST TOOP L°6 L86E LT C66E ST O66E LT €B6E € 
L6 6800 cot O800 OTT 9LOP 8eT TLOp 8ET LSOD c 
96 VETP cot vCIP TUE Lev 6IT ZAG TEL LTtp T 
I Ns 5 ia 66 anaes aac cee cere cee Se 
sTow/wy = a ToOUl/ury wie dToOu/wy rd asTow/wy Fes oTou/uy as pow 
oW CET €=4 oY 786°c=a oY CEB°Z= oY €89°C=4 oY CEG° C= 


“petzea st ‘(O---0)Y Se0ueySTp oTUOWe 
—193UT 9Yy} Se TOWTp AesUTT |aYW AOFZ SeTITSUSUT pue SToquNUsAeM pezOTperd TI°p eTqeyL, 


105 


intermolecular distance) for the asymmetric stretch. However, 
the changes predicted in the intensities as the intermolecular 
distance is increased are towards the values calculated for 
the isolated electron donor. This indicates that as the inter- 
molecular distance is lengthened, the electron donor behaves 
more like an isolated electron donor molecule, and takes on 
the properties of the isolated molecule. This is what should 
be expected based on the fact that as the distance between the 
two molecules is lengthened, the effect of the intermolecular 
interaction on the spectral properties of the molecule should 
be reduced. In summary, the predicted frequencies and inten- 
sities for the electron donor of the linear dimer do not change 
much at interaction distances shorter than equilibrium. At 
longer interaction distances, the predicted frequencies and 
intensities approach the isolated electron donor values. 
However, as the intermolecular distance is varied, the 
frequencies and intensities for the electron acceptor show 
considerable changes (see Table 4.11). These changes are 
primarily isolated in the symmetric stretching mode, and 
therefore we will limit our discussion to this vibration. 
The frequency and intensity for this mode were plotted against 
the interaction distance and this plot is shown in Figure 4.1. 
Both parameters exhibited some dependence on the interaction 
distance, but the nature of the dependence of the frequency 
and nature of the dependence of the intensity were quite 
different. The frequency exhibited a linear dependence, 


while the intensity fit equally well to an exponential and 
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quadratic function. The exponential fit was chosen over the 
quadratic fit because the quadratic function does not cor- 
rectly describe the dependence of the intensity at very long 
intermolecular distances. As the two molecules are infinitely 
separated, it is expected that the intensities of their vibra- 
tions would approach a limiting value, namely the value for 
the isolated molecules. The quadratic function would predict 
that at long intermolecular distances the intensities of the 
isolated molecules would begin to rise again, and eventually, 
at distances long enough, would reapproach the values calcu- 
lated at short intermolecular distances. On the other hand, 
the exponential function predicts that the intensities would 
approach a limiting value at extremely long intermolecular 
distances, and is therefore preferred over the quadratic fit. 
The plot in Figure 4.1 indicates that both absorption param- 
eters for the symmetric stretch of the electron acceptor are 
dependent upon the intermolecular interaction distance, but 
that the intensity is much more sensitive to the interaction 
distance than the frequency. The intensity exhibits larger 
and more drastic changes than the frequency for the same change 
in interaction distance. This suggests that the intensity 
changes and the frequency shifts as the interaction distance 
is varied have different origins. ‘This was analyzed by inter- 
preting both the intensity and frequency parameters as a 
function of R. 

The intensity parameters were transformed to APTs for 


each case. Using the arguments given above, our attention 
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was focused on the APT of the bonded hydrogen. It was found 
that this atom contained the only significant change in its 

APT elements relative to the same terms in the isolated monomer. 
This change was isolated in the Py term in all cases, similar 
to the theoretical geometry described above. By using the CCFO 
analysis most of the predicted change for the Poy term, in each 
case, can be attributed to an increase in charge flux relative 
to the isolated monomer, and not to change in either atomic 
charge or overlap. The total APT and the constituent tensors 
for the bonded hydrogen are listed for all the interaction 
distances calculated, along with the corresponding invariants, 
in Table 4.12. 

The effective invariants, described in Chapter 2, for the 
constituent CCFO tensors are useful in analyzing the various 
contributions to the effective invariant of the total APT 
(effective charge). According to the sum rule (see equation 
54), the square of the effective charge is directly proportional 
to the intensity sum for a molecule. If we recall the data pre- 
sented earlier, the APTs for all of the atoms in the iinear 
dimer with the exception of the bonded hydrogen were essentially 
unchanged from the isolated monomer tensors. Therefore, we 
believe that the changes in intensity predicted for the linear 
dimer relative to the isolated monomer, as the interaction 
distance is varied, can be monitored by observing the changes 
in the effective charge for the bonded hydrogen atom. In turn, 
the relative contributions from the CCFO tensors to the effec- 
tive charge can be visualized by using the effective invariants 


for the CCFO tensors. 
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The effective charge, the effective atomic charge, the 
effective charge flux, and the effective overlap were calcu- 
lated for the bonded hydrogen atom at each interaction 
distance. The values of these invariants for all of the com- 
ponent tensors as well as the total APT for the bonded hydro- 
gen atom at all of the interaction distances calculated are 
listed in Table 4.12. These values were then plotted as a 
function of interatomic distance in Figure 4.2. The effective 
charge exhibited an exponential dependence on the interaction 
distance, similar in form to the dependence that the calcu- 
lated integrated intensity for the symmetric stretch of the 
electron acceptor exhibited with these same changes in inter- 
molecular distance (see Figure 4.1). While the effective 
charge flux also showed an exponential dependence on the in- 
teraction distance, the effective atomic charge and the 
effective overlap both exhibited a linear dependence. In 
addition, the slopes of these lines were very small, indicat- 
ing that the effective atomic charge and the effective overlap 
changed very little as the interaction distance was varied. 
This further indicates that the contributions from the effec- 
tive atomic charge and the effective overlap to the change 
in effective charge as the intermolecular distance is varied 
are also minimal. However, the fact that the effective charge 
flux followed an exponential dependence on this distance indi- 
cates that it is primarily responsible for the predicted in- 
crease in effective charge for the bonded hydrogen. Therefore, 


we conclude that the charge flux tensor. is mainly responsible 
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for the predicted intensity increase in the linear dimer, at 
distances both shorter and longer than the equilibrium distance. 
As shown in Figure 4.1, the frequency of the symmetric 
stretching mode of the electron acceptor exhibited a linear 
dependence on the interaction distance. The frequency param- 
eters were transformed to internal force constants via equation 
33 in order to analyze exactly the origin of this linear 
dependence. The force constants in md/A° for the four intra- 
molecular stretching coordinates, the two intramolecular bend- 
ing coordinates and the hydrogen bond stretching coordinate 
are given, for all the interaction distances calculated, in 
Table 4.13. There appears to be no obvious dependence on a 
single force constant on the interatomic distance. [In fact, 
while some of the expected trends in the force constants are 
found, others are not. For example, it would be expected 
that the hydrogen bond stretching constant F(R 5R.) would de- 
crease as the interaction distance increased. This is not 
the case as indicated by Table 4.13. However, we see in 
Table 4.13 that this force constant does increase as the 
intermolecular distance decreases from the equilibrium value. 
Hence, the expected trend in the hydrogen bond force constant 
is observed for shorter intermolecular distances but not for 
longer distances. However, the predicted changes in the 
hydrogen bond stretching constant are certainly not of the 
same order of magnitude at shorter distances compared with 
those at longer distances. At shorter distances there is 


a definite increase in this force constant; in fact it is 
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almost eight times greater (or an increase of almost 0.600 
md/A°) for an increment in R(O---0O) of only -0.15 A°. At 
longer intermolecular distances this increase is not as pro- 
nounced, being less than 0.100 md/A° for an increment of 
0.15 A° in R(O---0O). We believe that the predicted increase 
in the hydrogen bond force constant as the interaction dis- 
tance becomes longer is actually within the limits of error 
for these types of calculations, and that if this distance 
were lengthened even further, then this force constant would 
indeed actually decrease from the theoretical equilibrium 
value. 

It is interesting that even though the errors inherent 
in the force constant calculations create somewhat unexpected 
trends in the predicted force constants, the calculated fre- 
quencies (for the symmetric stretch of the electron acceptor) 
do seem to follow the expected trends. This suggests that 
there is another effect present which is of greater importance 
than the limits of error in the force constant calculations. 
This effect is the intermolecular coupling that is present 
between the two molecules as they are allowed to interact. 
Using the program VIBRATE (see Appendix A), the symmetric 
stretching vibration for the electron acceptor at three dif- 
ferent intermolecular distances R were frozen and dumped to 
the printer. These vibrations are pictured in Figures 4.3, 
4.4, and 4.5 for distances of R = 2.532 A°, 2.832 A°, and 
3.132 A° respectively. At R = 2.832 A°, or the theoretical 


equilibrium, the symmetric stretching vibration for the 
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electron acceptor is predominately of electron acceptor char- 
acter, with only a small amount of mixing with the motion of 
the two hydrogen atoms of the electron donor. Essentially this 
vibration is localized on the electron acceptor, and there is 
very little coupling between the two molecules. For R = 

3.132 A°, the symmetric stretching vibrations of the electron 
acceptor and the electron donor are still essentially uncoupled. 
The slight increase in the frequency of this vibration relative 
to the frequency of this same vibration at R = 2.832 A° is due 
instead to the slight increase in the O-H stretching force con- 
stant as seen in Table 4.13. This increase in both the force 
constant and the frequency would be expected based on the argu- 
ment that as the molecules are separated further, they should 
behave more like independent monomers, and the corresponding 
properties of these molecules should approach those of the 
isolated molecules. Therefore the bonded O-H stretching force 
constant increases, and is tending towards the monomer value 

of 9.14 md/A°®. Finally, at R = 2.532 A°, the frequency of 

the symmetric stretching vibration of the electron acceptor 

is slightly less than the frequency of the corresponding 
vibration at R = 2.832 A°. The vibration for R = 2.532 A° 

is pictured in Figure 4.3, Here, as evidenced by the slight 
increase in the amplitudes of motion for the hydrogen atoms 

in the electron donor relative to these same amplitudes of 
motion for this same vibration pictured in Figure 4.4 for 

R = 2.832 A°, we can see that the intermolecular coupling 


becomes more important as the distance is shortened. 
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This increase in coupling results in a further decrease in 
frequency (relative to R = 2.832 A°) due to the decrease 
in local O-H oscillator stretching character of this 
vibration. 

Physically, these results appear to make much sense. 
As the two molecules are separated, they would be expected 
to begin to regain the properties of the totally separated 
and isolated molecules. This fact is evidenced by the in- 
crease in the frequency of the symmetric stretching vibration 
of the electron acceptor and the increase in the O-H stretch- 
ing force constant, both towards the monomer values. At 
shorter distances we would expect that the coupling between 
the two molecules would become more important. This is sup- 
ported by the increased coupling present between the electron 
acceptor and the electron donor in the symmetric stretching 


mode of the electron acceptor, as depicted in Figure 4.3. 


CHAPTER 5 
THE CYCLIC WATER-WATER COMPLEX 

The cyclic configuration of a pair of interacting water 
molecules was believed by some workers (20,27) to be the cor- 
rect interpretation of the experimental spectrum for the water 
dimer. Behrens and Luck (27) argue that the available exper- 
imental and theoretical data do not exclude the presence of 
cyclic dimers in matrices, and that there may be both linear 
and cyclic structures present. 

In an ab initio SCF calculation using a 4-31G basis set, 
Kerns and Allen (40) have determined the stabilization energy 
of the cyclic complex to be 7.22 kcal/mole. This value can 
be compared with the value of 8.2 kcal/mole calculated by 
Curtiss and Pople (30), which is in exact agreement with the 
value calculated by Kerns and Allen (40) for the linear dimer. 
While the calculations of the binding energy for the cyclic 
complex by Kerns and Allen (40) was carried out using only 
a restricted geometry optimization (intermolecular parameters 
R and 6 in Figure 1.2), the predicted order of stability 
(linear is more stable than cyclic) is in agreement with 
other theoretical predictions (31,35-37). However, this 
increased stabilization predicted for the linear form over 
the cyclic form is not more than 1.0 kcal/mole. Therefore, 


is it of interest to calculate the infrared spectrum 
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expected for the cyclic structure to compare these results 
with those predicted for the linear form in case both forms 
might be present in the matrix. 

The assumed equilibrium geometry of the pair of water 
molecules in a cyclic configuration is given in Figure 1.2. 
The internal geometries of the two water molecules were 
chosen to be the same as for the monomer at its calculated 
geometry, and the intermolecular parameters were taken from 
the SCF calculation by Morokuma and Pedersen (36). In their 
calculation Morokuma and Pedersen (36) have used a smaller 
basis set than the 4-31G basis set (3 s-type orbitals on H, 
5 s-type with 3 p-type orbitals on 0) and chose to optimize 
the cyclic structure with respect to only three intermolecular 
parameters R, 6 and 4 (see Figure 1.2). Therefore, our as- 
sumed structure does not have the totally optimize geometry 
for this configuration that we might calculate using the 4-31G 
basis set, but it is believed that the optimized geometry would 
not differ by much. The calculation of Morokuma and Pedersen 
(36) on the linear structure gave a distance of R = 2.68 A?®. 
This distance can be compared to the value calculated by Curtiss 
and Pople (30) for the linear structure using the 4-31G basis 
set of R = 2.832 A°. The agreement between these two results 
is fair, and suggests that the parameters calculated by Morokuma 
and Pedersen (36) for the cyclic configuration are a good esti- 
mate of the geometry that we might calculate with the 4-31G 
basis set. However, it should be emphasized here that the 
intermolecular distance (R(O---O) = 2.74 A°) calculated by 


Morokuma and Pedersen (36) for the cyclic complex appears 
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to be shorter than the value that we might predict using the 
4-31G basis set. 

The cyclic configuration pictured in Figure 1.2 has C; 
symmetry, with the center of inversion marked by the asterisk. 
Both molecules are involved in electron donor-acceptor action, 
and therefore neither one can be classified strictly as one or 
the other. For this discussion, we will examine the calcula- 
tions done at the optimized intermolecular distance of Morokuma 
and Pedersen (36) of R = 2.74 A°, where R is defined as the 
distance between heavy atoms (in this case O atoms). This 
distance is probably shorter than the value that would be pre- 
dicted using the 4-31G basis set and we would expect that this 
will overemphasize the importance of any interaction which may 
be indicated by the calculation. The principal cartesian coor- 
dinates for the cyclic configuration at this intermolecular 
distance are given in Table 5.1. The 3N = 18 degrees of free- 
dom for this complex are 9A, + 9A.) and the 3N-6 or 12 vibra- 
tional degrees of freedom are Sh... + 6A. The 6 symmetric A, 
vibrations are not infrared active, but only Raman active, 
while the 6 asymmetric By vibrations are infrared active, 
but not Raman active. The net result is that three high 
frequency modes and three low frequency modes are infrared 
active, while corresponding numbers of modes are Raman active. 

The frequencies and intensities predicted using the pro- 
cedure described earlier for the cyclic pair of water molecules 
at an intermolecular distance of 2.74 A° are listed in Table 


5.2. The corresponding assignments were made with the aid of 
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Table 5.1 Principal cartesian coordinates for the assumed 
equilibrium geometry of the cyclic complex (A°). Atom number- 
ing as given in Figure 1.2. 


X Y Z 
07 -0204 -.0606 -1.3692 
Ho -4294 - 4603 -.6855 
03 -.0204 - 0606 1.3692 
Hy -.4294 -. 4603 -6855 
He -.3590 - 4902 -2.0438 
H - 3590 -.4902 2.0438 
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the graphics program VIBRATE (see Appendix A). The negative 
eigenvalues arise due to the fact that the assumed geometry 
was not the optimized equilibrium positions. From our earlier 
calculations on the linear dimer we would not expect the low 
frequency modes to exert much influence (certainly almost no 
influence on the stretching modes, with only a slight influ- 
ence expected on the bending modes) on the high frequency 
modes due to the large difference in energy. In addition, 

it has been our experience from earlier calculations on iso- 
lated water monomers that the intensities are almost indepen- 
dent of the assumed geometry. 

The essential characteristics of the high frequency vibra- 
tional modes are such that the individual modes of the two mol- 
ecules are coupled either in phase (A) or out-of-phase (A) 
to form the normal mode for the complex. One of these vibra- 
tions, the infrared active symmetric stretching mode, is 
pictured in Figure 5.1. This figure depicts the magnitude 
and phase of the coupling present in this vibration, both of 
which are a direct consequence of the symmetry of the molecule. 
From this figure it is seen that the vibrational displacements 
contain equal contributions from both molecules, and that the 
phase of the coupling between the molecules determines the 
activity of the vibration in the infrared spectrum. In this 
case, the two symmetric stretching vibrations are coupled 
out-of-phase, and therefore this vibration is infrared active. 
Similar descriptions of the remaining vibrations are given in 


Table 5.2. 
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Table 5.2 Predicted wavenumbers and intensities for the 
assumed equilibrium geometry of the cyclic complex at 
R(O---0O) = 2.74 A°, compared to the corresponding pre- 
dicted wavenumbers and intensities for the monomer. 


Complex Monomer 
Mode Description vies > A(km/mole) Sons} A (km/mole) 
A Out of phase 4115 260 4104 54 
asymmetric 
stretch 
2 In phase asym- 4115 0° 4104 54 
metric stretch 
3 Out of phase 3957 46 3983 369 
stretch 
4 In phase 3929 o* 3983 3.9 
symmetric 
stretch 
5 In phase bending 1771 o* 1676 125 
6 Out of phase 1740 270 1676 125 
bending 
7 Out of phase 776 , 756 
twisting of both 
molecules 
8 In phase OH. and 339 o* 
04H. wagging 
9 Ring breathing 208 o* 
10 b 
1l b 
12 b 


“Zero by symmetry. 


Dimaginary frequency. 
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Figure 5.1 Pictorial representation for the infrared active 
symmetric stretching vibration of the cyclic configuration 
of a pair of water molecules with an assumed intermolecular 
distance of 2.74 A°. The solid circles represent the equi- 
librium positions of the atoms while the dotted circles are 
the positions of the atom during one phase of the vibration. 
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The frequencies for these high frequency modes are pre- 
dicted not to be shifted very much from their corresponding 
monomer values even at this very close interaction distance. 
The vibrations corresponding to Ay and AL asymmetric stretch- 
ing remain almost constant, while the symmetric stretching 


4 and 54 cm7} towards lower energy for the 


modes shift 26 cm 
infrared and Raman active modes respectively. These predicted 
shifts are considerably less than the value predicted for the 

shift of the symmetric stretch of the electron acceptor in the 
linear case. The bending modes are predicted to shift upwards 


1 for the infrared active mode and 95 em”! For the Raman 


64 cm 
active band, again less than the prediction for the correspond- 
ing modes of the dimer in the linear configuration. The pre- 
dicted intensity pattern is very interesting, since each of 

the three infrared active bands are predicted to be more in- 
tense than the corresponding bands in the monomer. The 
asymmetric stretch increases by almost a factor of four from 

54 km/mole to 260 km/mole, the symmetric stretch by a factor 

of 11, from 3.9 km/mole to 46 km/mole, and the bending mode 

by a factor of two from 125 km/mole to 270 km/mole. The 
predicted changes for these individual bands are rather large, 
but nevertheless still much smaller than the predicted changes 
for the same bands in the linear configuration. Let us just 
point out here that the intensity sum for these three infrared 
active bands is 576 km/mole, and this sum is comparable to 


twice the intensity sum for an isolated water molecule 


(366 km/mole). 
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Sak Intensity Analysis for the 
Cyclic Complex 


Before beginning, it is necessary to determine exactly 
which effects are responsible for the predicted intensity 
changes; is it the change of the normal coordinates (mech- 
anical effect) or the electronic effects of the interaction 
which would be reflected in the APTs? In the linear case we 
saw that the mechanical effect was very small, while the elec- 
tronic effects were mainly responsible for the predicted 
intensity changes. In order to determine if the same situ- 
ation is true for the cyclic configuration, the program PXPQ 
(see Appendix A) was used. The intensities calculated with 
the normal coordinate transformation for the cyclic configura- 
tion and the corresponding calculated monomer APTs are given 
in column three of Table 5.3. These intensities reflect the 
values that would be expected from two water molecules that 
were not allowed to interact electronically, but were only 
allowed to mix their vibrations to form new normal vibrations. 
It can be seen that the changes in the normal coordinates upon 
complex formation are responsible for almost 93% of the total 
increase in the bending intensity, 42% of the total increase 
in the asymmetric stretching intensity, and 18% of the total 
increase in the symmetric stretching intensity. These results 
indicate that some of the predicted changes in the band inten- 
sities in the cyclic configuration (relative to the isolated 
monomer) are due to a mechanical effect. This result was not 


found to be true in the linear case. However, this can be 
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Table 5.3 Calculated contributions to the predicted intensi- 
ties of the six intramolecular vibrations for the assumed 
equilibrium geometry of the cyclic complex. 


Intensities (km/mole) 


v(em™) Calculated Calculated ay “ Add Monomer 
Cyclic sas amass 2 P(er) pay * pd) 
ete ert en enna re kiutiownhaste sated au Bee EL 
4115 260 110 184 222 54 
4115 o* - - 
3957 46 833 24 35 349 
3929 o* - sé 
1771 0? 
1740 270 250 253 255 125 


"Zero by symmetry. 
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partially explained by remembering that the normal coordinates 
for the cyclic configuration are simply the in phase and out-of- 
phase coupled vibrations of the individual monomer units. Based 
only on these normal coordinates, this would lead us to believe 
that the intensity for each of the three infrared active bands 
would be approximately twice that of the corresponding monomer 
value. The intensities calculated taking only this effect into 
account are almost in exact agreement with twice the intensity 
for the corresponding band in the monomer (see columns three 
and six of Table 5.3). The remainder of the predicted inten- 
sity increase must be due to some electronic effects of the 
interaction between the two molecules, and its origin must lie 
in the polar tensors. 

The predicted intensities for the cyclic configuration 
were then analyzed by transforming the intensity parameters 
to APTs using equations. 20 and 22. The APTs for three of the 
atoms of the cyclic configuration are listed in Table 5.4 in 
the coordinate system shown in Figure 1.2, along with the cor- 
responding monomer APTs rotated to the same coordinate system 
for comparison. Only three of the six atoms are reported here 
because the center of inversion of the molecule requires that 
each atom have an equivalent partner, and therefore the APTs 
for each pair of atoms must be exactly the same (i.e., 01, O33 
H 


Since the normal coordinates are mass weighted motions of 


9.4 


all the atoms in the molecule, the hydrogen atoms are the most 


important atoms to consider in analyzing the intensities. The 
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Table 5.4 Calculated APTs for the assumed equilibrium geometry 
of the cyclic complex compared to ghe corresponding APTs for 
the mononer (e) (le = 1.602 x 10-19 c). Atom numbering and 
coordinate system as in Figure 1.2. 


Cyclic Monomer 
-.7464 -.0642 22257 -.686 -.034 .238 
(@) 
af 
Px, = +.0177 -.6635 -.1188 -.034 -.517 -.029 
- 2006 -.1022 -.8411 .238 -.029 -.760 
-4094 -0406 -.0495 -407 .032 -.064 
H2 
r= .0470 .2926 .1361 .068 .148 -058 
-.0705 .0760 .5137 -.064 .027 -427 
«3370 .0236 -.1762 280 .002 -.173 
Hs 
a -.0647 »371L0 - = 5.174 -.034 369 -.029 


-.1301 -0227 ~3274 =-.173 -002 +333 
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hydrogen atoms Ho and Hy exibit the most significant changes 
in their APT elements. Some changes are predicted in the APTs 
for these atoms; the largest change occurs in the Py term, 
where there is a predicted increase of 0.1446 e or almost a 
98% increase. The predicted changes are relatively small 
compared to the changes predicted for the linear configuration. 
The APTs of the remaining hydrogen atoms are not predicted to 
change much from the corresponding polar tensors in the monomer. 
The CCFO analysis was performed on the APTs for the cyclic 
configuration, as described earlier. The atomic charge, charge 
flux and overlap tensors for the three non-equivalent atoms of 
this geometry are given in Tables 5.5, 5.6, and 5.7 respectively. 
The CCFO analysis is helpful in determining which effects are 
primarily responsible for the increase in the APT element, and 
hence in ihe intensity. For the atomic charge tensor and the 


2 


from the monomer to the cyclic pair. However, the charge flux 


overlap tensor, the Py term of H, is almost constant on.going 


term for this atom increases by 0.1018 e or 70% of the total 
change in the APT element. It becomes clear then that the 
main contribution to the increase in the APT element is due 
to an increase in the charge flux term, and not due to a 
change in either the atomic charge or overlap terms. 

Column four of Table 5.3 gives the calculated intensities 
for the cyclic configuration using the calculated normal coor- 
dinates for this complex and the calculated monomer APTS, 
except now we have adjusted the Poy term of the monomer APTs 


to reflect the predicted increase in charge flux in this term 
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Table 5.5 Calculated charge tensors for the assumed equilibrium 
geometry of the cyclic complex compared to the corresponding 
tensors in isolated water (e) (le= 1.602 x 10-19 c). Atom 
numbering and coordinate system as in Figure 1.2. 


Cyclic Monomer 
-.8381 0 0 -.8034 0 0 
07 
P(C) = 0 -.8381 0 0 -.8034 0 
0 0 -.8381 0 0 -.8034 
-4312 0 0 -4017 0 0 
Ho 
P(C) = 0 ~4312 0 0 -4017 0 
0 0 - 4312 oD 0 4017 , 
-4069 0 0 -4017 0 0 
Hs 
P(C) = 0 - 4069 0 0 -4017 0 


0 0 - 4069 0 0 ~4017 | 


Table 5.6 
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Calculated charge flux tensors for the assumed 


equilibrium geometry of the cyclic complex compared to the 
corresponding tensors in the monomer (e) 


10-19 ¢). 
in Figure 1.2. 


0; -0102 

P(CF)= | -.0634 
-.0132 

H, -.0551 

P(CF)= -0223 
-.0428 

H, -0449 

P(CF)= ~0411 


-0559 


Cyclic 


-.0902 
-. 1583 
-.1361 


-0050 
- 1908 
«O57T7 


-0852 
-.0325 


-0803 


-.0265 016 
-.1648 -.070 
-.0622 -014 
= 0173 -.046 
-1442 |: 041 
-0458 =, 039 
-0437 - 030 
-0206 -030 
-0164 026 


(le = 1.602 x 
Atom numbering and coordinate system as given 


Monomer 
-.070 . OL3 
-.046 -.060 
-.060 SOREL 
-.009 -.039 
-089 -035 
-.008 -.034 
-080 .026 
-.043 025 
-069 .022 
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Table 5.7 Calculated Overlap tensors for the assumed equi- 
librium geometry of the cyclic complex compared to the 
corresponding tensors in the monomer (e) (le = 1.602 x 
10-19 C). Atom numbering and coordinate system as given 

in Figure 1.2. 


Cyclic Monomer 

0 -082 -026 <253 ye ee -035 w22e 
Hl 

P(O) = -O81 ~332 -046 -035 -329 -030 

<2L4 -036 059 ~224 -030 031 

H -033 -036 -.033 -052 -041 -.024 
2 

P(O) = -0258 -.329 -.008 S027 -.343 023 

-.027 021 036 -.024 «035 -060 

H -.115 -.062 -.220 -.153 -.077 -.200 
5 

P(O) =] -.106 -.003 -.038 -.065 01:3 -.055 
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upon complex formation. These calculated intensities indicate 
that the increase in charge flux in the Py term of atom Ho 
is responsible for 50% and 42% of the intensity increase above 
the mechanical effect, for the asymmetric and symmetric stretch- 
ing modes respectively. In addition, these two effects com- 
bined also account for 71% and 52% of the total intensity in- 
crease for these two modes. Column five of Table 5.3 gives 

the intensities calculated for the cyclic configuration by 
adding the differences from the complex to the monomer in 

both charge and overlap. The differences in the charge and 
Overlap tensors are small, and it is seen that their influence 
on the intensities is also smal. In fact, these intensities 

are not much different from those in column four, where we 

have only included the change in charge flux and the mechani- 
cal effect. The intensities in column five are not quite 

equal to the predicted intensities for the cyclic configura- 
tion (column two) because of some minor contributions from 

other APT elements besides the Py term. 

To summarize, the results in Table 5.3 indicate that the 
predicted intensity increase for the cyclic pair of water mole- 
cules has two contributions. Without any changes in electrical 
properties due to the interaction, the mechanical coupling of 
the vibrations of the two individual molecules to form the 
normal coordinates of the complex is responsible for almost 
half of the predicted intensity changes. However, the change 


in electrical properties from the charge flux contributions 
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the Poy terms in the APT of H are also responsible for 
approximately an equal amount for the stretching modes. 


The contributions from charge and overlap are minor. 


5.2 Frequency Analysis for the Cyclic 
Complex 


The predicted frequencies for the cyclic configuration 
are in general not perturbed much from the corresponding fre- 
quencies in the isolated monomer. The frequencies correspond- 
ing to the stretching vibrations of the complex are all within 
65 wavenumbers of the corresponding monomer values. However, 
the bending frequencies are in a bit worse agreement with the 
monomer values, being predicted about 100 cm} higher in the 
cyclic complex. These differences can be analyzed by calcu- 
lating the frequency parameters in internal coordinates. 

The transformation from force constants in S$ coordinates 
to internal force constants is given by equations 18 and 33. 
This transformation requires the A matrix, or the "inverse" 
of the B matrix, derived using equation 32. This requires 
that the G matrix, defined in equation 29, be non-singular. 
The G matrix is non-singular only when a complete non- 
redundant set of internal coordinates have been defined. 

For cyclic molecules it is necessary to define 3N "typical" 
internal coordinates to completely define the kinetic energy 
of the molecule (10,11). This introduces six cyclic redun- 
dancies, which must first be removed before the A matrix can 
be calculated. Various workers have discussed this problem 


in the past (84-90), and have proposed different solutions, 
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any one of which can be used. In essence, the solution to 
this problem involves choosing the proper combinations of 
internal coordinates such that these combinations are exactly 
the redundant coordinate. These redundant coordinates can 
then be ignored, provided the remaining vibrational coordin- 
ates are described by taking symmetrized linear combinations 
of these internal coordinates such that these combinations 
are orthogonal to the redundant coordinates. In this way we 
can reduce the number of original internal coordinates (by 
the number of redundancies) to the appropriate number of 
symmetry coordinates (usually 3N-6). More detailed discus- 
sions concerning the actual procedures involved can be found 
elsewhere (84-90). 

After reducing the number of internal coordinates (by 
the number of redundancies) in order to form the symmetry 
coordinates, the potential energy in terms of these symmetry 
coordinates can be uniquely determined. However, it can be 
shown that the transformation of the potential from the set 
of non redundancy symmetry coordinates to the original redun- 
dant set of internal coordinates is not unique (see Appendix 
C, equation 114). Instead, the force constants in the origi- 
nal set of redundant internal coordinates can only be deter- 
mined by making some assumption concerning the form of the 
potential in these internal coordinates. In these calcula- 
tions we have assumed that the force constants corresponding 
to the interaction between the intramolecular and the inter- 


molecular vibrations are all equal to zero. We have seen 
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from previous results that this approximation is almost true, 
but not entirely correct. The vibrations corresponding to 

the bending motions were shown to mix with the low frequency 
intermolecular vibrations. Therefore, using this assumed 
potential, the internal bending force constants are determined 
as "effective" bending constants ("effective" with respect to 
the fact that they now include the interactions with the low 
frequency vibrations and are no longer pure bending force 
constants), and are not strictly comparable to the bending 
constants calculated in other isolated molecules. We would 
expect that the stretching constants calculated in this way 
are almost pure stretching constants though because these vi- 
brations do not mix with the low frequency intermolecular modes. 

The set of intramolecular internal coordinates used in 
these calculations are defined in Table 5.8, and the corre- 
sponding calculated force constants using the approximation 
previously discussed are given in Table 5.9. The frequencies 
calculated with these force constants are given in Table 5.10, 
column three. These frequencies agree very well with the ones 
given in Table 5.10, column two, which were calculated without 
any assumptions about the form of the potential. 

The intramolecular frequencies given in Table 5.10 can now 
be analyzed in terms of these calculated force constants. The 
diagonal force constants are not much different from the corre- 
sponding monomer values (see Table 3.6), which is to be ex- 
pected since the individual monomer geometries in the cyclic 


configuration were not changed from the isolated monomer 
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Table 5.8 Intramolecular internal coordinate definitions 
used in the cyclic complex. Atom numberingas in Figure 1.2. 


Coordinate No. Definition 
Ri R (0,-H,) 
Ro R (0, -H,) 
R3 R (0,-H,) 
Ry R (0,-Hg) 
Re dé (H,-0,-H,) 
R d6@ (H¢-0,-H,) 
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Table 5.9 Calculated internal "effective" force constants 
(see text) for the six intramolecular internal coordinates 
(md/A°) of the cyclic complex. Internal coordinate defini- 
tions as given in Table 5.8 (weighting is done by d, where 
d== 1 ‘A°). 


Ri Ry R3 Ry R, R¢ 
Ri 9.029 -.1646 -.1264 =.0203 - 2289 -.0143 
Ro 9.123 -.0203 -0602 - 3459 =<5002 
R3 9.029 -.1646 -.0145 - 2289 
Ry 9:..123 -.002 - 3459 
Re - 7743 -0081 
R tas 
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geometry. The bending constant is the only one which is pre- 
dicted to change, increasing slightly by 0.06 md/A° from 0.714 
md/A° in the monomer to 0.774 md/A° in the cyclic complex. How- 
ever, this bending force constant is an "effective" bending con- 
stant and cannot be strictly compared to the bending constant 
calculated in the isolated monomer. However, if we assume that 
this "effective" bending constant is the same in the cyclic com- 
plex as for the monomer (0.714 md/A°) the resulting predicted 
frequencies are given in Table 5.10, column four. Both bending 
frequencies are then predicted to shift down from 1759 cm - and 


: and 1668 cm, close to the expected mon- 


omer value of 1673 cm=+, This suggests that the origin of the 


1739 cm ~ to 1689 cm 
predicted increases in the bending frequencies is due to coup- 
ling between these bending motions and the low frequency inter- 
molecular vibrations. The predicted stretching force constants 
for all four O-H bonds of the cyclic complex (Table 5.9) are not 
very different from the monomer values. The differences are 
small, predicted to be less than 0.5%. These changes do not 
affect the predicted frequencies very much. The differences 
between columns four and five of Table 5.10 indicate the 

effect of these changes, which are less than 0.25%. The change 

in the stretching frequencies upon complex formation can be 
explained by the predicted change in the interaction constant 
between the two stretching coordinates of one molecule. This 
force constant (Fy>5 and P34) is predicted to be -.1646 md/A° 

in the cyclic complex, and only -0.065 md/A° in the monomer. 


The frequencies calculated after adjusting this force constant 
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to its monomer value are listed in Table 5.10, column six. 
These frequencies are now in good agreement with the monomer 
values, the average of the two frequencies being less than 0.53 
different from the monomer value of 3997 cmt, This interaction 
constant determines the amount of splitting between the asym- 
metric and symmetric stretching modes within an individual 
molecule. The large negative value (-0.1646 md/A°) in the 
cyclic complex indicates that in the formation of the complex 
these two modes are split further than they are in the iso- 
lated monomer. The splitting of these modes in the isolated 
monomer is 121 cmt, while these same splittings are 158 cm=2 


and 187 cm7- 


in the cyclic complex. If the value of this 
interaction constant is changed from its predicted value of 
-0.164 md/A° to the predicted monomer value of -0.065 ke 
then the splittings predicted for the cyclic complex are re- 
duced to 110 cm Ne and 138 cm * (see Table 5.10, column Six)-. 
In addition the stretching frequencies predicted with this 
interaction constant are now in better agreement with the 
predicted monomer values. 

The predicted increases in the splitting between the 
stretching modes in the cyclic complex can best be explained 
by referring back to the structure of the cyclic configuration 


given in Figure 1.2. The Symmetry is reduced from c ‘symmetry 


2v 
in the monomer to C; Symmetry in the complex. With this reduc- 
tion in symmetry the two 0-H bonds within each molecule are no 
longer equivalent, and the coupling between these two bonds 


Changes. Asa result, the predicted stretching frequencies 
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are split further in the complex than in the isolated monomer. 
This seems to point to the origins of the predicted frequency 
shifts in the cyclic complex relative to the isolated monomer 
as being related to symmetry arguments. We have shown that 

the predicted shifts in bending frequencies are due to the 
increase in the "effective" bending constant, which in turn 

is probably due to the presence of intermolecular coupling 
between the bending vibrations and the low frequency vibrations. 
In addition, the predicted shifts in the stretching frequencies 
have been shown to be due to the presence of increased intra- 
molecular coupling, reflected in the interaction force con- 
stants RR, and R3R,- These results are completely different 
from the results predicted for the linear dimer, where we saw 
that the frequency shift of the symmetric stretching mode for 
the electron acceptor was almost entirely a result of the 
weakening of the bonded O-H bond upon dimer formation. 

Finally, it is also interesting to note that the intermolecular 
coupling does not seem to play an important role in determinine 
the frequencies of the six intramolecular vibrations. The 
intermolecular coupling constants between the stretching coor- 


dinates of the two different molecules (F is -0.12644 md/A°, 


13) 
a reasonably large amount. However, the effect of changing 

this constant to -0.05 md/A° is minimal. The frequencies cal- 
culated using this value are listed in column six of Table 5.10. 


The changes from the predicted intramolecular frequencies of 


the cyclic complex are insignificant. 


148 

5.3 Further Calculations on the Cyclic Complex 

at Different Interaction Distances 

In the linear dimer it was shown that the predicted spec- 
tral changes upon complex formation were strongly dependent 
upon the assumed geometrical orientation of the two molecules. 
In the cyclic configuration the spectral changes are predicted 
to be much less significant than those predicted for the linear 
dimer. This suggests that for less energetically stable com- 
plexes the dependence of the spectral properties of the inter- 
molecular distance might be expected not to be as strong as in 
more energetically favorable complexes. Therefore, additional 
calculations were performed on the cyclic configurations with 
varying intermolecular distances. The infrared spectrum of 
the cyclic pair was calculated at a shorter intermolecular 
distance (2.65 A°), a longer distance (2.99 A°), and also with 
one molecule rotated 180 degrees so that the center of inver- 
sion was destroyed. The predicted frequencies and intensities 
for these different configurations of the cyclic pair of water 
molecules are given in Table 5.11. The predicted frequencies 
for an intermolecular distance of 2.65 A° are strongly per- 
turbed from the expected monomer values. This is not what 
would be expected for such a weak energetically stabilized 
system. A calculation at a distance of 2.51 A° predicted even 
more highly perturbed results. The large predicted changes in 
frequencies at these shorter intermolecular distances are not 
due to the interaction of the two molecules, but are believed 
to result instead from the presence of extremely strong mech- 


anical effects. The effect of forcing two molecules together 
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is that we no longer have two interacting molecules, but rather 
one unrecognizable entity which retains none of the properties 
of the individual molecules. For example, the predicted high 
frequency normal vibrations for the cyclic configuration at an 
intermolecular distance of 2.51 A° do not even begin to resemble 
the corresponding vibrations in isolated water. This is shown 
in Figure 5.2, where we have pictured the normal mode for the 
vibration which should correspond to the out-of-phase combina- 
tion of the asymmetric stretching modes of the two molecules. 
While the bonded O-H bonds are indeed stretching, the non- 
bonded O-H bonds are wagging. Clearly, the combinations of 

the molecular vibrations that one would expect (for example 

see Figure 5.1 and the descriptions of these vibrations in 
Table 5.2) to form the normal vibrations are not present at 
this intermolecular distance. Instead, the predicted vibra- 
tions are unrecognizable and uninterpretable. For this reason 
the frequencies predicted for the intermolecular distance of 
2.51 A° are not believed to carry much significance and are 

not listed in Table 5.11. 

Therefore, disregarding the configuration without the 
center of symmetry for now, we are left with only two addi- 
tional distances to consider. The predicted frequencies for 
the molecular separation of 2.74 A° and 2.99 A° are listed 
in Table 5.11. The predicted wavenumbers for the cyclic com- 
plex at both of these distances are almost identical, with the 
largest difference being only 12 cm7+ for the Raman active bend- 


ing mode. Because of the near equivalence of the two sets of 
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Figure 5.2 Pictorial representation for the infrared active 


asymmetric stretching vibration of the cyclic configuration 
of a pair of water molecules with an assumed intermolecular 
distance of 2.51 A°. The solid circles represent the equi- 
librium positions of the atoms while the dotted circles are 
the positions of the atoms during one phase of the vibration. 
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predicted frequencies, it is not worthwhile to try to analyze 
these differences based on changes in force constants as the 
distance is changed. It suffices to say here that the absorp- 
tion frequencies for weak energetically stabilized complexes 
are almost independent of the interaction distance (at reason- 
able values for this distance). This statement is further sup- 
ported by a series of calculations performed on the weak van 
der Waals complex H,0--CO at different interaction distances 
(91). 

On the other hand, the predicted intensities for these 
complexes discussed above show more significant changes. In 
general, the predicted intensities for the cyclic configura- 
tion exhibit a slight increase as the interaction distance is 
decreased. These results are similar to those that were found 
earlier for the linear dimer. Even at the interaction distance 
of 2.65 A°, where the predicted frequencies are physically un- 
reasonable, the intensities follow the expected trend, in that 
they are predicted to show. a slight increase over those pre- 
dicted at longer distances (2.74 A° and 2.99 A°). In addition, 
these predicted values seem to be reasonable quantities based 
on the predicted intensities at an intermolecular distance of 
2.74 A°. It appears from these results that the prediction of 
infrared intensities is not as sensitive to the assumed geometry 
as the prediction of infrared absorption frequencies. This may 
be due in part to the fact that the frequencies are directly 
dependent upon the geometry due to the contributions from the 
kinetic energy (G matrix) (see equation 6), while the intensi- 


ties are not directly dependent on the geometry in any way. 
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The predicted changes in intensities in the cyclic complex, 
which are due to electronic effects as the interaction distance 
is changed, can be analyzed by interpreting the predicted APTs 
for this complex at varying interaction distances. We have 
previously shown that the bonded hydrogen atom is the most 
important atom to consider when trying to explain the pre- 
dicted electronic intensity changes in the cyclic complex 
relative to isolated water. Therefore, the APTs for the 
bonded hydrogen atoms at various intermolecular distances 
were calculated and are listed in Table 5.12. The calcula- 
tion at R = 2.65 A° was also included in this table because 
while the predicted frequencies seem unreasonable, the pre- 
dicted intensities are reasonable values for what we might 
expect. The CCFO analysis for these predicted APTs is also 
gives in this table. These tensors are all reported in the 
coordinate system given in Figure 1.2. Upon inspection of the 
APTs, it becomes obvious that there are changes as the inter- 
action distance is varied, but that these changes are not very 
drastic. In addition, when the CCFO tensors are analyzed, it 
appears that the majority of these changes in the APTs are due 
to differences in charge flux. The validity of this statement 
can be tested by analyzing the effective invariants of the APT 
and its constituent tensors, the CCFO tensors. 

The effective invariants for the APT and its constituent 
tensors have been previously defined in Chapter 2. These four 
effective invariants are plotted versus the interaction dis- 


tance in Figure 5.3. From this figure it becomes obvious that 
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the total effective charge exhibits a small rise as the inter- 
action distance is decreased, thereby reflecting the predicted 
slight increase in intensity as the interaction distance is 
shortened (see Table 5.11). Also, the effective atomic charge 
and the effective Overlap remain almost constant with the change 
in interaction distance; in fact, the two linear least squares 
lines drawn through these two sets of points are both almost 
exactly horizontal. The effective charge flux exhibits a 
Slight rise with decreasing interaction distance indicating 
that it is the tensor that is responsible for the changed pre- 
dicted in the behavior of the total APT, 

The predicted frequencies and intensities for the cyclic 
configuration without a center of inversion, and with an inter- 
molecular distance of 2.74 a° are listed in Table 5.11. The 
frequencies are Similar, but significantly different, from 
those predicted for the configuration with the same intermol- 
ecular distance but with the center of inversion. Of course, 
the predicted intensities for the individual bands in these 
two complexes are somewhat different because of the different 
symmetries in these two complexes. In the complex with the 
center of inversion only three intramolecular vibrations are 
infrared active. However, in the complex without the center 
of inversion all six intramolecular vibrations are infrared 
active. Therefore we would expect that the form of the normal 
vibrations would be different in these two complexes, and that 
this would have an effect on the intensity distribution among 


the normal modes (i.e., the mechanical effect previously 
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discussed). If the intensity sums for the two configurations 
are compared they are found to be almost equal (576 km/mole 
with the center of inversion compared to 554 km/mole without 
the center of inversion). This indicates that the predict 
intensities for this complex are indeed distributed differ- 
ently as the orientation of the two molecules is changed, 

but that there is no predicted change in total intensity. 

This emphasizes the importance of the mechanical effect (form 
of the normal coordinates) in determining infrared intensities. 
Finally, this calculation also reinforces the notion that the 
predicted frequencies are very sensitive to the assumed geom- 
etry, while the total predicted intensities are less sensitive 


to these parameters. 


5.4 Comparisons Between Linear and 
Cyclic Complexes 


The cyclic configuration of a pair of interacting water 
molecules is energetically less favorable than the linear 
configuration (40). In addition, the predicted spectra for 
these two different configurations are quite different from 
each other, and also from the predicted spectrum of isolated 
water. 

The predicted frequency shifts in the cyclic complex are 
considerably smaller than those predicted in the linear dimer. 
In addition, the origin of these shifts is different in the 
two complexes. In the linear dimer we have shown that the 
large predicted shift in the symmetric stretching mode of the 


electron acceptor was due primarily to the predicted weakening 
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of the bonded O-H bond. In the cyclic configuration, the pre- 
dicted shifts in the stretching modes are due to increased 
intermolecular coupling. Additionally, while the predicted 
shift in the symmetric stretching mode of the electron accep- 
tor in the linear dimer was dependent upon the assumed inter- 
molecular distance, the frequencies of the stretching modes 
in the cyclic complex were almost independent of this distance. 
In the case of the predicted intensities, the linear dimer 
was predicted to have a large intensity increase, relative to 
the monomer, in the symmetric stretching mode of the electron 
acceptor. This increase was shown to be due almost entirely 
to electronic effects present in the polar tensor of the bonded 
hydrogen atom, specifically isolated to the charge flux tensor. 
In the cyclic complex the predicted intensities exhibited less 
significant differences from the predicted intensities for iso- 
lated water. The origin of these differences were also found 
to be a bit different from the linear dimer. In the cyclic 
configuration the predicted intensity increases were due to 
two effects in almost equal contributions. The first effect 
was responsible for almost half of the predicted intensity 
increases and is called the mechanical effect, or the chang- 
ing of the form of the normal coordinates in the complex rela- 
tive to the monomer. This effect was not found to be of any 
Significance in the linear dimer. The second effect was the 
electronic effect reflected in the polar tensor of the bonded 
hydrogen, specifically in the charge flux tensor. This effect 


is similar to the electronic effect in the linear dimer, but 
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its magnitude and importance were considerably reduced for the 
cyclic configuration. Finally, the electronic effect was shown 
to be rather strongly dependent upon the assumed interaction 
distance in the linear dimer. In the cyclic complex the 
electronic effect exhibited only a very subtle dependence 


on the interaction distance, 


CHAPTER 6 
THE IN-PLANE BIFURCATED WATER-WATER COMPLEX 

The third orientation for the water-water complex which 
has received some attention by past experimental (26) and 
theoretical (34,36,40) work is the in-plane bifurcated struc- 
ture pictured in Figure 1.3. This structure is predicted theo- 
retically (40) to be the least energetically favorable 
(predicted binding energy of 5.96 kcal/mole, compared to 7.22 
kcal/mole and 8.2 kcal/mole in the cyclic and linear forms 
respectively), and has all but been eliminated as a possible 
structure for the water dimer in an inert gas matrix (26). 
However, the experimental conditions in a carbon monoxide 
matrix can be considerably different from the conditions in 
an inert gas matrix. Hagen and Tielens (26) have suggested 
that in a CO matrix the bifurcated structure for the water- 
water complex may be the most stable form. Therefore, from 
the experimental point of view it becomes of interest to pre- 
dict the spectrum for the bifurcated structure. Additionally, 
it is of interest to this study because it offers the Opportu- 
nity to investgate the effect on the predicted spectrum of yet 
another orientation for the water-water complex. 

The intermolecular distance used in these calculations 
was taken from the SCF calculation of Morokuma and Pedersen 


(36), where the same calculational procedure was followed 


Louk 
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as was the case of cyclic complex. The geometries of the 
individual monomer units were fixed at the predicted theoret- 
ical values (R = 0.941 A°, and 6 = 111.2°). The complex was 
assumed to have Coy symmetry, and both molecules were con- 
strained to lie in the same plane. The calculated 0O--0 dis- 
tance of Morokuma and Pedersen (36) was 2.89 A°. The 
calculated principal cartesian coordinates for this 
interaction distance are listed in Table 6.1. 

The irreducible representations for the 3N = 18 degrees 
of freedom for this molecule are 6A, + 2A, + 6B, + 4B, and 


af 1 1 2 


the 3N-6 = 12 vibrational degrees of freedom are 5A, +A 


1 2 


+ 4B, + 2B. 

In this particular orientation, both hydrogens of the 
electron acceptor molecule (Hy and He in Figure 1.3) are able 
to interact with the electron donating oxygen (A, in Figure 
1.3) of the second molecule. These two interaction vectors 


are not orientated along the same direction, but are at an 


angle of 39 degrees to each other for R(O---O) = 2.89 A°. 


6.1 The Predicted Spectrum for the 

Bifurcated Water-Water Complex 

(R = 2.89 A°) 

The calculated wavenumbers and intensities for the bifur- 
cated complex at this interaction distance are given in Table 
6.2, along with the corresponding predicted monomer values. 
The vibrational characteristics of the calculated normal modes 


listed in this table were assigned with the aid of the program 


VIBRATE (see Appendix A). The predicted frequencies for this 
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Table 6.1 Principal cartesian coordinates for the in-plane 
bifurcated water-water complex using an assumed intermolecular 
distance R(O---0O) = 2.89 A°. Atom numbering as given in 
Figure 1.3. (A°) 


xX Y Z 
Hy “1.9272 0 0.7847 
05 #1.3899 0 0 
H3 -1.9272 0 -0.7847 
Hy 0.9728 0 0.7847 
0. 529104 0 0 
H 0.9728 0 -0.7847 
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bifurcated configuration exhibit only minute changes from the 
predicted frequencies for the isolated monomer. The largest 
predicted frequency shift occurs for the bending modes, where 
these modes for both molecules in the complex are predicted to 
increase about 90 cm in frequency from the isolated monomer. 
This predicted increase in the frequencies of the bending 
modes appears to be typical of these types of calculations. 
Approximately the same increase was predicted in the linear 
and cyclic configurations. Since there is a marked difference 
between these three interactions, there does not appear to be 
any direct correlation between the predicted increase in the 
bending frequencies and the nature of the interaction. Instead, 
these shifts are predicted to occur because of the coupling be- 
tween these bending vibrations and the intermolecular vibrations 
of the complex. This effect is present in all three configura- 
tions, and therefore appears to be mechanical in nature. It was 
reflected in the linear dimer by the calculated coupling con- 
stants between the low frequency vibrations and the bending 
vibrations and was seen in the cyclic complex by the predicted 
increase in the effective bending constant. In the bifurcated 
complex we have used an approximate potential similar to the one 
used in the cyclic configuration to describe the intermolecular 
vibrations. This potential contains only the six intramolecular 
coordinates defined in Table 6.3, and does not include any coup- 
ling constants between the intramolecular high frequency and 
intermolecular low frequency vibrations. The force constants 


calculated in this way are listed in Table 6.4. This table 


predicts an increase in the effective bending constant from a 
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Table 6.3 Internal coordinate definitions used for the 
assumed equilibrium geometry of the bifurcated complex 
(R(O---O) = 2.89 A°). Atom numbering as given in Figure 
1.3 (d is a weighting factor where d = 1 A®). 


Coordinate No. Definition 
Ri R (Hy-O.) 
Ro R (H_-O5) 
R3 R (H,-0) 
Ry R (H,-0,) 
Rs dy 8 (H,-O.-H,) 


do 9 (H)-O,-H) 


Table 6.4 


a 


9.062 


Ro 


=<L373 


9.062 
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Calculated "effective" 

(see text Chapter 5) (md/A°) for th 
geometry of the bifurcated 
Internal coordinate def 
(Weighting is done by d 


Re 


-.0029 


-.0051 


9.1208 


complex (R(0---0) 
initions as given in Tab 
where d 


2 Bees 


Ry 


-.0051 
-.0029 
= 25 


9.1208 


internal force constants 
e assumed equilibrium 
289 A°). 

le 6.3. 
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value of 0.714 md/A° in the monomer to the predicted value 

of 0.802 md/A° in the bifurcated complex. When both effec- 
tive bending constants are lowered to the predicted monomer 
value, the calculated bending frequencies in the bifurcated 
complex are 1667 om and 1679 cmt, almost in exact agreement 
with the predicted value of 1676 om7+ for the isolated monomer. 
This suggests that the predicted increase in the bending fre- 
quencies of the bifurcated complex are a result of the pre- 
dicted increases in the effective bending constants. These 
effective force constants are predicted to increase as a re- 
sult of the mechanical intermolecular coupling present in the 
complex. 

In summary, the predicted wavenumbers for the bifurcated 

complex agree very well with the predicted monomer value. 
The small differences are not believed to be Significant, or 
can be explained on the basis of mechanical coupling between 
the bending vibrations and the low frequency intermolecular 
vibrations. 

The predicted integrated intensities for the bifurcated 
configuration are also given in Table 6.2, along with the 
corresponding monomer values. Here again, the intensities 
for the four stretching modes are predicted to be practically 
the same as the corresponding values for the isolated monomer, 
while the bending modes are predicted to be slightly more 
intense than in the monomer. The small differences for the 


stretching modes are not considered to be significant. 
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However, the predicted increase of almost a factor of two in 
the intensity for the bending mode of the electron donor does 
need to be explained. 

The APTs, and the charge, charge flux and overlap tensors 
for the bifurcated complex are listed in Tables 6.5, 6.6, 6.7, 
and 6.8 respectively. The coordinate system used is defined 
in Figure 1.3. The predicted values for the monomer in the 
same coordinate system are also listed for comparison. There 
does not appear to be any significant difference between the 
APTs of the complex and those of the monomer, perhaps indi- 
cating that the predicted changes in the intensity of the 
bending vibrations are not an electronic effect, but rather 
a mechanical effect. The charge, charge flux and overlap 
tensors are also very similar in the two systems. The CCFO 
analysis appears not to offer any reason for the apparent 
increase in intensity for ss bending mode of the electron 
donor. This leaves the explanation with the second contrib- 
uting factor in the calculation of infrared intensities, 
namely changes predicted in the form of the normal coordin- 
ates as the complex is formed, or the mechanical effect. 

The increase in the frequency of the bending vibration 
of the electron donor was explained on the basis of the 
presence of intermolecular coupling between this vibration 
and the low frequency vibrations of the complex. It is 
expected that this coupling affects not only the predicted 
vibrational frequency, but also the form of the normal 


vibration. Any change in the form of the normal vibration 
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Table 6.5 Calculated APTs for the assumed equilibrium geometry 
of the bifurcated complex (R(O---0O) = 2.89 A°) (e) (le = 1.602 


x 107°" Cc). Atom numbering and coordinate system as given in 
Figure 1.3. 
Bifurcated Monomer 
.4969 0 0 .4820 0 0 
Hy 
Px = .1995 -.0638 0 .1482 -.0887 
-.0411 .3434 0 -.0147 .3518 
.9856 0 0 -.9640 ) 0 
fe) 
: 
Py = -.535 +.0010 0 =.5714 .0448 
+.0010 -.5348 0 .0448 -.4826 
.4969 0 0 4820 0 0 
H3 
fy = .3599 +.0211 0 3692 .0440 
-.0016 1830 0 -.0030 1308 
.4900 0 0 
Hy 
ey * .1867 -.0232 
-.0330 23971 
-.9882 0 0 
Os 
ella 0 -.5992 +.0195 
0 +.0195 -.5841 
4900 0 0 
He 
Py = 0 3886 +.0453 
0 +.0551 1953 
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Table 6.6 Calculated charge tensors® for the assumed 
equilibrium geometry of the bifurcated complex (R(0---0) = 
2.89 A°) (e) (le = 1.602 x 107-19c). Atom numbering and 
coordinate system as given in Figure 1.3. 


Bifurcated Monomer 
Hy -4129 - 4017 
05 -.8180 -.8034 
H, 4129 .4017 
Hy - 4097 -4017 
O. -.8272 -.8034 
He - 4097 - 4017 


“The quantities listed here are the calculated Mulliken 
charges, Qar where the charge tensor is simply PA(C) = 
ORs 

A= 
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Table 6.7 Calculated charge flux tensors for the assumed 
equilibrium geometry of the bifurcated complex (R(0---0) = 
2.89 A°) (e) (le = 1.602 x 10-19 @. Atom numbering and 
coordinate system as given in Figure 1.3. 


Bifurcated Monomer 
0 0 0 0 0 0 
Hy 
P(CF)= 0 1279 -.0353 0 .0888 -.0546 
0 -0063 -.1025 0 .0124 -.0788 
05 0 0 0 0 0 0 
P(CF)= 0 -.0392 -0699 0 -.0451 .0932 
0 -0699 0151 0 -0932 -0272 
0 0 0 0 0 0 
H3 
P(CF)= 0 -.0626 -.0675 0 -.0426 -.0386 
0 -.1091 -0880 0 -.1056 -0526 
/ 
0 0 0 
H4 
P(CF)= 0 - 1089 -0263 
0 -0143 -.0418 
0 0 0 
Os 
P(CF)= 0 -.0991 -0484 
0 .0484 -.0616 
0 0 0 
He 
P(CF) = 0 -0358 -.0418 
0 -.0298 - 1028 
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Table 6.8 Calculated overlap tensors for the assumed equi- 
librium geometry of the bifurcated complex (R(O---0) = 

2.89 A°) (e) (le = 1.602 x 10719 c). Atom numbering and 
coordinate system as given in Figure 1.3. 


Bifurcated Monomer 


-.1610 0 0 
°5 
P(O) = 0 3271 -.0289 


.0840 0 0 .0803 0 0 
Hy 
P(O) = 0 =.3413 -.0285 0 =~. 9423- 5. »i0gee 
0 -.0474 .0330 0 -.0541 .0289 
-.1676 0 0 -.1606 0 0 
02 
P(O) = 0 3127 -.0689 0 53) -.0484 
0 -.0689 .2681 0 -.0484 . 2936 
. .0840 0 0 .0803 0 0 
3 
P(O0) = 0 .0096 .0886 0 .0101 .0826 
0 «1695. -=.3179 0 , 1026... .~,3235 
.0803 0 0 
Hy 
P(O) = 0 -.3319 -.0495 
0 ~.0473°~. . 2092 
0 -.0289 “° 
.0803 0 0 
H¢ 
P(O) = 0 .0147 .0871 
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upon complex formation should reflect a corresponding change 
in the predicted intensity for that mode due to the dependence 
of the predicted intensity on the normal coordinates (see 
equation 24). The intensities calculated with the monomer 
APTs and the predicted normal coordinate transformation for 
the bifurcated complex are listed in Table 6.9, along with 
the intensities predicted for the bifurcated pair. The pre- 
dicted intensities using the monomer APTs are in good agree- 
ment with the values predicted for the bifurcated complex. 
This indicates that the change in the form of the normal 
coordinate is responsible for almost all of the predicted 
intensity increase in the bifurcated complex compared to 

the isolated monomer. Finally, it also indicates that in a 
complex that is only slightly energetically stabilized (40) 
the predicted frequencies and intensities should be very sim- 
ilar to the isolated monomer values. Any difference between 
the two sets of values are erobabiy not due to an electronic 


interaction, but rather a mechanical effect. 


6.2 Additional Calculations on the 


Bifurcated Water-Water Complex 


The calculated wavenumbers and intensities for the bi- 
furcated pair of water molecules at different interaction 
distances and with the two molecular planes perpendicular 
instead of parallel are listed in Table 6.10. The predicted 
values listed here for all three cases are not much different 


from each other or from the isolated monomer numbers. The 
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Table 6.9 Calculated mechanical effect (see text) on the pre- 
dicted intensities for the six intramolecular vibrations of 
the assumed equilibrium geometry of the bifurcated complex 
(km/mole). 


Intensity (km/mole) 


Mode Calculated Calculated bifurcated Calculated 
Monomer using monomer APTs bifurcated 
1 54 57 79 
2 54 SL 79 
3 3.9 2.6 6.5 
4 3.9 4.7 14 
5 125 192 234 


6 125 58 LL7 


Lv7 


Table 6.10 Predicted wavenumbers and intensities for the 
six intramolecular vibrations of the bifurcated complex 
as the interaction distance, R(O---0), is varied. 


Mode R = 2.89 A° R = 2.89 aca R = 3.05 ae 


v(cm_ ) A(km/mole) viea>) A(km/mole) ven) A(km/mole) 


1 4119 79 4112 88 4113 72 
2 4102 79 4102 85 4101 81 
3 3960 6.5 3954 6.5 3955 6.1 
4 3949 . 14 3948 18 3948 14 
5 L777 234 1776 306 1779 216 
6 L757 117 1748 51 1756 i123 


a ; 
The electron donor molecule is rotated 90° so that the two molecules no 
longer lie in the same plane, but have their molecular planes now 


pendicular to each other. 
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intensities of the bending modes are the only vibrations that 
are predicted to differ from the monomer intensities. How- 
ever, the sum of the intensities for the two bending modes 
remains almost constant for the three different cases. 
Therefore, the distribution of the total intensity changes 

in each case, but the sum is constant. This indicates that 
the predicted differences in intensity between the bifurca- 
ted complex and the isolated monomer are due to a mechanical 
effect and not an electronic interaction. This mechanical 
effect is the coupling present between the bending modes 

and the low frequency vibrations, and is responsible for 
changing the character of the normal vibrations, and hence 
the distribution of the total intensity. The small differ- 
ences predicted in the frequencies of the stretching modes 
relative to the monomer in all three cases are not considered 
significant. Finally, the predicted shifts in the bending 
vibrations for the bifurcated complex are independent of the 
interaction distance assumed, or of the orientation of the 
two molecules. This is additional evidence that the predicted 
frequency shifts in the bifurcated complex are due to mechan- 
ical coupling and not due to an intermolecular interaction 
between the two molecules. If these predicted shifts were 
due to some type of intermolecular interaction then we would 
expect that they would change as some of the intermolecular 


parameters were varied. 


CHAPTER 7 
THE OM.H,0.0M COMPLEX, M=0, -l, -2 

Until now, we have been concerned with the calculation and 
interpretation of the infrared spectra of one water molecule 
interacting with another water molecule in various configura- 
tions and at different interaction distances. The predicted 
spectra for these species have varied considerably, and depend 
rather strongly on both the assumed geometry and the interac- 
tion distance. Analysis of these predicted spectra has pro- 
videc considerable insight into the mechanical and electrical 
contributions to the predicted changes in intensity upon com- 
plex formation, and also into the influence of factors such 
as geometry, intramolecular coupling and intermolecular cou- 
pling on the predicted frequencies. Naturally, it follows 
that it would be of interest to study the effect of changing 
the strength of the electron donor in these types of calcula- 
tions, and to determine if the effects predicted in the water- 
water complex are dependent upon the physical presence of the 
second water molecule or if they might be a general character- 
istic of intermolecular interactions. 

One way to do this is to investigate the interaction of 
water with two single oxygen atoms of varying charge. It 
would be expected that as the charge on the oxygen atom 


increased (i.e., became more negative), then that atom would 
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became a better electron donor. Therefore, we have made cal- 
culations of the spectra of water in three such complexes, 
with M, the charge on the O atom varying from 0 to -1 to -2. 
The oxygen atoms with M = 0 are expected to be the weakest 
electron donors, and those oxygen atoms with M = -2 are 


expected to be the strongest electron donors. 


Tel Intensity Theory for Charged Species 


Before beginning a discussion of the predicted spectra 
for these complexes it is first necessary to present some 
additional theory pertaining to the prediction of infrared 
intensities of charged species. The main points will be 
reviewed briefly here; more details have been given else- 
where (55,56). 

The problem of calculating and interpreting APTs is 
slightly complicated for a charged species. First, the 
value of the dipole moment is determined by the choice of 
the origin, due to the charge on the molecule (92). We have 
adopted for these calculations the same unique coordinate sys- 
tem as for the normal coordinate calculations, namely the space 
fixed (not molecule fixed) principal cartesian coordinate 
system (with the origin fixed at the center of the mass by 
the Eckart conditions). The value of the dipole moment cal- 
culated for the charged species at the equilibrium geometry 
is arbitrary because of its dependence on this arbitrary 
choice of origin. However, we are interested in calculating 
dipole derivatives using the method of finite differences; 


hence the arbitrary choice of origin affects each of the 
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calculated dipole moments for each displacement in the same 
manner; therefore the calculated change in the dipole is 
meaningful. For example, the dipole derivative is approxi- 
mated using the method of finite differences with equation 5, 


rewritten here for convenience: 


oP dP P(X _+AX) - P(X _-AX) 
a a: 3. a." 16 a oO (5) 
aX aX 2AX 


where all the terms have been previously defined (Chapter 2). 
The values calculated for P  (Xo+4X) and P (X57 4X) are arbi- 
trary due to the choice of the origin. However, the displace- 
ments are constrained by the Eckart conditions to keep the 
center of mass fixed. The calculated values of PY (X+4X) and 
P (X74) are consistent with’ this constraint so that the 
dipole derivative aP,/9X is meaningful. 

The second complication in calculating APTs for charged 
species arises because there is a contribution to the dipole 
derivative as a result of the charge on the molecule. We have 
already seen that the dipole moment for a charged species is 
dependent on the choice of origin. Therefore, as a charged 
molecule is translated through space, its position relative to 
the origin changes and the calculated value of its dipole mo- 

ment changes. The subsequent change in the dipole moment as 
the molecule is translated through space represents the trans- 
lational correction to the dipole derivative. The translational 
contribution to the dipole derivative is automatically included 


in the value of the atomic polar tensors that we have calculated. 
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However, it is separate from the vibrational contribution to the 
dipole derivatives and makes no contribution to the intensities, 
so that the translational contribution should be explicitly 
separated if one is interested in analyzing dipole derivatives 
in order to explain vibrational intensities. It will be seen 
that the translational correction to the atomic polar tensors 
is, in most cases, not small. The APTs calculated for a molec- 
ular ion with this correction included cannot be compared with 
the APTs of the corresponding neutral molecule. Therefore, 
the translational contribution for the ion must be removed 
from the calculated APTs before comparing them with the APTs 
for the corresponding neutral molecule. 

The translational contribution to the APT of atom A in 


the ion is a diagonal tensor, given by 


Translational Correction = ——| q E (87) 


Here qe is the total charge on the molecular species in ques- 


tion (either positive or negative), M, is the mass of atom A, 


A 
and Mo is the mass of the molecule. Before these corrections 
are applied to the predicted polar tensors, the null condition 


(see equation 47) is 


rE Pha qe (88) 
N = 
where the sum runs over all the atoms in the molecule and the 
other terms are as previously defined. After subtracting the 
translational corrections to the APTs (given in equation 87), 


the null condition again returns to its original form (55) 
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because we have now corrected the dipole derivatives for the 


presence of the charge on the molecule 


E.Py Go = OE (89) 
yn a 
where Py Cc represents the "corrected APTs." 
La 


7.2 The oM+H,0:0" Complex 


The geometry for the oM-H,0-0% complex is given in Fig- 
ure 1.4. The principal cartesian coordinates for the assumed 
geometry for this molecule are given in Table 7.1. The same 
principal cartesian coordinates were used in the three differ- 
ent calculations. The charge on the oxygen atoms was varied 
from 0 to -1 to -2, but the geometry assumed for the complex 
was the same (Figure 1.4) for each calculation. The symmetry 


of the molecule is C.., with the 3N-6 = 9 vibrations corre- 


2v 
sponding to the irreducible representations 4A, + Ay + By - 
3B. and the three translations and three rotations belonging 


to A. + A. + 2B, + 2B In this orientation of atoms in the 


aL 2 1 2. 
complex, both of the free oxygen atoms are electron donating 
in nature, and both hydrogen atoms in the water molecule are 
electron accepting in nature. 

The predicted wavenumbers and intensities for the neutral, 
-2 and -4 charged species are given in Table 7.2 along with the 
corresponding monomer values. For all three cases the predicted 
wavenumbers and intensities of the vibrations of the complexed 


water are considerably different from the predicted values for 


the isolated monomer. The large discrepancies in the predicted 
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Table 7.1 Principal cartesian coordinates for the ener o-o! 


complex (M = 0, -1l, -2) (A°). Atom numbering as in Figfre 
: 


Xx Y Z 
0, 0 -.3010 2.4645 
Ho 0 -0549 - 8320 
O, 0 35952 0 
Hy 0 -0549 = 8320 
(e) 0 =, 3010 -2.4645 
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frequencies of the complexed water versus those of isolated 
water caused us some concern about whether or not our 
assignment was correct. However, we could easily investi- 
gate this question using the graphics program VIBRATE to 
display the normal modes. 

The general pattern for the frequencies is that the 
two stretching modes of the complexed water are shifted to 
considerably lower energy relative to the isolated monomer, 
while the bending mode is predicted to be unchanged from the 
monomer value. The corresponding predicted intensities for 
complexed water range from values approximately equal to the 
monomer values (for the neutral species), to values increas- 


ingly more intense (in the -2 and then the -4 charged species). 


7.3 Intensity Analysis for the o't+H,0-0" 
Complex 


The intensities calculated for these water complexes have 
been summarized previously in Table 7.2. The "corrected APTs" 
corresponding to these calculated intensities are given in 
Table 7.3. 

The "corrected APTs" reported in Table 7.3 are corrected 
for the translational contribution to the dipole derivative 
as described previously. However, these corrections have been 
applied in a manner unique to these complexes. We are calcu- 
lating the infrared spectra for neutral water complexed with 
charged oxygen atoms. Therefore, the charge present in the 


complex is almost entirely localized in these oxygen atoms. 


187 


€L40" Hz 0 heso"- — 10bT 0 yhoo" ROL" 
6461° E10" 0 chal: 100" 0 S610" 1s00°~ aw Fy 
7 
0 0 Bz" 0 0 LLLz" 9 0 O00" 
zee" 994 0 Tylt’ ots 0 6600° OLI0° 
gozt* 9 490"- 0 62790" 69 *- 0 €4lo’ —- Lolo*- . ate 
0 0 8CEC* 0 0 LLLe° 0 0 7000" - . 
BUeI’  E00°- 0 998° BSZh"- 6199" oles 0 9sat* 040"- 
Ob%0" 2696" 0 "S9S°- 9hE9° ISZe"- S465" 9 19IN"= 992° = ot 
0 0 oz" 0 0 cole’ 0 BELY" 0 0 "SLY" 
9784°- BhhO" 0 Zhe I- 6Izy" 0 9S80°I-  469z" 0 6EZS°- zhB0° 
B40" = WLIS*- sg 61Z%" 818" I- 0 "69%" ZZEI- 0 Z¥80" GBB *- = a 
“O 
0 0 0496"- 0 0 907° 1- 0 0 8Z0S *I- 0 0 6656" 
Ise’ Ltw'- 0 Bye’ ISTE’ 0 7696" —- BY6T*- 0 ozo" LLLO- 0 
£880°- zBrI" 0 SSLI°-  8961'T 0 9L8I°- 2068" 0 Ozol'- —-zeez" 0 “ "a 
0 0 0z8h" 0 0 Cole" 0 0 BELY" 0 0 SLY" 
z z z 
aIom oO . . 
ouoy 7-0.0°H, 70 1-0:0°H: 0-0°H-0 


“p’T *bta ut oases se we3sAs ojeUTpIOOD pue Hbutrequnu 


wozW °“(4X89 89S) peAOoWser useq Sey SLdWY eSey} OF UOTIOSeTTOO TeuoTAeTsues} oUuL “(9 67-OT 
xX Z09°T = 8T) (8) Z- ‘T- ‘0 = W ‘sexeTduos pO * OCH yO ey} JOFJ SLdV peyeATOD E°L STqeL 


188 


Hence the translational corrections are applied by assuming 
that the water molecule in the complex is neutral, and that 
the charge present in the complex results entirely from the 
charge on these oxygen atoms. The APTs for the charged oxygen 
atoms are "corrected" using this assumption while the APTs for 
the neutral water molecule are not "corrected." The resulting 
APTs are those listed in Table 7.3. 

Upon careful inspection of Table 7.3, the similarites and 
differences between the APTs of the atoms in the complexes and 
in the isolated water molecule can be seen. For the neutral 
complex, the predicted polar tensors for the atoms in the com- 
plex are very similar to the values predicted for the monomer. 
This is to be expected since the predicted intensities for 
this complex are very similar to those values predicted for 
the monomer (see Table 7.2). Based on these results one 
would expect that this interaction is relatively weak. 

The situation for H50 in the -2 and -4 charged complexes 
is entirely different. For the stretching mode of the charged 
complex we predict a rather drastic increase in intensity (al- 
most a factor of 60 and 30 in the symmetric and asymmetric modes, 
respectively). For the -4 charged complex these factors are 
even greater, being close to 94 and 50 respectively. 

The CCFO analysis of these APTs for these species offers 
some insight into the nature of these drastic predicted inten- 
sity increases. The predicted charge, charge flux, and over- 
lap tensors for the three complexes, and also for isolated 


water, are given in Tables 7.4, 7.5, and 7.6 respectively. 
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Table 7.4 Calculated charge tensors® for the a o-o! com- 
plexes, M = 0, -1l, -2, with and without the transl&tional 
corrections applied (see text) (e) (le = 1.602 x 10719 q). 
Atom numbering and coordinate system as given in Figure 1.4. 


0+H,0-0 o-*-H,0.-+ 0-*-H,0-077 Monomer 

With translational contribution 

O71 -.0062 -.8873 -1.7340 

Hy -4292 »4257 3340 - 4017 
03 -.8460 #150767 “1.1930 -.8034 
Hy, 4292 -4257 - 3310 -4017 
0. -.0062 =.8873 -1.7340 

With translational contribution removed 

Oo; -.0062 “1,27 - 2660 

Ho -4292 -4257 - 3310 - 4017 
03 - 8460 =1",,0'7.6 7 -1.1930 -.8034 
Hy - 4292 -4257 - 3310 -4017 
oO. -.0062 «L127 - 2660 


“The charge tensor is obtained by multiplying these Mulliken 
charges by the (3x3) identity matrix E; P (c) = Q,E. 
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The predicted APTs for the neutral complex will not be dis- 
cussed further because they are essentially unchanged from 
those of the corresponding monomer. Upon careful inspection 
of Tables 7.4, 7.5 and 7.6, it becomes obvious that the 
largest single contribution to the difference between the 
predicted APTs for the complexed H50 and the predicted APTs 
of the monomer is from the charge flux terms. However, there 
are also some significant differences predicted in the charge 
and overlap tensors. It is impossible to determine the exact 
contributions from each of these three tensors to the total 
predicted intensity increase of the complexed water by simply 
examining the APTs and the CCFO tensors alone. Therefore 
Table 7.7 becomes much more useful. 

Table 7.7 lists the calculated contributions to the total 
predicted intensity increase from each of the charge, charge 
flux, and overlap changes from monomer to complex. The second 
column gives the calculated intensities for the vibrations of 
the complexed water using the unchanged predicted monomer APTs 
and the normal coordinate transformation calculated for HO 
in the complex (the APTs for the charged(or free) oxygen 
atoms are included as the null tensor). This column 
therefore gives the relative contribution to the total in- 
tensities that occurs from the perturbed normal coordinates 
of the H,0 in the complex (as opposed to the monomer), or 
what has been termed the "mechanical effect" on the 
intensities. We see here that the mechanical effect on 


the predicted intensities is very small. For the case of 
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the -2 charged species this effect is almost nonexistent. 

For the -4 charged species this effect is a bit larger, but 
nevertheless still not large enough to be significant. 
Column three gives the calculated intensities obtained by 
including the contributions from the calculated APTs of the 
free oxygen atoms. This column indicates that the APTs for 
these atoms have very little effect on the total predicted 
intensities of the HO vibrations in the complex, and there- 
fore that these normal vibrations must not involve much 
motion of these oxygen atoms. This further indicates that 
the normal vibrations of the complexed water molecule are 
perturbed very little from the corresponding vibrations in 
isolated water. The next three successive columns show the 
relative contributions to the predicted intensity increase 
from the differences in charge flux, charge, and overlap be- 
tween the complex and the monomer in the P term of the APT 


¥y 


for H and will be described below. Finally, the last col- 


2" 
umn lists the total calculated intensities for the complex. 
The discrepancies between column six andsevenoccur because 
there are other minor contributions to the total intensity 
increase due to changes in other APT elements besides the 
Poy term of H5- 

Columns four, five and six of Table 7.7 were obtained 
in the following manner. The predicted APTs and the charge, 
charge flux, and overlap tensors are given for the complexes 


and for the monomer in the coordinate system shown in Figure 


Ted. The light hydrogen atoms are considerably more important 
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than the heavy oxygen atoms and therefore we focus our atten- 
tion on these two atoms. Using this "bond" coordinate system 
(Figure 1.4) it is easy to see the predicted changes in the 
APTs of complexed water relative to isolated water (see Table 
7.3). By examining Table 7.3 it is obvious that the Poy term 
of Ho is drastically different from the monomer value, while 
the remaining terms in the APT for this atom are similar to 
those in the monomer. Similarly, there are differences in 


the APT of Hy, but not isolated to one term as in H This 


2° 


result arises because if we place H., in the "bond" coordinate 


2 
system Hy is now oriented differently from Ho with respect to 
this coordinate system. However, by symmetry, the two atoms 
are identical. In other words, the APT for one hydrogen 

atom is exactly the same as the APT for the second atom as 
long as the proper rotations are performed. 

Columns four, five and six in Table 7.7 are obtained 
by adding the differences between the complex and the monomer 
in charge flux, charge, and overlap in the Poy term of Ho to 
the corresponding APT of the monomer. The APT for Hy is then 
obtained by properly rotating the adjusted APT for Ho. The 
APT for 0, is then obtained by insuring that the combined 
changes in both hydrogen atoms are cancelled by the changes 
in the oxygen atom (i.e., the null condition still holds). 
As each separate contribution (CCFO tensors) is added to 


the monomer APTs the intensities are recalculated using these 


APTs and the normal coordinate transformation for the complex.) 
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In this fashion each contribution (CCFO contributions) to 
the total predicted intensity is investigated separately. 

The predicted intensity for the vibrations of the com- 
plexed water molecule are considerably greater than those 
predicted for‘ isolated water (see Table 7.2). We have al- 
ready seen that the mechanical effect and the contributions 
from the "free" oxygen atoms are both negligible. The drastic 
predicted increase in intensity must therefore be a result of 
electronic effects reflected in the APTs. 

The -2 charged complex will be discussed first. The 
majority of the differences between the APTs of the complexed 
water molecule and isolated water are in the charge flux 
tensors. The Poy term of Hy (and the corresponding terms for 
Hy) are responsible for most of the predicted intensity in- 
crease in complexed water relative to monomer value (59% of 
the total increase in the asymmetric stretching mode). There 
is also a smaller but significant contribution from overlap 
(19%), and a relatively insignificant contribution from the 
change in charge (5%). 

For the -4 charged complex the results are similar. Both 
the mechanical effect and the effect of the "free" oxygen atoms 
are very small. The primary source of the intensity increase 
comes from the large increase in the charge flux as the complex 
is formed. This large increase is responsible for 73% of the 
total predicted intensity increase in the asymmetric stretching 
mode. The change in overlap offers a smaller contribution 
(16%), with the change in charge actually decreasing the pre- 


dicted intensities (-8%). 
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Qualitatively, the results for the -2 and -4 charged 
species are exactly the same. Namely that most of the in- 
tensity increase in the complex with respect to the monomer 
is due to a large increase in charge flux of the Poy term of 
the two hydrogen atoms. In addition, there is a smaller but 
nevertheless significant contribution from overlap in both 
cases. Finally, there is essentially no contribution (as the 
complex is formed) from the perturbation of the normal coordin- 
ates, from the "free" oxygen atoms, or from the change in charge. 

Quantitatively, there are some differences between the re- 
sults for the two species. The predicted intensity changes in 
the -4 charged species are definitely greater than the changes 
niedtonsd for the -2 charged species. When the numerical re-.- 
sults for the two species are compared, it becomes immediately 
apparent that the intensity increase in the complex with the 
stronger electron donor is considerably greater than in the 
complex with the weaker electron donor (by almost a factor of 
two). This suggests that a stronger electron donor would re- 
sult in larger predicted changes in intensities than a weaker 
electron donor. In addition, while the relative contributions 
from the change in overlap remain almost constant (16% vs. 19%), 
there is a rather substantial increase in the contribution from 
charge flux (59% vs. 73%) as the donor strength is varied. This 
would indicate that as the strength of the electron donor in- 
creases, the importance of charge flux in the predicted inten- 


sities increases also. 
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7.4 Frequency Analysis for the oM-H,0-0" 
Complexes 


The predicted absorption frequencies for the vibrations of 
the water molecule in these complexes have been summarized in 
Table 7.2, and compared to the predicted monomer values. In 
general, the frequencies for the stretching vibrations are con- 
siderably lower than those of the isolated molecule. In order 
to analyze these changes, the frequency parameters were trans- 
formed to force constants in internal coordinates. A potential 
function was assumed which does not include any interaction be- 
tween the intramolecular vibrations of the water molecule and 
the intermolecular vibrations of the complex. This seems rea- 
sonable because it appears from the intensity data that the 
normal vibrations of the neutral complexed water molecule are 
not much different from the corresponding vibrations in iso- 
lated water. Therefore, the coupling predicted between the 
intramolecular vibrations of complexed H50 and the intermol- 
ecular vibrations of the complex must be small. With this 
assumption a potential function which included only the intra- 
molecular internal coordinates of the complexed H50 molecule 
was used. The calculated force constants corresponding to 
these intramolecular internal coordinates are given in Table 
7.8. The frequencies calculated using this approximation are 
listed in column one of Table 7.9. A comparison to the fre- 
quencies from the "total" (i.e., including all possible inter- 
actions among vibrations; see Table 7.2) calculation shows 


that the agreement between the two calculations is very good. 


199 


Table 7.8 Calculated "effective" (see text) internal force 
constants (md/A°) for the internal coordinates of the water 
molecule in the OM-Ho0- complexes, M = 0, -l, -2. Internal 
coordinates are defined the same as for isolated water (see 
Table 3.6). Atom number and coordinate system as given in 
Figure 1.4, 


Force constants (md/A°) 


Monomer OH.0-0 07" +HL0+0 ~1 0-*-H,0-077 
Fi17F 55 9:14) 6.778 5.984 5.028 
F55 714 .874 -9437 .815 
Fi -.065 -.162 .064 .195 
F -300 -340 -090 -.083 
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Table 7.9 Predicted wavenumbers for the intramolecular 
vibrations in the 0”-H,0-0M complexes, M = 0, -1, -2, 
as several internal fofce constants are adjusted (see 
text). Force constants as defined in Table 7.8. 


Wavenumbers (em7 1) 


Mode Calculated Adjust Adjust Adjust Monomer 
Complex Fil and Fo> Fio F33 
nce ples a EN SNES Se 5 Se Td a OE oes NS PR 
0°H50-0 
an 3402 3957 3978 3978 3983 
1766 1778 1778 1601 1676 
3568 4130 4109 4109 4104 
O H50-0 
-1 -1 
O H50 oO 
3248 3454 3397 3397 3402 
1873 1872 1872 1802 1779 
3297... 3509 3568 3568 3570 
-2 -2 
O H,0 O 
1 3035 3475 3410 3411 3402 
1730 1734 1733 1795 1779 


2997 3500 3568 3568 3570 
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The discrepancies between the two sets of wavenumbers differ at 
most by 13 cm + with most values within 5 cm! of each other. 
We will begin our discussion of the predicted changes in 
absorption frequencies with the neutral complex. The predicted 
intensities for the neutral complex exhibited only slight 
changes from the monomer values, but the predicted stretching 
frequencies are considerably different from the corresponding 
monomer values. The stretching frequencies are predicted to 
be approximately 500 em lower in the complexed water than 
in the monomer, while the bending frequency is almost constant. 
Therefore, it is important to understand why such large changes 
are predicted in the stretching frequencies. First, it must be 
remembered that the assumed geometry for the water molecule in 
the complex (Ry = R, = 0.991 A®, and 6 = 114.0°) is consider- 
ably different from the calculated geometry of isolated water 
(Ry = Ry = 0.951 A°, and 6 = 111.2°). Therefore, before a 
direct comparison can be made of the calculated frequencies 
of the neutral complex and the isolated monomer, we must first 
determine what one might expect for the predicted spectral par- 
ameters (including both frequencies and intensities) of an 
isolated water molecule with the assumed geometry of the 
water molecule in the complex. 

‘If we recall the data presented in Table 3.3, the effect 
of the assumed geometry on the predicted spectral parameters 
can be seen. While the predicted intensities remain relatively 
insensitive to the assumed geometry, the predicted frequencies 


are very sensitive to the assumed geometry. The plot given in 
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Figure 7.1 of the O-H stretching force constant versus the O-H 
bond length shows that there is an inverse linear relationship 
between the two quantities (i.e., as the bond is lengthened, 
the force constants decreases). If this plot is extrapolated 
to 0.991 A° for the O-H bond length, the predicted stretching 
force constant is 6.634 md/A°. The frequencies calculated 
using this value in place of the value calculated for monomer 
water at the theoretical geometry (0.141 md/A°) are 3392 om, 


1665 cm 


and 350 cm for Vj" Yor v3 respectively. These fre- 
quencies are close to the values that one might expect for an 
isolated water molecule with the assumed geometry of the water 
molecule in the O-H20-O species. The agreement between the 
values predicted for the vibrations of H,0 in the neutral com- 
plex (see Table 7.2) with these values calculated for an iso- 
lated water molecule with this same assumed geometry indicates 
that the predicted spectrum for the neutral complex is what one 
might expect from an isolated water molecule with the same 
assumed geometry of the water molecule in the neutral complex. 
This suggests that in this complex the interaction between the 
water molecule and the neutral oxygen atom is weak, with little 
effect on the vibrational properties. 

The -2 charged complex exhibits a further decrease of the 
stretching frequencies relative to the neutral complex of about 
200 cm? for both stretching vibrations. Since the assumed 
geometry of the water molecule in this complex is exactly the 
same as in the neutral complex, this decrease can be attributed 


to the presence of an intermolecular interaction between the 
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Figure 7.1 Plot of O-H stretching force constant versus 
the assumed 0-H bond length for HO. 
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water molecule and the two "free" oxygen atoms. The calculated 
internal force constants for monomeric water, and for water com- 
plexed in the neutral, -2, and -4 charged species are listed in 
Table 7.8. From this data, the predicted decrease in the 
stretching frequencies can be attributed primarily to a de- 
crease in the O-H stretching force constant from 6.778 md/A° 
in the neutral complex, to 5.985 md/A° in the -2 charge complex. 
This change alone is responsible for 134% and 78% of the fre- 
quency shifts in the symmetric and asymmetric stretching modes 
respectively. The slight increase in intramolecular coupling 
between the two stretching coordinates (reflected in Fio pa 
Table 7.8) is responsible for the remaining shifts in stretch- 
ing frequencies, -37% and 21% respectively. The increase in 
the bending constant from 0.8744 md/A° in the neutral complex 
to 0.9437 md/A® in the -2 charged species is responsible for 
75% of the increase in the bending frequency (see Table 7.9). 
In the -4 charged complex, the results are qualitatively 
Similar. The predicted decrease in the stretching freqencies 
are greater than for the -2 charged complex, as would be expec- 
ted in the case of a stronger electron donor. These predicted 
shifts from the néutral complex are 367 cmt and 573 cmt for 
the symmetric and asymmetric stretching modes respectively. 
The predicted decrease in the stretching force constant from 
the neutral complex) from 6.778 md/A° to 5.028 md/A° is re- 
sponsible for 120% and 88% of these predicted shifts, 


respectively. The predicted increase in the intramolecular 


coupling between the two O-H stretching coordinates 
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reflected in Flo of Table 7.8) is responsible for the remaining 
-81% and 12%, respectively (see Table 7.9). 

Quantitatively, the frequency shifts predicted in the -2 
and -4 charged complex are considerably different, with the 
shifts of the symmetric and asymmetric stretching modes greater 
in the latter complex by 213 cm and 300 cmt, respectively 
(relative to the noninteracting neutral complex). Qualita- 
tively, the results of the frequency analyses for both com- 
plexes are similar. The asymmetric stretching mode is 
predicted to undergo a larger shift than the symmetric 


i. 


stretching mode in both cases (273 cm~ and 154 cm in the 


tN and’ 367" cn = iri Eievcd 


-2 charged complex, versus 573 cm 
charged complex). (It is worthwhile to note here that the 
predicted shift of the asymmetric stretching mode in the -4 
charged complex is such that this vibration is now of lower 
frequency than the symmetric stretching mode.) Most of these 
predicted frequency shifts are due to the large decrease pre- 
dicted in the O-H stretching force constant in both cases 
(over and above the decrease expected as a result only of 

the geometry change.) This data suggests that while the 
natures of the predicted frequency shifts are similar, a 
stronger electron donor is responsible for larger predicted 


shifts in the stretching frequencies than a weaker electron 


donor. 
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7.5 Model for the Infrared Spectra of 

Beta-Alumina 

Water, as an impurity in the structure of beta-alumina, 
is important in the mechanical strength and electrical conduc- 
tivity properties of the beta-alumina. This in turn is impor- 
tant in the use of beta-alumina as a solid electrolyte in 
batteries. Some of the kinetic and thermodynamic properties 
of the hydration reaction of water ina crystal of beta- 
alumina might best be studied by using the infrared spectrum 
as a probe (44,45). Hence, the crystal structure and infrared 
spectrum of water isolated in beta-alumina has been investi- 
gated recently by Bates and coworkers (44). Figure 7.2 gives 
the orientation of the water molecules in the crystal. The 
orientation is such that the nearest neighbors of the water 
molecule are two oxygen atoms belonging to the aluminate ions. 
The oxygen atoms in the aluminate ion are most assuredly nega- 
tively charged, but the exact numerical charge on these atoms 
is unknown. We have previously calculated in the spectrum of 
one water molecule interacting with two neutral, -l, and -2 
charged oxygen atoms. [In our calculations, the assumed geom- 
etry for the water molecule (see Figure 1.4) is the same 
geometry determined by Bates et al. (44) to be the experimental 
geometry for water isolated in a crystal of beta-alumina (see 
Figure 7.2) In addition, the assumed orientation of the 
charged oxygen atoms in our calculations is the same as deter- 
mined by the crystallographic studies on the crystal by Bates 
et al. (44). The structure for our simulation has been pre- 


viously given in Figure 1.4, where the charge on the oxygen 


All 


Figure 7.2 Orientation of the H.o molecule in a hydrated 


crystal of beta-alumina. (Figuré reprinted from reference 
44.) 
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atoms (O, and O05) varies from 0 to -l to -2. We have seen 
previously that the spectra calculated using this model are 
quite different from the spectrum calculated for the isolated 
water monomer. As a first approximation, we might also expect 
the experimental spectrum for water isolated in beta-alumina 
to exhibit some of the predicted changes that were seen in 

our model; but it must be remembered that our calculation is 
only a model for what might be expected in the crystal and 
that the agreement between our calculated results and the 
experimental results cannot be expected to be exact. 

The experimental wavenumbers and intensities for water 
isolated in beta-alumina have been given in Table 7.2, along 
with our calculated results for comparison. The best agree- 
ment with our calculated results and the experimental results 
eis when M = ~ © or when the charge on the oxygen atoms in 
the aluminate ions is -l. The calculated frequencies for the 
symmetric and asymmetric stretching vibrations are predicted 
83 cmt and 153 cm=+ higher than the experimental values, while 
the bending frequency is predicted 359 cmt too high. The pre- 
dicted intensities for these three modes are predicted 203 
km/mole too low for the symmetric stretch, and 61 km/mole and 
47 km/mole too high for the asymmetric stretching and bending 
vibrations, respectively. While the agreement between these 
calculated results and the experimental results is only fair, 
the differences are greater as the charge on the oxygen atoms 
is increased to -2. However, as the charge on the oxygen 


atoms is increased to -2, most of the predicted differences 
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between the calculated and experimental results change sign 
relative to the differences predicted when M = -l. This sug- 
gests that in the crystal we might expect the charge on these 
oxygen atoms to be somewhere between -1l and -2, but closer to 
-1. In our model calculation, the vibrational bands predicted 
for the complexed water molecule were shifted to lower energy 
and increased in intensity relative to the isolated monomer. 
It was also shown that these predicted spectral shifts were 
due to an electronic interaction bewteen the water molecule 
and the charged oxygen atoms. The agreement in both frequen- 
cies and intensities with the experimental spectrum suggests 
several important points. First, the assignment of the infra- 
red spectrum made by Bates et al. (44) appears to be correct. 
Secondly, the large spectral shifts observed experimentally 
for the vibrations corresponding to the complexed water mole- 
cule in the crystal indicate the presence of strong hydrogen 
bonding in the crystal. Finally, in the crystal, the nearest 
neighbor oxygen atoms of the aluminate ions appears to be the 
most important atoms in trying to understand this intermolecular 
interaction. As a first approximation, these atoms alone are 
responsible for most of the observed spectral changes for the 
water molecule in the crystal compared to isolated gaseous 


water. 


CHAPTER 8 


IONIZED WATER, H,0" 


Ionized water, H,0", represents a different challenge for 
these types of calculations in that it is not an intermolecular 
complex, but rather an ionic species which may be 
stable in this form. The calculation of its infrared spectrum 
is not directly related to the intermolecular complexes pre- 
viously discussed, but it does represent a spectrum that one 
might expect for water in a related chemical environment 
(complete electron donation), and therefore is of interest 
here. 

The geometry, and therefore the principal cartesian co- 
ordinates, were taken to be the same as the equilibrium theo- 
retical geometry for neutral water using the 4-31G basis set. 
The two O-H bond lengths are assumed to be equal and with a 
value of 0.951 A°, while the H-O-H angle is assumed to be 
111.2°. The predicted wavenumbers and absolute intensities 
for H,0" at this geometry are listed in Table 8.1. In gen- 
eral, the frequencies are predicted to be more or less the 
same as the frequencies for the neutral molecule, but the 
intensities are predicted to be considerably larger in the 


ion. 
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Table 8.1 Predicted wavenumbers and intensities for H,0" 
compared to those predicted for H,0. 


Mode Description Calculated H,0" Calculated H50 
viem™*) A(km/mole) v(cm™4) a(km/mole) 


ees 


IL Symmetric 4012 350 3983 3.9 
stretch 

2 Bending 1545 303 1676 125 

3 Asymmetric 4033 642 4104 54 


stretch 
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Since this molecule is charged, the value of the dipole 
moments depends on the choice of the origin. Since we are 
interested in calculating dipole derivatives using the method 
of finite differences, we may fix the origin at the center of 
mass for consistent estimation of the dipole derivative (see 
discussion in Chapter 7). As discussed also in Chapter 7, 
the charge on the molecule requires a "translational correc- 
tion" (56) to the dipole derivative (APT) before we attempt 
to compare the APTs of an ion with the APTs of the correspond- 
ing neutral molecule. The APTs for ionized water both with 
and without the translational contribution are reported in 
Tables 8.2 and 8.3, respectively, together with the corre- 
sponding APTs for the neutral molecule in the coordinate 
system as defined in Table 8.2. By comparing Tables 8.2 
and 8.3 it becomes immediately apparent that the APTs re- 
ported with and those reported without the translational con- 
tribution are considerably different. The "corrected" APTs 
given in Table 8.3 should be compared with the APTs from 
neutral water. There are considerable differences between 
the corrected APTs of the ion and the APTs for the neutral 
molecule. These differences were analyzed in terms of the 
CCFO analysis. The charge, charge flux and Overlap tensors 
for the ion and neutral molecule are listed in Tables 8.4, 

8.5 and 8.6, respectively. 

The charge contribution to the APTs of the ion are listed 

in Table 8.4 with and without the translational contribution. 


The values that include the translational contributions are 
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Table 8.2 Calculated APTS (e) for ionized water, H»07, com- 
pared with the calculated APTs for monomeric neutral water 
(le = 1.602 c 10719 c). Atom numbering and coordinate system 
as defined below. (Includes translational contribution.) 


y- a 
ee 
Be as 
+ 
H,0 H,0 
479 0 0 482 0 0 
pl. 0 484 044 0 5370 -- =n 
0 .115°> .730 0 003 .130 
043 0 0 S964" 0 0 
02 
p*s 0 =.235°>™-.055 0; =,5E2. <.04s 
0 <3088. =,19 OG. a beelh acana 
479 0 0 .482 0 0 
p> 0 i Ly a el 0 148 089 


0 -.060 -462 0 -042 +352 
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Table 8.3 Corrected APTs_for H o* compared with the APTs for 
H,O (e) (le = 1.602 x 1071 C). “Atom numbering and coordinate 
system as defined in Table 8.2. (Translational contributions 
have been removed, see text, Chapter 7.) 


H0 H0 
423 0 0 482 0 0 
pil = 
x7C 0 428 044 0 370 =. 044 
0 115 674 0 003 .130 
=, 848 0 0 -.964 0 0 
0, 
eo ies 0 agi 22. ~,056 0 =,817 ~=coas 
0 /,955-- *=,1080 0 w/045 . =. 48 
f 423 0 0 482 0 0 
hea = 0 696 O11 0 148 089 
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Table 8.4 Calculated charge tensors® (e) (le = 1.602 x 
10-19 ¢c) for H,0* compared to H20. Atom numbering and coor- 
dinate system as defined in Table 8.2. 


With translational contribution 


Hy 595 - 402 
05 ~.196 -.804 
H3 595 - 402 


With translational contribution removed 


Hy F530 -402 
05 -1.078 -.804 
H3 3939 - 402 


“The charge tensors are given by Q E, where Qn is the Mulliken 
charge for atom A as listed here. 
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Table 8.5 Calculated charge flux tensors for H.O* compared 
to H20 (e) (le= 1.602 x 10719). atom numbering and coordin- 
ate system as defined in Table 8.2. 


+ 
H,0 H,0 
0 0 0 0 0 0 
Hy 
P(CF) = 0 .060 200 0 -.043 .039 
0 .186 4.12 0 106 052 
0 0 0 0 0 0 
0, 
P(CF) = 0 -.556 -.216 0 -.045 -.093 
0 S716. a 38 0 -.093 027 
0 0 0 0 0 0 
Hy 
P(CF) = 0 .496 .016 0 089 054 
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Na 8.6 Calculated Tort tensors for H o* compared to 
O--(e).2(Le:=- 10602. x<10™ Atom numbering and coordin- 
ata soe tan as defined in vais 8.62). 


+ 
H,0 HO 
x -.116 0 0 .080 0 0 
1 
P(O) = 0 -.171 -.156 0 O12 -.083 
0 -.071 <-.277 0 -.103 -.324 
.233 0 0 -.160 0 0 
02 
P(O) = 0 511 .161 0 332 048 
0) 161 . 386 0 048 295 
-.116 0 0 .080 0 0 
H3 
P(O) = 0 -.339 -.005 0 -.343 035 
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simply the Mulliken equilibrium charges on the atoms in the 
ion. These Charges represent the actual distribution of 
electron density on each of the atoms in the ion, and repre- 
sent the total contribution of these static charges to the 
APTs. The charges reported with the translational contribu- 
tion removed represent the vibrational-rotational contribution 
to the APTs from only the static Mulliken charges (56). 

The charges with the translational corrections included 
indicate that approximately 0.60 e is removed from the oxygen 
atom and 0.20 e from each of the two hydrogen atoms as the 
ion is formed. This Suggests that the ion forms primarily by 
loss of a lone pair electron on the oxygen. The contribution 
from the uncorrected charge tensor to the APT for the oxygen 
atom in the ion is therefore less (closer to zero) than the 
charge contribution in the neutral molecule (-.190 e vs. 
-.804 e). The charges reported with the translational con- 
tribution removed (or the corrected charge tensors) actually 
show that the charge contribution to the Apt of the oxygen 
atom is greater (more negative) in the ion than in the neu- 
tral molecule (-1.078 e vs. -.804 e). This means that the 
charge contribution to the corrected APT (or vibrational- 
rotational APT) of the oxygen atom is significantly greater 
in the ion than in the neutral molecule. This result is 
important because it is the first time that we have seen a 
situation where the predicted changes in the APTs contain 
such a large contribution from anything other than the 


charge flux tensor. 
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The charge flux and overlap tensors for the ion and the 
neutral molecule in the coordinate system defined in Table 8.2 
are listed in Tables 8.5 and 8.6 respectively. The two hydro- 
gen atoms are identical by symmetry, but the tensors reported 
here differ because in the bond coordinate system used here 
the orientation with respect to the coordinate system is dif- 
ferent for these two atoms. We chose to use the bond coordin- 
ate system because in this coordinate system the differences 
in charge flux between the ion and the neutral molecule are 
isolated to the Poy term for H (and correspondingly to the 
rotated tensor for Hy). There are almost no differences in 
the overlap tensors for the ion and the neutral molecule. 
Therefore, it seems that the biggest contributions to the 
changes in the total APTs between the ion and the neutral 
molecule are due to changes in charge and charge flux. 

The contribution to the intensity change from the change 
in normal coordinates as the ion is formed (mechanical effect) 
and the contributions from each of the three constituent ten- 
sors of the APT for ionized water are given in Table 8.7. 

The differences between the normal coordinates of the ion 

and the neutral molecule are negligible, and do not account 
for any significant differences in the intensities of the 

ion. The difference in charge flux between the ion and the 
neutral molecule is responsible for the majority of the in- 
crease in the asymmetric stretching mode (79%) and for 38% 

of the predicted intensity increase for the symmetric stretch- 


ing mode. The overlap contribution is negligible, while the 
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Table 8.7 Calculated contributions to the predicted inten- 


sities of H,0* (km/mole) . 


Mode Using Monomer Add Add Add Calculated 
APTs H H H + 
P ACF P xO P xc H,O 
yyicF) P40) Bc) 9 
1 ae 136 140 Zs 350 
2 126 155 156 167 303 


3 54 516 525 788 642 
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charge contribution actually overestimates (by 146 km/mole) 

the predicted intensity increase for the asymmetric stretch- 
ing mode. The predicted intensities in column five of Table 
8.7 do not agree with the values predicted for the ion (col- 
umn six) due to some contributions to the APTs other than 

the Pyy term of H3- However, the large contributions from 
both charge flux and charge to the predicted intensity increase 
in the ion relative to the neutral molecule are important. This 
will be discussed in more detail in Chapter 9. 

The force constants for ionized water compared to those 
for the neutral molecule are given in Table 8.8. The major 
differences between the two are the interaction constant be- 
tween the two stretching coordinates (Fi 5) and the interaction 
constant between the stretching coordinate and the bending 


coordinate (F These differences are relatively minor, 


de 
and are responsible for the minor differences predicted in 
frequencies between the two species. The increased inter- 
action constant would serve to increase the frequency of 

the symmetric stretching vibration, and decrease the fre- 

quency of the asymmetric stretching mode. This is exactly 
the result shown in Table 8.1. The slight decrease in the 
bending constant is responsible for the slight decrease in 
the bending frequency seen in Table 8.1. The fact that the 
stretching constants are not affected very much would indi- 
cate that the electron removed when forming the ion is of 


nonbonding character. But again, this discussion will be 


given in more detail in Chapter 9. 
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Table 8.8 Calculated internal force constants for HO" com- 

Sey to facies Internal ae are defined as Rj(09-H)) = 
(O5 H3 =d (Hy -O 5-H - Atom numbering as in Tabte 

(Ret gnezas is donk with? : = 1 A°, see text, Chapter 3). 


Force Constant (md/A°) 


H,0" H50 
Fil = F5> 9.03 9.14 
F33 640 -714 
Fid -144 -.065 
F 469 - 300 


CHAPTER 9 
GENERAL COMMENTS AND FINAL SUMMARY 

Ab initio self consistent field calculations using the 
4-31G basis set have been used to predict and interpret the 
infrared spectra of water molecules interacting in different 
chemical environments. These environments have varied with 
respect to orientation of the electron donor/acceptor pair, 
with respect to the assumed intermolecular distance, and with 
respect to the strength of the electron donor. The predicted 
spectra are strongly dependent upon all these parameters, and 
by analyzing these spectra we have been able to make some com- 
ments concerning the nature of intermolecular interactions. 

The predicted spectrum for the linear water dimer (sum- 
marized in Table 4.2) compares favorably with experimental 
spectra (19-22). The symmetric stretching vibration of the 
electron acceptor molecule is predicted to shift to lower 
energy and increase in intensity, as is characteristic of 
hydrogen bonded systems (43) and as is found in the experi- 
mental spectra. The frequency and intensity predicted for 
this vibration are also seen to be rather sensitive to the 
assumed intermolecular distance. The origin of the predicted 
frequency shift is the predicted change in the strength of 
the O-H bond as indicated by the predicted change in O.-Hy 


bond length in the electron acceptor molecule of the dimer 
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relative to the O-H bond length in the isolated monomer. The 
origin of the intensity enhancement is predicted to be due 
entirely to the increase in charge flux for the bonded hydro- 
gen atom of the electron acceptor molecule. 

The cyclic and bifurcated forms of the water dimer are 
predicted to be less stable energetically (40), and the spectra 
for these two forms are not predicted to exhibit the same dras- 
tic spectral changes as for the linear dimer. In both the 
cyclic and bifurcated configurations previous calculations 
(36) have indicated that the geometries of the individual mon- 
omers in these complexes are not perturbed from the geometry 
of the corresponding isolated monomer. Therefore, we have 
fixed the geometries for the individual water molecules in 
these configurations and calculated the spectra with varying 
intermolecular distances. In the cyclic configuration, the 
frequencies of the stretching vibrations shifted only slightly 
from the corresponding monomer vibrations, and these shifts 
were due to the presence of increased intramolecular coupling 
in the complex relative to the monomer. The bending vibra- 
tions were predicted to shift slightly more than the stretch- 
ing vibrations, and these predicted shifts were seen to be 
due to intermolecular coupling between these vibrations and 
the low frequency vibrations. The intensities for the cyclic 
configuration were predicted to be slightly higher than those 
of the corresponding vibrations in the monomer, but the pre- 
dicted increases were not nearly as drastic as those for the 


linear dimer. These predicted changes in intensity were seen 
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to be a result of two effects in almost equal contributions; 
first, the mechanical effect due to the symmetry of the cyclic 
configuration (C;) and the perturbation of the form of the 
normal vibration in the complex relative to the same mode in 
the isolated monomer, and second, the electronic effect (or 

the increase in charge flux for the bonded hydrogen atom). 

The predicted changes in the spectral properties (both fre- 
quencies and intensities) of the cyclic configuration rela- 
tive to the isolated monomer were seen to be almost independent 
of the assumed intermolecular distance. Finally, in order to 
determine if perhaps these spectral changes were induced by 

the symmetry of the cyclic complex, the spectrum for the cyclic 
configuration without a center of inversion was also predicted. 
This spectrum was predicted to be quite similar to the predicted 
spectrum for the cyclic configuration with the center of inver- 
sion (if the different selection rules for these two symmetry 
groups are properly accounted for). 

The predicted spectrum for the bifurcated form of the 
water-water complex was found to be very similar to the spec- 
trum one might expect from two noninteracting water molecules. 
The predicted frequencies in the bifurcated complex agree very 
well with the predicted frequencies of the corresponding vibra- 
tions in isolated monomeric water. The intensities were pre- 
dicted to change in the complex relative to the monomer, but 
the predicted changes were due to the mechanical effect cou- 
pling the vibrations of the two molecules, and not. the result 


of the electronic effect or change in charge flux. The 
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predicted spectrum for the bifurcated form of the complex was 
also almost entirely independent of the assumed intermolecular 
distance. 

The predicted spectra for the water molecule in the 
o'-H,0-0™ complex, with M varying from 0 to -l to -2, indi- 
cated that the spectral properties of a hydrogen bonded 
electron acceptor molecule are strongly dependent upon the 
strength of the electron donor: the stronger the electron 
donor the larger the predicted changes in the spectrum of the 
complex (at a particular fixed geometry) relative to the iso- 
lated monomer. This was found to be true for both the pre- 
dicted shifts in frequencies and for the predicted increases 
in intensities. The frequency shift occured because the O-H 
stretching force constant was predicted to decrease, while 
the predicted intensity increases were due primarily to pre- 
dicted increases in charge flux. The prediction that the 
intensities of intermolecular complexes are determined by 
charge flux seems to be ubiquitous in these calculations. 

The importance of the charge flux tensor will be briefly 


discussed below. 


9.1 Comments on the Charge Flux, Charge, 
and Overlap Tensors 


Previous attempts to parameterize infrared intensities 
have used various parametric methods such as APTs (7,8) and 
electrooptical parameters (EOPs) (5,6,9), which include contri- 


butions from static charges, charge flux, and overlap. 
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However, the intensities for intermolecular complexes appear 
to be determined by the distribution of electrical charge at 
equilibrium and the changes in these charges during atomic 
motions (i.e., vibrations). Therefore, it appears to be of 
some use to attempt to parameterize charge flows in intermo- 
lecular complexes. The following discussion of charge flows 
will focus its attention on the bonded hydrogen atom because 
at this point it seems clear that this is the most important 
atom to consider. 

Gussoni, Jona, and Zerbi (93) have attempted recently to 
understand the infrared intensities for vibrations of C-H bonds 
by measuring what they have termed "charge undeformability." 
They have defined this term with respect to the invariants of 
the APTs. This quantity is defined (93) as the ratio of the 
mean dipole moment, D, to the anisotropy Bae This parameter 
measures the resistance to charge flow; it is infinite for 
completely undeformable charges, and small for highly deform- 
able charges (93). 

In the limiting case of a heteronuclear diatomic molecule 
AB and the EECF model of Decius (68), this ratio takes the 


following form (93): 


—A a 
iE | = Be + L (90) 
Ba a 3 

qa ~° 

ol ap AB 


where the dipole moment of a heteronuclear diatomic is given 


by 
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Pp = qn (91) 


so that the charge flux and overlap terms of the CCFO model 
: : fe) fe) 

are incorporated into the term (8q,/3r,_) Tap: 

For the same heteronuclear diatomic AB, the CCFO model 


gives this same ratio in the following form 


—A 
Pine ieee ae = 
3 ery + 3 (92) 
A 
where the APT for atom A is defined below along with the 
definitions of the terms a and b 
o A 
a 0 0 Qn + P(O) 0 0 
A 5 PA A 
Py = 0 b 0 = 0 G,tP (0) +P (CF) 0 
6 A 
0 0 a 0 0 q,tP (0) 
(93) 


A 
where dp is the equilibrium charge on atom A, P(0O) is the 


overlap contribution along the two degenerate axes perpen- 
dicular to the AB bond, P (0) is the overlap contribution 
along the molecular axis, and P (CF) is the charge flux con- 
tribution along the molecular axis. Substituting into 


equation 92 gives the ratio as 
A 


° 
—A q, + P(O) 
ee ret. ree ae (94) 
' 
By A A A 3 


P(CE) +2 (0) =P (0) 


By representing the ratio |p*78, | in terms of the CCFO model 
we see that this ratio is not actually a pure measure of the 


charge undeformability because it contains contributions 
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from overlap also. It would be more accurate to represent 
charge undeformability as 


(eo) 
qa 
A 
P (CF) 


Charge undeformability = (95) 


because this ratio represents the actual resistance to the 
change in equilibrium charge on atom A without any contribu- 
tion from overlap. 

We attempted to correlate our predicted increased charge 
flux in the complexes discussed here, and also with other com- 
plexes (91), with predicted decreased values of charge unde- 
formabilities using [p*7a, | as defined by Gussoni, Jona and 
Zerbi (93). These data are presented in Table 9.1; unfortu- 
nately they do not seem to correlate well. (For example, 
the predicted increased charge flux either at shorter inter- 
molecular distances or with stronger electron donors is 
expected to be reflected in the smaller values of [P7853 |.) 
We believe that these discrepancies arise in part due to the 
presence of the overlap contribution in the ratio |P"/e,|- 

As an alternative definition for charge undeformability we 
propose to use the effective invariants for the CCFO tensors 
to give a parameter analogous to the ratio expressed by 
equation 95. This ratio, which is defined in the general 


case for any polyatomic molecule is given by 


(96) 


Charge undeformability = 
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Table 9.1 Correlation between different charge undeforma- 
bility parameters and increased charge flux for different 
intermolecular complexes. (Smaller ratio should follow 
larger charge flux.) 


H 
H == 5 X(c | 
Complex x (2) — ce 
: H b 
(e) b X (CF) 
Linear Dimer 
R = 2.532 aA° .3820 1.4232 1.280 
R = 2.682 A° .3158 1.593 1.509 
R = 2.832 A’ . 2668 2.085 1.739 
R = 2.982 A Bye ip, 2.55 1.972 
R = 3.132 A° .1998 2.657 2.203 
a 
H,0---CO 
R = 2.832 A .1359 1.885 2.927 
Riess. 132 A a BLS. 1.543 3.390 
R = 3.532 A .1021 1.399 3.923 
R = 3.732 A .0969 1534] 4.133 
Cyclic Dimer 
R= 2.74 A .1510 1.346 2.855 
Bifurcated Dimer 
R = 2.89 A° .0695 1.305 5.890 
-1 -1 
0 ~*H,0-0 .4205 1.065 1.108 
-2 ~2 
0 “+H,0-0 26545 .7279 .6283 


@Rreference 91. 


Zak 


where Xe) and Yer) are the previously defined effective in- 
variants (see equations 71 and 78) for the charge tensor and 
the charge flux tensor respectively. From analysis of these 
data presented in Table 9.1, it appears that this ratio does 
correlate better with our predicted increased charge flux 
(i.e., as the charge flux increases, the charge undeforma- 
bility decreases). Therefore, we suggest that the ratio 
Cae eee is a good parameter to measure charge flows in 
these complexes. 

While the charge flux tensor appears to be the most 
important contribution to the predicted intensities of inter- 
molecular complexes, we have seen that in ions the contribu- 
tions from the charge and overlap tensors cannot be completely 
overlooked (when comparing the intensities of these species 
with the parent or corresponding neutral molecule). Recent 
calculations on other ions (18) and radicals (94) also point 
out the importance of these two component tensors in infrared 
intensities. These results will be discussed below. 

The overlap tensor represents the change in the dipole 
moment due to the rehybridization of charge, without a change 
in charge on a particular atom, as an atom in the molecule 
vibrates. The charge tensors are represented by the static 
equilibrium Mulliken charges (70) on the atoms in a molecule, 
complex or ion. In general, the charge and overlap contribu- 
tions to the predicted increases in intensities in an inter- 
molecular complex relative to the isolated monomer are 


negligible. Therefore, these calculations offer little 
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insight into the nature of both of these tensors. However, 
calculations reported recently on molecular ions (18) and 
radicals (94), and the comparison to the corresponding 
neutral species have been helpful in characterizing both 
the charge tensor and the overlap tensor, and their impor- 
tance in vibrational intensities. 

The calculations were all performed in the same manner 
as previously described (see Chapter 2), and the results are 
presented elsewhere (18,94). These results will be briefly 


summarized here. In general, the intensities of the asym- 


+ + 
2 a 


predicted to be less than the intensities in the correspond- 


metric stretching vibrations of CO (18) and CS (18) are 


ing neutral molecules, while those of H o* (see Chapter 8) 


2 


are predicted to be greater than those of H,O. The differ- 


2 
ences in charge flux between the ion and the neutral mole- 
cule are primarily responsible for these predicted differences 
in intensity; however, there are also important contributions 
from both changes in charge and overlap. It is these latter 
changes which lead to some important conclusions concerning 
the nature of these tensors in vibrational intensities. 

The asymmetric stretching vibrations for these ions have 
been analyzed in terms of the CCFO theory, and the various 
contributions to the differences in intensities have been 
analyzed using the program PXPQ (see Appendix A). For the 
ions ad and H,0” there is a considerable contribution to 


the predicted changes in intensity (relative to the corre- 


sponding neutral molecule) from the charge tensor (21% and 
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45% respectively). For od the contribution from the charge 
tensor is minimal, but there is instead a significant contri- 
bution (45%) from the overlap tensor. 

Presumably, the differences in intensity between an ion 
and a neutral molecule (for these positive ions) are due to the 
absence of an electron in the ion. Therefore, the differences 
in the contributions to the changes in intensity would be found 
logically in the nature of the electron removed. The contri- 
bution from charge flux appears to be ubiquitous and not re- 
lated to the nature of the electron removed. In the cases of 
go,” and sy" it was found that the differences in the charge 
flux could be explained on the basis of different resonance 
structures in the ion relative to the neutral molecule (18). 
However, the contributions from charge and overlap appear 
only in certain situations, and therefore comments on the 
nature of these tensors can be offered. 

es. and H,0" are both ions which are formed by the re- 
moval of a lone pair electron. In os,” this electron is 
removed from the S atoms, and in H,0" from the O atom. Con- 
sequently, there are subsequent changes in the static Mulli- 
ken charges on these atoms in the ions. Therefore, there is 
a significant contribution to the changes in intensity from 
the charge tensor. oo. is formed by the removal of a bond- 
ing electron from co, (18). Therefore, the overlap contri- 
bution to the APTs in the ion is different from that in the 


neutral molecule. Based on these results we believe that 


the importance of the charge tensor lies in the lone pair 
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effects, and the importance of the overlap tensor more in 


the bonding electrons. 


9.2 Other Applications for These 
Calculations 


While we have used these calculated spectra to investi- 
gate some of the more important aspects of intermolecular 
interactions, these types of predictions may also be useful 
in other areas. For instance, by comparing the predicted 
spectral pattern for the linear, cyclic and bifurcated water 
dimers with experimental spectra (19-21), it is obvious that 
the spectral pattern for the linear configuration is in best 
agreement with the experimental spectra. While the agreement 
between experiment and theory is not perfect, the lack of 
agreement between experiment and theory for the cyclic and 
bifurcated forms seems to indicate that the interpretation 
of the experimental spectrum in terms of a linear dimer is 
indeed correct. While these calculations do not exclude the 
presence of the cyclic (or bifurcated) dimers, they strongly 
Support the idea that the linear configuration dominate the 
spectrum. 

As another example, one can determine the relative con- 
centrations of the two species in a given sample (for example 
in a matrix) by using the predicted values of the integrated 
intensities for the monomer and for the dimer. This deter- 
mination is made using equation 97 below where ay is the 
integrated absorbance of one of the monomer absorption bands 


(obtained from the spectrum), A is the calculated integrated 
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molar absorption coefficient for that band, 2 is the path- 


length, and C is the concentration 


ay = AuCcu? (97a) 
or 

“> _ *u@p (97b) 

Cc. Bes 

M DM 


Using these relative concentrations from equation 97b, the 


equilibrium constant for the reaction 2 H50 = (H50) » is 


[dimer ] “1 Cp (98) 


eq [monomer ]? (e.,) 


under the experimental conditions of that particular sample 
can be roughly determined (95). 

The calculation of the vibrational spectrum for the 
o-H,0-0" complex as a model for the infrared spectrum of 
water isolated in a crystal of beta-alumina indicated that 
it is possible to model complicated chemical systems using 
rather simple theoretical models, at least to a first 
approximation. This is important because it offers hope 
that vibrational spectroscopic problems that have in the 
past not been investigated by theoretical calculations be- 
cause they were too complicated or too large for economically 
reasonable and reliable calculations may now confidently be 
attacked using simpler models. 

Finally, we should like to offer some additional comments 


on the calculations of infrared spectra which have not been 
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discussed in this thesis but which seem promising for future 
applications of this work. Vibrational spectroscopy is often 
used to characterize species which cannot be obtained ina 
bottle and isolated in pure form (96). In these cases, 

where standard spectra are often absent, the prediction of 
these spectra would appear to be of great help in trying to 
identify these species (for instance in flames, in planetary 
atmospheres, or on catalytic surfaces). The predictions of 
the integrated intensities could also be used to estimate 
roughly the relative concentrations of these species in a 
particular sample, or if the pathlength is known, the abso- 
lute concentration. For example, this is important in the 
spectroscopy of astrophysical sources (such as evolved stars), 
where the absolute concentrations of species believed to be 
important in these sources are important in the development 
of thermo- and chemical models for these sources (97). Pre- 
liminary work has recently been done in this area (18,94), but 
more detailed investigations are necessary to full understand 
the capabilities for predictions in this field. 

Understanding of the infrared spectra of large (> 10 
atoms) biologically important molecules may be an important 
tool in reaching an understanding of the myeries of living 
systems. However, the size of these molecules and the even 
larger ones in living systems often means the spectra are 
very complicated, with many vibrational infrared absportion 
bands close in energy to each other. This makes the inter- 


pretation of these spectra extremely difficult, and sometimes 
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impossible, even for isolated 12 atom molecules (17). Ab 
initio calculations of vibrational force fields (and hence 
vibrational frequencies and the normal coordinate transfor- 
Mation matrix) can be used in conjunction with ab initio 
APTs to predict both infrared absorption frequencies and 
integrated intensities. These predicted spectra must match 
the experimental spectra (including both frequencies and 
intensities) if the assignment is correct. Once a reasonable 
theoretical prediction has been obtained, the corresponding 
normal modes can then easily be assigned from the calculated 
normal coordinate transformation and the visual vibrating 
Program VIBRATE (see Appendix A). This has just recently 
been done for the 12 atom molecule uracil (17), and the re- 
sults are very encouraging. Work in this area should be 


continued, as the future Prospects seem extremely bright. 


PRINCIPA 


APPENDIX A 
SUMMARY OF PROGRAMS 


Disk Input - none 


Operator Input - 


Terminal Output- 


Disk output - 
Description - 


ORTHO 


masses and cartesian coordinates in any 
arbitrary coordinate system 

principal cartesian coordinates, princi- 
pal moments of inertia, transformation 
matrix from principal cartesian coordin- 
ates to the original coordinates. 
principal cartesian coordinates 
calculates the moment of inertia tensor 
in the original coordinate system and 
then diagonalizes that tensor to find 
the transformation to principal carte- 
Sian axes. Also calculates the center 
of mass, and then uses this as the 
origin and the transformation to prin- 
cipal cartesian axes to find the 
principal cartesian coordinates. 


Disk Input - principal cartesian coordinates 


Operator Input - 


Terminal Output- 
Disk Output - 
Description - 


masses, principal moments of inertia, 
basis vectors 

U transformation matrix (3Nx3N) 

U transformation matrix (3Nx3N) 

sets up matrix of basis vectors, where 
each row corresponds to a different 
vector. These vectors correspond to 
the linear combinations of principal 
cartesian coordinates which reduce as 
the irreducible representations of the 
point group of the molecule. Each set 
of three columns represent the triplet 
X, Y, Z for each of the N atoms in the 
molecule. The column coefficients for 
each row represent the amount of each 
coordinate in that vector. The first 
six vectors are set equal to zero since 
the translations and rotation are ex- 
plicitly defined and calculated by the 
program. There are 3N remaining vectors 
to be defined, and the basis matrix is 
therefore ((3N+6) x 3N). These vectors 
are then orthogonalized to the 
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translations and rotations, with the six 
redundant coordinates becoming zero. 
These orthogonal vectors are then normal- 
ized, giving the orthonormal U matrix. 
The first six rows of this matrix are the 
explicit definitions of the translations 
and rotations as defined by the Eckart 
conditions, and the remaining 3N-6 ortho- 
normal rows represent pure vibrational 
coordinates. 


PUNCH 

Disk input - U matrix 

Operator input - masses, symmetry block information, 
routing cards for Gaussian 76 

Terminal output- cartesian coordinates of displaced 
vibrational coordinates 

Disk input - cartesian coordinates of displaced 
vibrational coordinates in format suit- 
able for input into Gaussian 76 (includ- 
ing routing cards). 

Description - calculates $;, then displaces along 
+AS; for vibrational coordinates. Uses 
syths;, then -AS; until all i's in sym- 
metry block are done. Then does off 
diagonal elements © [ (S;+AS;, S.+AS,) 

ij>i st 

in same manner until all complete. .Does 
similar procedure on all symmetry blocks. 


FORCERIT 

Disk input - none 

Operator input - upper right hand corner of force constant 
matrix Fo, in row order. 

Terminal output - none 

Disk output - force constant matrix, Fo (3Nx3N) 

Description - simply takes force constant Matrix Fo, 
beginning at row 7 column 7 (first vi- 
brational coordinate) and writes entire 
3Nx3N matrix to disk. Automatically 
fills constants corresponding to trans- 
iarions and rotations to zero. (isen: 


6 N 
Dy x F,.and £ 2 Bas 
i=l j=6 13 i=l j=7 IJ 


APTWRITE 

Disk input - none 

Operator input - dipole derivatives in S coordinates, Pg 
(3x3N), in row order. 

Terminal output - none 


Disk output - Ps matrix (3x3N) 

Description - takes dipole derivatives (in e) in § co- 
ordinates, Po, and writes them to a disk 
set. 


6. 


240 


FLUXRITE 

Disk input - none 

Operator input - charge derivatives in S coordinates (Nx3N), 
row order. 

Terminal output - none 

Disk output - charge derivatives in S$ coordinates (Nx3N). 

Description - takes charge derivatives and writes them 
to disk. Each row of this matrix corre- 
sponds to one of the N atoms in the mole- 
cule, and each column represents one of 
the 3N S coordinates. Therefore, the 
first six columns are all zero. The 
Operator input begins at atom 1 column 7 
(S7) and is entered in row order. 


FREQ 

Disk input - U (3Nx3N); Fo (3Nx3N) 

Operator input - masses, transformation from principal 
cartesians to original cartesians 

Terminal output- frequencies (cm-l) and normal coordinate 
transformation in original cartesian 
coordinates, L~1 (3Nx3N) in mass weighted 
principal cartesian coordinates. 

Disk output - normal coordinate transformation Matrix, 

, 2-1 (3Nx3N) in mass weighted principal 

cartesian coordinates 

Description - using U, Fg are transformed to K, where K 
is the force constrant matrix in mass 
weighted cartesians. This matrix is then 
diagonalized, giving the frequency param- 
eters and normal coordinate transformation. 
The frequency parameters are converted to 
wavenumbers, and the normal coordinates 
are expressed in the original system of 
cartesian coordinates. 


INTENSIT 

Disk input - U (3Nx3N), Po (3x3N), 271 (3Nx3N) 
Operator input - none 

Terminal output- intensities of normal coordinates in 


km/mole. 
Disk output - none 
Description - transforms Ps to Py, in principal carte- 


sian coordinates by using U, then calcu- 
lates Po by multiplying Py and 2-1, 
Finally, Po are converted to intensities 
by using the appropriate constants. The 
numbering of the coordinates corresponds 
to the numbering of the coordinates 

from FREQ. 
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ATOMPOLT 

Disk input - U (3Nx3N), Pg (3x3N) 

Operator input - masses and transformation matrix from prin- 
cipal cartesians to original coordinates 

Terminal output- Atomic Polar Tensors in original coordin- 
ates with units of e. 

Disk output - APTs in original coordinates 

Description - transforms P, to tensors in mass weighted 
principal cartesians by using U. Then 
uses the masses and transformation matrix 
to transform these tensors into APTs inthe 
original cartesian coordinate system. 
Note that the APTs in any coordinate system 
can be calculated by simply entering the 
appropriate transformation matrix. 


FLUX 
Disk input - charge derivatives in S coordinates (Nx3N) , 
U (3Nx3N), principal cartesian coordinates 

Operator input- masses and transformation matrix from 

principal cartesian coordinates to 
_ Original coordinates 

Terminal output-atomic charge derivatives in original car- 

tesian coordinates and charge flux tensors 
. in original cartesian coordinates 

Disk output - none 

Description - transforms atomic charge deriatives from $ 
coordinates to mass weighted principal 
cartesian coordinates by using U. Then, 
using the transformation matrix and the 
masses, calculates atomic charge deriva- 
tives in original cartesian coordinates. 
The charge flux tensors are finally calcu- 
lated by multiplying the charge derivatives 
by the appropriate cartesian coordinates. 
These tensors can be calculated in any co- 
ordinate system by entering the appropriate 
transformation matrix from principal carte- 
sians to the coordinate system desired. 


PXPQ 

Disk input - Py (3x3N), 27) (3Nx3N), U (3Nx3N), Pg (3x3N) 

Operator input- transformation matrix from principal car- 
tesians to original coordinates, adjust- 
ment of APT elements 

Terminal output-intensities using adjuted APTs in km/mole 

Disk output - none 

Description - uses Uand P, to calculate translational and 
rotational Corrections to intensities. 
Transforms 2-1 from principal cartesian 
coordinates to original coordinates using 
the transformation matrix. Calculates in- 
tensities using Py and transformed g-i, 


zjle 


2, 
es 
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L2. 


bie 


14. 


242 


subtracting translational and rotational 
corrections so that the intensities corre- 
sponding to these motions are zero. Allows 
the operator to adjust any element of any 
APT in the molecule and then recalculate 
the intensities. This ma be done as many 
times as is desired. 


BMAT 

Disk input - principal cartesian coordinates 

Operator input- internal coordinate definitions 

Terminal output-internal coordinate definitions, B matrix 

Disk output - B matrix 

Description - the little S vector techniques are used to 
calculate the B matrix elements for the 
internal coordinates; bond stretching, angle 
bending, wagging, torsional bending and 
linear bending. The operator defines the 
number of the coordinate, the type, and 
the atoms involved in the coordinate. A 
weighting factor may be applied, but if 
none is specified, 1.0 is assumed. If two 
coordinates are given the same number, 
then that coordinate is the sum of the two 
coordinates. A coordinate number of @ sig- 
nifies that no more coordindtes are to be 
defined. 


AMATRIX 

Disk input - B matrix 

Operator input - masses 

Terminal output - A matrix 

Disk output - A matrix 

Description - using the technique described in Chapter 2, 
a nonredundant set of internal coordinates 
described by the B matrix is used to calcu- 
late the G matrix. G is then inverted and 
used to calculate the A matrix such that 
AB=E and BA = E., 


FORCE 

Disk input - Fo (3Nx3N0, U (3Nx3N), A 

Operator input - masses 

Terminal output - force constants in internal coordinates 

Disk output - force constants in internal coordinates 

Description - Fg are transformed to K (force constants 
in mass weighted principal cartesians) with 
U, then further transformed to force con- 
stants in principal cartesians with the 
masses. Finally, A is used to transform 
into force constants in internal coordin- 
ates (Fj). 


L3. 


16. 


Lie 


243 


PERTURB 

Disk input - Fy, B matrix 

Operator input - masses, adjustment of force constants 

Terminal output - adjusted force constants, and result- 

ing frequencies 

Disk output - 271 matrix for adjusted frequencies 

Description - allows the operator to adjust any number 
of any of the internal force constants, 
then uses the masses and B matrix to 
transform these adjusted constants to 
mass weighted principal cartesian force 
constants. This matrix is then diagonal- 
ized to give the frequencies calculated 
with these adjusted constants. These 


adjustments can be made as often as desired. 


Finally, the final calculated 2-1 matrix 
corresponding to the final calculated fre- 


quencies may be written to disk. A "return" 


signifies that no more force constants are 
to be adjusted. 


VIBCORDS 

Disk input - principal cartesians, g-+ matrix (3Nx3N) 

Operator input - masses, displacement increment 

Terminal output - none 

Disk output - cartesian coordinates of displaced normal 

coordinates, Q; + AQ; 

Description - the values of Q are calculated by multiplying 
the principal cartesians, the masses and 
2-1, The Q's are then displaced by +AQj 
where the increment is input by the 
operator. Using these displaced Q's, 
the displaced cartesians corresponding 
to these displaced Q's are calculated. 
This is repeated (both + and - AQ;) for 
each Q. 


VIBRATE (written in BASIC) 
Disk input - principal cartesians, cartesian coordinates 
of displaced normal coordinates 
Operator input - radii of atoms, rotation matrix 
Terminal output- XY projection of molecule, option to 
vibrate molecule along any normal 


coordinate 
Disk output - none 
Description - draws a ball and stick model of the xy 


projection of the molecule with the radii 
input, then draws bonds between the spec- 
ified atoms. Allows rotation of the mol- 
ecule by entering the appropriate 
transformation matrix. Furthermore, the 
picture can be labelled. Finally, the 
molecule can be made to vibrate along any 
of its normal coordinates. This vibration 
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can be frozen in either +AQ; phase, or 
another normal coordinate can be chosen 
to vibrate along by pressing "C". The 
bonds are no longer drawn by entering 
-l, 2, labelling is ended by typing 
“Control E.” 


SUBROUTINES (all contained in library MYLIB. REL) 


1. FILEREAD - reads a file from disk 

2. FILERITE - writes a file to disk 

3. MULTIP - multiplies two matrices 

4. JACOBI - diagonalizes a square symmetrix matrix 

5. SQINUT - inverts a square symmetrix matrix 

6. STRET - calculates B matrix elements for stretching 
coordinate 

7. BEND - calculates B matrices elements for bending 
coordinate 

8. OPLA - calculates B matrix elements for wagging coordinate 

9. TORS - calculates B matrix elements for torsional bending 
coordinate 

10. LIBE - calculates B matrix elements for linear bending 
coordinates 

ll. VECTOR - calculates the vector between two atoms given 

cartesian coordinates : 

12. DOTP - calculates the dot product of two vectors 

13. CROSS - calculates the cross product of two vectors 

14. SINE - calculates the sine of an angle from the cos of 
an angle 


APPENDIX B 
NORMAL COORDINATE TRANSFORMATION MATRICES 


The normal coordinate transformation matrices used in 
these calculations have all been stored on 8 inch floppy 
disks. These disks are formatted single side double den- 
sity for use with the CP/M operating system, version 2.2 
(73). These data sets were created with Microsoft FORTRAN- 
80 (74) source code. These matrices have been stored as 
vectors, in row order, one element at a time, and should be 
read into memory in the same fashion. This appendix con- 
tains a list of filenames corresponding to the various 
normal coordinate transformation matrices used in these 


calculations, for all of the systems studied here. 


1. Theoretical Water Monomer (_MONOT) (N = number of atoms) 

Principal cartesian coordinates XMONOT (3xN) 

U matrix UMONOT (3Nx3N) 

Force constants (K) S coordinates FMONOT (3Nx3N) 

Dipole derivatives (Pg ) S coordinates DMONOT (3x3N) 

Normal coordinate transformation (2-1) principal car- 
tesians LMONOT (3Nx3N) 

Atomic Polar Tensors (Py) original cartesians PMONOT 
(3x3N) 

Atomic charge derivatives for calculating charge flux 
tensors, S coordinates CMONOT (Nx3N) 

B matrix (B) principal cartesians BMONOT (# internal 
coordinates x 3N) 

A matrix (A) principal cartesians AMONOT (3N x # internal 
coordinates) 

Internal force constants (F;) IMONOT (# internal coordin- 
ates x # internal coordinates) 

Coordinates for picturing normal vibrations NMONOT (dim- 
ensions not needed) 
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The filenaming convention employed here is the same in every 
case. The system being calculated is given a code name (up 
to seven letters and/or numbers) and then the appropriate 
Matrices are defined by the one letter prefix, as given 
above. In some instances (i.e., in the case that redundan- 
cies exist in internal coordinates), two additional matrices 
may have to be defined. These matrices are the transformation 
matrix from the redundant set of internal coordinates to the 
nonredundant set of symmetry coordinates, and then the symme- 
trized B matrix. These matrices are given the prefix codes 
Q and S respectively. 
The code names for the systems studied here are listed 

below. 
2.832 A° (_LINDIMT) 

-982 A° (_LINDIM1) 

-682 A° ( LINDIM2) 


2 

2 —_ 

3.132 A° (_LINDIM3) 
2.532 A° (_LINDIM4) 


2. Linear Dimer 


DAWA DAD 
nuuw we 


«74 A*~( CY@DIM1) 
-65 A° ( CYCDIM2) 
-99 A° ( CYCDIM3) 
.74 A° No inversion center (_CYCDIM4) 


3. Cyclic Dimer 


ADAD 
iouu i 


4. Bifurcated Dimer R = 2.90 A® (_BIFDIM1) 
R = 3.05 A®° (_ BIFDIM2) 
R = 2.90 A° perpendicular molecules 
(_ BIFDIM3) 
+ 
5. H,0” (_H20P) 
6. O+H,0-0 (_H203N) 
7. o-+-H,0-07! (4203-2) 


2 


8. 07%-H,0-072 (_H203-4) 


APPENDIX C 
REDUNDANT INTERNAL COORDINATES 

In the case where more internal coordinates than vibra- 
tional degrees of freedom are defined, a redundancy condition 
exists. This means that the set of internal coordinates are 
dependent, and the redundancy is the linear dependency be- 
tween some of the coordinates chosen. The function which 
describes the redundancy condition does not change through- 
out vibrational motion and therefore makes no contribution 
to the solutions of the secular equation. If one chooses to 
simply eliminate the dependent coordinates then it is often 
the case that the doventiel function cannot be properly de- 
scribed in terms of symmetry coordinates. Therefore one 
usually constructs as many symmetry coordinates as there 
are internal coordinates where one (or more) of these sym- 
metry coordinates is exactly the redundancy condition and 
can be easily removed. However, the set of appropriate 
symmetry coordinates is not unique, and depends in part on 
both the geometry and the original internal coordinates 
defined. In addition, it can be shown that while the poten- 
tial function in symmetry coordinates uniquely defines the 
vibrational properties of the molecule, the potential func- 
tion in the original set of internal coordinates is not 
uniquely defined. 
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If R is the set of redundant interanl coordinates (where 
in this case we use an example with one redundancy), then the 


transformation to symmetry coordinates S is written 
S(N) = U(N x (N+1)) R((N+1) x 1) (99) 


where U is a transformation matrix, R contains one redundancy 
condition and is of dimension (N+l), S is of dimension (N) 
because the redundancy has been explicitly defined and re- 
moved, and U is of dimensions (N x (N+l)). If we assume 

that U represents a unitary transformation matrix, then the 


following are true: 


U(N x (N+1)) U*((N+1) x N) = E(NXN) (100a) 


U" ( (N41) XN) U(N x (N+1)) = E((N+1) x (N+1)) (100b) 


Using this assumption, equation 99 may be pre-multiplied by 


+ ‘ 
U to give 


R=u'S (101 
or by substituting with equation 99 


R= (U U)R (102) 


The set of internal coordinates R were originally defined to 
contain a redundancy condition. Equation 102, in conjunction 
with equation 100b directly contradicts this statement. If 
equation 100b is correct, then equation 102 indicates that 
there is a one-to-one correspondence in the internal coordin- 
ate space, or that the internal coordinates are linearly 
independent. However, the presence of a redundancy condition 


necessitates that the internal coordinates by linearly 
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dependent. Thereefore, we May conclude that the assumption 
that U is a unitary transformation matrix is not correct, 


and subsequently, equation 100b is not true, 
U" ( (N+1) x N) U(N x (N+1)) # E( (N41) x (N+1) ) (103) 


Therefore we define 


C((N+1) x (N+1)) = uty (104) 
Using equation 100, we note write 
ce= uw y= uU=¢ (105) 
c? = € idr ese = OE((N+1) x (N+1) (106) 
C(C-E) = 0 (107) 
where if we premultiply equation 107 by the inverse of c 
we obtain 
(C-E) = 0 or C=E (108) 


but equation 108 directly contradicts equation 103, so we must 


conclude that C has no inverse, or in other words 
u| = Ic] = 0 (109) 
The potential energy of this system can be expressed as 
—+_+ — a 
2v=S FS = R°UF_LUR=R FF (110) 
where we have made use of equation 99 and where 


P, =U" reg (111) 


1 
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Alternately, we can write the potential energy as 


2v=R'F,R=S‘uF.U S=S'"FS (112) 


where we make use of equation 101. Therefore, we obtain 


F. = UF.U (113) 


By substituting equation 113 into equation 111 and using equa- 


tion 104 we obtain 


F.O'U jor F.. = CF;C (114) 


Equation 114 points out the fact that Fs cannot be unique, 


but rather the transformation given by equation 114 can be 


successively applied an infinite number of times, and each 
' 
PF; is different, but correctly determines F via equation 113. 


Therefore, we have shown that while the values of F are 


determinable and can be used to solve the vibrational secular 
equation, the transformation of these parameters back to Fs 


in the original redundant coordinates is not uniquely 


determined. 
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